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Lt. General Bernard A. Schriever 


Bernard A. Schriever was born in Bremen, Ger- 
many, on September 14, 1910. He came to the United 
States with his parents in 1917 and was naturalized 
in 1923. He received his early schooling in San 
Antonio, Texas, and was graduated from Texas 
A & M, College Station, Texas, in 1931, with a B.S. 
degree in engineering. He received the M.S. degree 
in aeronautical engineering in 1942 from Stanford 
University, Stanford, Calif. 

During World War II, he organized material and 
maintenance supply lines for the 5th Air Force Serv- 
ice Command, and flew 63 combat missions in the 
South Pacific. At the war’s end, he was Commander, 
Advanced Headquarters, Far East Air Service Com- 
mand. 

In 1946, he was assigned to Air Force Headquar- 
ters as Chief, Scientific Liaison Section, Deputy 
Chief of Staff, Material. In 1949, he entered the Na- 
tional War College. Following his graduation, he re- 
turned to Headquarters, USAF, for assignment as 
Assistant DCSD for Development Planning. 

He became Assistant to the Commander, Air Re- 
search and Development Command, in May, 1954. In 
this position, he organized the Ballistic Missile Divi- 
sion (then Western Development Division) and was 


its first commander. This organizational arrangement 
provided a clear line of authority to top echelons in 
the Air Force and the Department of Defense. 

In August of that year, he assumed the additional 
responsibilities of Commander, Air Force Ballistic 
Missile Division, ARDC. In this position, he imple- 


mented the unique management concepts by which. 


the ARDC/AMC/SAC team develops the Atlas, 
Titan, Thor, and Minuteman ballistic missile weapon 
systems, along with their complementary space sys- 
tems and early warning and communications satel- 
lites, and brings them to initial operational status. 
One of the most far-reaching of these management 
concepts is that of concurrency, wherein all phases of 
the normal design-to-deployment cycle are pursued 
simultaneously rather than sequentially. Concurrency 
was an important factor in enabling the scientific- 
industrial-military team to cut roughly in half the 
normal weapon system development cycle. 

In April, 1959, General Schriever was appointed 
Commander, ARDC, with the rank of Lieutenant 
General. ARDC includes eleven major research, de- 
velopment, and test centers, and two major develop- 
ment branches: the Ballistic Missile Division and the 
Air Defense Systems Integration Division. 
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Guest Editorial 


OR many years the Air Force has manifested 
an intense interest in science and technology. 
Because of this interest, it has initiated or sup- 
ported a large number of programs which have con- 
tributed in a major way to aeronautics and space. Ex- 
amples of some positive results are the development of 
the first large liquid rocket engines in this country ; the 
development of the precision-floated gyro capability 
at the Instrumentation Laboratory, Massachusetts In- 
stitute of Technology ; and numerous contributions to 
aviation medicine, to the basic knowledge of stability 
and control, and to systems engineering. Major im- 
provement of airborne computers by the Air Force 
has provided a firm basis for later missile and space 
applications. Light-weight components developed and 
tested under Air Force sponsorship have likewise 
found many applications. Our stellar inertial and ter- 
minal guidance techniques are likely to find applica- 
tion in space technology as they already have in mis- 
siles. 
In the past five years, the state-of-the-art has re- 


‘ceived an even greater impetus because of the accel- 


erated ballistic missile programs. These programs 
have benefited from the assistance of some of the 
best scientific, technical, and managerial talent in the 
country. The keynote of the entire effort has been 
teamwork and cooperation among science, industry, 
and the military. The result has been a compression 
of the time between initial planning and operational 
availability of important weapon systems, with nu- 
merous corollary results such as major additions to 
production and test facilities, and to the experience of 
personnel. 

Air Force space programs such as Discoverer, 
Atlas in orbit, and the Thor-Able series came as a 
natural consequence of the vehicles, facilities and ex- 


perience developed in our ballistic missile programs. 
The result of these programs has been some out- 
standing successes in space exploration, with a broad 
capability for further important work in space by the 
Air Force and by other agencies. A little thought is 
sufficient to convince one that certain space applica- 
tions are a necessary adjunct to strategic weapon sys- 
tems such as ballistic missiles. It is therefore our in- 
tention to continue to pursue new programs with 
great vigor. Without doubt, a sizeable portion of our 
research and development budget will be devoted to, 
advanced systems and techniques related to space ap- 
plications. It is obvious that a substantial part of this 
will be in electronics, guidance, communications, and 
related areas. 

The advances in technology which I have men- 
tioned have resulted in a number of significant tech- 
nical papers, some of which have already been pub- 
lished. Examples of the subjects are tropospheric 
propagation, accelerometers, gyroscopes, testing, long 
and short baseline radio-guidance systems, inertial 
guidance, computer circuitry, and reliability. The 
paper on “Baseline Guidance Systems,” by Grisetti 
and Mullen, published in the December, 1958 issue of 
these TRANSACTIONS, was based upon work per- 
formed in one of our programs. Scientists and engi- 
neers working with us have been encouraged to 
publish the results of their work and to participate in 
lectures and in symposia. The accompanying papers 
are exemplary of some recent work by scientists and 
engineers on our teams. Publication of such articles 
is consistent with my belief that the maximum dis- 
semination of unclassified information is important 
and is very much in the national interest. 

B. A. SCHRIEVER 
Lieutenant General, USAF 
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Scientific Objectives of the Able-3 Program* 


APPLIED PHYSICS DEPARTMENT, SPACE TECHNOLOGY LABORATORIES, INC. 


Summary—The purpose of the Able-3 program is to place a 
scientific observatory in a highly elliptic satellite orbit. A scintil- 
lation counter, ion chamber, Geiger counter, and proportional 
counter telescope in the payload permit comprehensive mapping 
of corpuscular radiation over a large volume of the geomagnetic 
region, including the Van Allen radiation belts. Simultaneously, 
the fundamental static quantities, electron density, and vector 
magnetic intensity are to be measured. A VLF experiment is 
included to measure upper atmosphere VLF noise and to study 
the propagation characteristic of VLF between the ground and 
the satellite. A radio scintillation experiment using radio tele- 
scopes will measure the effects of ionospheric irregularities on 
amplitude and phase. Micrometeorite detection equipment in the 
payload provides a two-point momentum spectrum of micro- 
meteorites encountered, and an image scanner is included to 
examine reflection characteristics of the earth. 


I. INTRODUCTION 


HE goals of the Able-3 program in placing a scientific 
observatory in a large elliptic satellite orbit (nominal 
apogee 20,000 miles, perigee 150 miles) are to pro- 
vide scientific data concerning the earth and its adjacent 
space and to test the equipment for application to later 
space flights. It is the purpose of this paper to establish the 
framework of scientific background and specific objectives 
and instrumentation for the Able-3 scientific measurements. 
The most important general knowledge which has 
emerged during the past decade or so, and which is strongly 
supported by the IGY findings, is the great degree of 
intimacy between solar and ionospheric weather. It is no 
longer adequate to envisage the sun as a remote glowing 
sphere radiating into space at 6000° K. Actually, the coro- 
nal plasma (at a temperature of order one million degrees) 
spreads out over all of interplanetary space. The density of 
ionized gas in interplanetary space is everywhere greater 
than previously supposed. Photographic studies of the mo- 
tion of comet tails indicate a continuous outpouring of 
coronal gas with velocities of order 500 to 1500 km/second. 
From these observations, it has been possible to estimate 
a plasma density at the orbit of earth of 500 ions per cubic 
centimeter, and it is expected that much higher densities 
occur during violent solar outbursts. Theories of whis- 
tling atmospherics also require densities of this same mag- 
nitude. Considerable additional evidence also supports the 
same conclusion. 
This highly conductive wind can be expected to draw 
out the solar magnetic field radially. However, this wind 


* Manuscript received by the PGMIL, July 21, 1959. The Able-3 
vehicle was launched on August 7, and has since been designated 
Explorer VI. 

+ Los Angeles, Calif. 


cannot penetrate freely across the geomagnetic field. Cal- 
culations of the interaction energies indicate that an every- 
day solar wind may penetrate to about 25,000 km above 
the surface of the earth, while a strong solar wind is ex- 
pected to penetrate to about 3500 km. The motion of this 
tenuous, ionized gas is controlled by electrical and magnetic 
forces which are themselves strongly perturbed by the 
motion of the gas itself. We can therefore expect this com- 
plicated electrodynamical interplay to define a magnetically 
erratic region which may properly be deemed the end of the 
earth’s magnetic field. This region is of paramount impor- 
tance to study experimentally since it is probably the 
source of those systems of currents and associated geo- 
magnetic disturbances that bear a close relationship to the 
morphology of the aurora and, indeed, the entire atmos- 
phere. Although the complicated cooperative plasma dy- 
namics of this zone of interaction are yet to be worked 
out, some theoreticians appear to believe that conditions 
are favorable for accelerating charges to the high energies, 
100 kev, observed in auroral arcs. One viewpoint concern- 
ing the recently discovered, trapped, high-energy particles, 
referred to as the Van Allen radiation belt, is that they 
too are the result of direct impact of the solar wind. To 
be sure, it has already been surmised that the auroral phe- 
nomenon may be a direct manifestation of instabilities in 
the radiation belt. However, the extent to which the aurora 
is under indirect solar control or is a manifestation of ter- 
restrial forces (radiation belt) is not as yet ascertained. 
Moreover, no explanation has as yet been offered as to why 
it is a night-side phenomenon. 

It should be noted that an alternative point of view for 
the origin of the inner radiation belt exists, that is, the 
neutron decay hypothesis. This theory assumes that neu- 
trons, back-scattered by nuclear collision produced by 
cosmic ray impact with molecules in the atmosphere, decay 
into protons and electrons, some of which become mag- 
netically trapped. Such an hypothesis does not, of course, 
require the interaction between the corona wind and the 
geomagnetic field. An understanding of how such magneto- 
dynamic processes store energy and produce high-energy 
corpuscles in the solar-terrestrial realm may go a long way 
toward explaining the still unanswered questions concern- 
ing the primary mechanisms which produce and accelerate 
cosmic ray particles. 

In addition to the corpuscular phenomenon, it is sup- 
posed that the solar-earth impact zone is a generator of a 
variety of electromagnetic, very low frequency CVEEED) 
signals which have been observed. This noise propagates to 
earth approximately along magnetic field lines in the man- 
ner of lightning-generated whistlers. This VLF emission 


130 IRE TRANSACTIONS ON 
appears to be strongly correlated with solar activity. For 
many, but not all, of these signals a recent theory suggests 
that the excitation mechanism is similar to the operation of 
a traveling-wave tube. That is, the electromagnetic noise 
is excited in the exosphere by streams of high-speed parti- 
cles traveling along the magnetic field lines. This electro- 
magnetic emission is a resonance effect which occurs when 
the velocity of the ionized stream is equal to the phase 
velocity of the electromagnetic wave. A model for most of 
the exosphere can be constructed in terms of this theory by 
detailed interpretation of several categories of received 
noises. This model requires that the ratio of electron den- 
sity to magnetic field strength be nearly constant. Unfortu- 
nately, this theory, which is related to the ordinary Cerenkov 
theory, does not appear advanced enough to predict signal 
intensities which can be tested against observed signal 
strength. Moreover, there appear to be several categories of 
VLF emission which cannot be explained in terms of the 
traveling-wave tube mechanism. The Cerenkov mechanism 
is a distinct possibility as a source of VLF radiation. There 
is a large range of low frequencies for which the phase 
velocity of the extraordinary propagation mode is less 
than the speed of light by as much as a factor of one 
hundred. According to the Cerenkov theory, electrons, pro- 
tons or any charged particle moving through the static 
plasma with velocities greater than the electromagnetic 
phase velocity ought to radiate, and the total energy of this 
radiation is propottional to the distance traveled by the 
charged stream as well as its density. It is the long mean 
free path for electrons and protons (approximately 107° 
cm) which provides the opportunity for accumulating sub- 
stantial signal strength. Moreover, unlike the applications 
in nuclear physics, Cerenkov emission is possible in the 
exosphere for nonrelativistic as well as relativistic par- 
ticles, a range of velocities likely to be present in the Van 
Allen radiation belt. Reasonable estimates indicate that 
streams of high-energy electrons, or protons, could provide 
measurable VLF signals if the electrons were to contribute 
incoherently to the Cerenkov emission. More realistic 
estimates accounting for the very long wavelengths in- 
dicate that aggregates of electrons when bunched into sizes 
considerably smaller than the wavelength may contribute 
coherently rather than independently with the result that 
relatively large VLF signal intensities are expected. How- 
ever, if the bunching in the high-energy streams extends 
over regions comparable or larger in size than the wave- 
length, no emission may occur in the sense that a steady 
current can produce no time-dependent polarization of the 
ambient medium. Consequently, the VLF signature bears 
a close relationship to the sizes of the electron aggregations 
in the high-energy streams and such bunching is clearly 
related to the injection mechanism by which such streams 
are created in the upper atmosphere. 

Additional sources of VLF background are “brem- 
strahlung,” the radiation produced by high-energy elec- 
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tron scatter from protons and the cyclotron radiation from 
the trapped Van Allen electrons. Estimates indicate that 
the latter mechanism is the more appreciable and that, now, 
coherent radiation of considerable intensity arises from 
bunches which are small in extent compared with the radii 
of their orbits rather than radiation wavelength while no 
radiation is expected for bunching sizes that are compar- 
able to orbital radii. 

It also appears that the signal strength should be con- 
siderably greater for ground reception, especially in polar 
areas where the strong geometric convergence factor effec- 
tively collects the radiation originating in a large area. 
This follows from the conservation of magnetic flux. 
It follows that the signal amplification on the ground 
would increase in proportion to the cube of the distance 
between the center of the earth and the excitation region. 
Conversely, it may be expected that VLF emission which is 
easily detected on the ground need not be as easily detected 
in a satellite. These geometric factors concerning signal 
amplification do not depend upon the particular VLF emis- 
sion mechanism. 

Again, regardless of the specific mechanism responsible 
for the VLF emission, there is the intriguing possibility 
that a trapping process exists which tends to contain and 
magnify the VLF energy in the exosphere. The contain- 
ment process envisaged here is simply one of selective re- 
flection, or partial reflection, produced by the gradients in 
the refractive index of the propagating mode, gradients 
which can look fairly abrupt in terms of the scale of wave- 
lengths characterizing VLF radiation. Such “discontinui- 
ties’ may occur upon approach toward the ionosphere 
from external regions. The increase in electron density as 
well as a corresponding increase in magnetic intensity 
produce the gradient in refractive index. It is known, for 
example, that whistlers do not always appear when there 
is a discharge of lightning. Thus, conditions for propaga- 
tion appear to vary in a manner which is presumably un- 
known, and present knowledge cannot rule out the con- 
sequent possibility of the existence of a VLF “corona” on 
earth and perhaps upon other magnetized planets. The im- 
portance of investigating this new propagating medium in 
the upper atmosphere and, in particular, of identifying 
the emission process or processes cannot be overempha- 
sized. If, for example, the VLF phenomenon were under- 
stood in terrestrial space, then VLF reception might prove 
valuable as an environmental probe of other planetary 
bodies. Apart from general scientific interest, there is the 
obvious possibility for a new means of long-range com- 
munication with low-power transmitters over natural, 
giant, magnetic propagation ducts. 

As is usual when great gaps in understanding exist, the 
discussion is bold and inconclusive. Nevertheless, it 
furnishes valuable guidelines for an experimental space 
program. It indicates the necessity for a coherent meas- 
urement program throughout a relatively enormous region 
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of space extending over perhaps 1 to 7 earth radii. This 
region includes the solar wind-geomagnetic interaction 
zone, which is inaccessible to conventional sounding tech- 
niques except by man-made VLF signals. This region also 
contains the provocative Van Allen radiation belt. At the 
present time very little detail is known about the energy 
spectrum, composition, distribution in space, and temporal 
behavior of these “‘belts.”” Indeed, the most valuable and 
interesting aspects of the radiation belts will be their spatial 
symmetries or asymmetries and their fluctuations as a 
function of time. 

At this point it is appropriate to enumerate a few of 
the principal lines of investigation of the Able-3 earth 
satellite. More specific description of the methods of in- 
vestigation appear in later sections. 


1) The satellite will be equipped with high and low 
energy cosmic ray detectors which will be operated 
continuously inside and outside the radiation belt. 
The instruments are a scintillation counter, Geiger 
counter, ionization chamber, and a_ proportional 
counter telescope. 

2) The vehicle will be equipped to receive a VLF keyed 
signal at 15,500+50 cps from a ground-based 
transmitter, and also measure the emission noise in- 
tensity in this same band. For the transmitted signal, 
the pulse travel time as well as the intensity will be 
measured. 

3) The third and perhaps central set of measurements 
is of the two vital, dynamically associated “ambient” 
space parameters, electron conceitration and mag- 
netic vector. These will be measured continuously 
throughout the interaction zone identified above. The 
basic instrument to be used for these studies are a 
Doppler-difference frequency technique for electron 
density measurement, and a dual magnetometer sys- 
tem, employing a search coil magnetometer and a 
flux gate magnetometer which will determine the 
magnetic field vector. 


It is to be noted that these experiments will be carried 
out over a highly elliptical orbit, for a nominal lifetime of 
a year, with an apogee in excess of 20,000 miles from the 
surface and with a perigee of 150 to 300 miles. Since the 
orbit is initially inclined at 45°, its precessional motion 
will permit measurements to be made over a more than 
adequate volume of the space about the earth. 

Since the upper regions of space which we have been 
discussing are geomagnetically coupled with the ionosphere, 
it is essential to the understanding of either region that we 
attempt to correlate some of the effects of their interaction. 
To identify the pertinent areas of study, it is first necessary 
to distinguish the two regions, both of which are, of course, 
ionized and magnetic. There is no sharp demarcation be- 
tween them, but the increase in neutral particle density in 
descent ultimately will define layers, the motion of which 
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is no longer controlled by magnetohydrodynamic forces 
owing to the increased electronic collision frequency. Con- 
sequently, an interim region exists where the macroscopic 
plasma motion becomes gradually disentangled from the 
geomagnetic field. The rate at which this occurs is, of 
course, influenced by solar photo effects as well as electron 
recombination processes. As a result, a variety of distrib- 
uted ionization irregularities as well as anisotropies as- 
sociated with the geomagnetic field occur throughout this 
region, including a large but unknown part of the iono- 
sphere. However, the disturbances through this important 
region can be revealed in depth by radio studies which 
utilize wavelengths that are small relative to the scale sizes 
of the irregularities. This affords a very important ad- 
vantage not possible with VLF measurements in such 
media, namely that the theory of radio propagation is well 
established for such media. Thus, the natural tool for re- 
lating the ionosphere with the higher realms already under 
investigation in the Able-3 program is electromagnetic 
radiation at frequencies which will register the integrated 
coupling effects. To accomplish this merely requires that 
the satellite be tracked at those frequencies in an appro- 
priate manner. 

It is now desirable to convey a brief description of the 
motion of the ionosphere. The “undisturbed” ionosphere 
exhibits a wide variety of motions ranging from large- 
scale tides, to east-west drifts, to small-scale irregularities. 
No composite account of these motions can as yet be given, 
notwithstanding the existence of a great amount of per- 
tinent but isolated data accumulated by radio observations 
of meteor trail configuration, auroral forms, and radio- 
star scintillation rate measurements. The latter have pro- 
vided useful data on line-of-sight irregularities as well as 
information concerning the drift motion of these irregu- 
larities, sometimes reffered to as blobs. It has also revealed 
puzzling phenomena such as the night-side maximum in 
scintillation amplitude, and the night-side phenomenon of 
spread F. However, these manifestations may only be 
superficial consequences of larger geomagnetic systematics 
perhaps involving the highly energized radiation belt. As 
previously mentioned, the aurora is also a “puzzling” 
night-side phenomenon. 

To gather useful data concerning these complex motions 
in the ionosphere, it is therefore important to include line- 
of-sight scintillation rate measurements for studies of ion 
cloud drifts and ionospheric irregularities simultaneously 
with the measurements made in the upper atmosphere. 
This is easily accomplished as part of the Able-3 program 
merely by tracking the telemetry beacon already on the 
vehicle. 

It is perhaps appropriate to elaborate further on the im- 
portance of measuring ionospheric fluctuation rates. The 
phase fluctuation imposed upon signals from distance 
sources produces frequency dependent angle-of-arrival 
errors, or source flicker. This introduces a natural limit to 
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tracking accuracy. Therefore, a thoroughgoing knowledge 
concerning expected phase and amplitude signal fluctua- 
tions for arbitrary line of sight is essential for the design 
and evaluation of long-range missiles and space vehicles. 
At the present time the data which are available are those 
accumulated by radio-star scintillation measurements. The 
data are inadequate in two important respects: 


1) The statistics, although probably adequate for deter- 
mining the variance of signal amplitude and direc- 
tion of arrival, are inadequate for establishing con- 
fidence in the power spectrum. Since star signals 
are weak, the poor signal-to-noise ratio causes poor 
time resolution of the signals and this interferes with 
extracting an adequate power spectrum from the 
records. This problem would not arise in the case of 
the satellite source. 

2) Little information is available concerning the effect 
of varying the distance between two antennas form- 
ing an interferometer. Such information is important 
for determining the scale, or size, of the irregularities 
of electron density in the ionosphere. The degree of 
anisotropy in this quantity is also important. For 
example, if it is elongated in the direction of the 
geomagnetic field, the region which would produce 
the major scintillation effect must be high enough so 
that collisions no longer dominate the motion of the 
ionized clouds, but rather magnetoelectrodynamic 
forces are determinate. The scale of electron density 
irregularities and the degree of isotropy can be easily 
measured by the addition of a mobile third receiver 
to an interferometer system, and such a unit has been 
included in the Able-3 research program. In addi- 
tion, such an arrangement will allow the determina- 
tion of drift motions in all directions. 


The instrumentation for the scintillation rate measure- 
ment involves two 40-foot diameter radio telescopes on an 
east-west baseline and a portable Yagi antenna. The Able-3 
108-me beacon signals will be tracked while scintillations 
in phase and amplitude as a function of baseline orienta- 
tion and antenna separation will be recorded. 

Before discussing the remaining scientific experiments 
of the Able-3 program, it is important to underline the 
significance of those which have already been discussed. 
Clearly this group of experiments involves measurements 
of quantities that, according to available facts as well as 
theoretical inference, are dynamically interrelated rather 
than isolated. The basic data gathered from such a co- 
herent set of experiments can be expected to be invaluable, 
especially when interpreted together with world-wide 
monitoring measurements that are continuously and in- 
dependently being carried out. Although it is impossible to 
predict the exact use of these comprehensive measurements 
in terms of what new theoretical or practical advances can 
ensue, it is safe to say that this new and interrelated in- 
formation will at least furnish a far more precise descrip- 
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tion of the important phenomena in the upper atmosphere 
and eliminate erroneous views. This aggregate of measure- 
ments certainly achieves the most coordinated simultaneity 
of observations in extraterrestrial space that has ever been 
attempted. 

Space research in the future will undoubtedly require 
instruments capable of providing pictures of celestial 
objects. The measurement of radiation ranging from the 
ultraviolet through infrared wavelengths can provide im- 
portant scientific data. It is particularly significant for 
solar system probes to the planets (or moon) because the 
data may be obtained from early missions, ones in which 
the vehicles do not make a soft landing. 

With regard to this photo range, the planetary bodies 
may be divided into three groups: first, bodies with little or 
no atmosphere, for example, Mercury or the moon; sec- 
ond, bodies with an atmosphere which is at least partially 
transparent, for example, earth and Mars; third, bodies 
with an atmosphere which is opaque, for example, Venus. 

Next, consider the measurements in more detail. Taking 
the moon as an illustration of the first class of bodies, one 
of the unsolved, but extremely important, problems is the 
determination of the structure and substance of the sur- 
face in local areas; for example, are such areas covered 
with dust or hard rock? A brightness scan of the surface 
giving reflection characteristics as a function of angle 
would help resolve this question. In addition, information 
on the appearance of the backside of the moon would help 
determine the origin of geophysical features of the lunar 
surface. 

Useful measurements of bodies with partially trans- 
parent atmospheres such as the earth would include meas- 
urement of cloud cover, movement of cloud masses, nar- 
row-band absorption measurements, and reflection char- 
acteristics of the surface at various wavelengths. The first 
two would be useful for predicting weather, the third for 
determining temperature structure and density of certain 
molecules in the upper atmosphere, and the first and fourth 
for determination of the albedo and the heat balance of the 
planet. 

For the case of a planet with an opaque atmosphere, 
such as Venus, brightness scans across the top of the at- 
mosphere would yield limb darkening data, which are use- 
ful in determining the temperature structure of the atmos- 
phere. Depending upon the data obtained, it may also be 
possible to ascertain the particle sizes in the upper atmos- 
phere, and also to determine whether or not they are mole- 
cules, vapor droplets, or dust particles. Local variations in 
these data may also help to estimate the amount of atmos- 
pheric circulation. In addition, energy measurements in or 
near a particular absorption band may be used to make 
quantitative estimates of the presence of certain atmos- 
pheric gases. 

In view of the wide variety of future applications for 
photo sensors, such a device has been placed on the Able-3 
vehicle. Because of the very stringent weight and volume 
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requirements, it is a relatively crude sensor. However, it 
does incorporate an automatic scanning system so that data 
will be obtained from the entire surface of the earth. Ex- 
perience gained in the design and use of this system will 
be invaluable for the eventual development of proficiency 
in acquiring and transmitting high-resolution pictorial 
images from space vehicles. 

Since the Able-3 studies are preparatory to future space 
probe programs, it was obvious that, quite apart from the 
considerable scientific interest attached to micrometeoritic 
data and which is discussed below, superior statistics as 
to their momentum flux are required for the practical pur- 
pose of evaluating their damage potential. Consequently, 
a micrometeoritic momentum spectrometer is included as 
an important part of the Able-3 program. 


II. RapIaATION EXPERIMENTS 


A. Introduction 


Three radiation experiments are carried in the Able-3 
payload: a scintillation counter prepared at Space Technol- 
ogy Laboratories, Inc., a Geiger counter and ionization 
chamber combination designed and built at the University 
of Minnesota, and a proportional counter telescope con- 
ceived by scientists of the Enrico Fermi Institute of 
Nuclear Studies of the University of Chicago and de- 
signed and built at the Chicago Midway Laboratories, a 
division of the University’s Laboratory for Applied Sci- 
ence. 

The primary objective of these experiments is to study 
the intensity, composition, energy spectrum, extent, and 
time variation of the terrestrial radiation belts discovered 
as a result of experiments carried on the Explorer satellites 
and Pioneer space probes. Previous experiments have 
shown that there are two distinct regions or belts of radia- 
tion surrounding the earth. The inner belt contains a large 
number of highly energetic particles believed to be protons 
produced by the decay of neutrons created in the earth’s at- 
mosphere by the primary cosmic radiation. In addition, it is 
believed that the inner belt also contains large numbers of 
low-energy particles, probably electrons of a solar origin. 
These low-energy particles may be much more numerous 
than the high-energy ones. The particles in the outer radia- 
tion belt have much lower energies than those in the inner 
belt and are believed to be electrons of solar origin released 
during solar flares. It now appears that a large percentage 
of these electrons have energies well below 500 kev. 

The Able-3 satellite passes through the radiation belt 
twice during each revolution. Due to the apsidal motion of 
the orbit, the vehicle will pass through nearly all of the 
belt between 45°N and 45°S latitude during its nominal 
one-year useful lifetime. Thus it will be possible to map 
nearly the complete radiation belt during this time. How- 
ever, since the apsidal motion is slow, it will also be possi- 
ble to study the time variations of the belt over both long 
and short periods. 
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B. Description of the Experiments 

1) Scintillation Counter (STL): The scintillation coun- 
ter is designed to measure the total flux of electrons of 
energy greater than 225 kev and protons of energy greater 
than 2 mev. The detector is a plastic cylinder 2.5 centi- 
meters in diameter and 0.64 centimeter thick. Light pulses 
are detected with a Dumont 6467 photomultiplier. The 
detector is surrounded by 10 mg/cm? of aluminum. There is 
a 5-centimeter aluminum window in the side of the vehicle, 
also 10 mg/cm? thick. To decrease noise in the output, 
the scaling circuit is biased to accept pulses equivalent to 
an energy loss in the plastic of at least 200 kev. The alumi- 
num window subtends a solid angle of about 0.5 steradian 
at the detector. 

Data are recorded and transmitted in both analog and 
digital form. The output of the photomultiplier is amplified 
and then fed to a scaling circuit with a scaling factor of 
2°. The analog data transmitted consist of the output of 
the 28, 21‘, and 2° scalers. These pulses are transmitted on 
the same channel and are distinguished by their amplitudes, 
which are respectively 12, 23, and 66 per cent of the maxi- 
mum bandwidth. In low flux regions, the 2° scale output 
furnishes most of the usable data. As the counting rate 
increases, these pulses blend into nearly a straight line, 
but the 214 scale pulses are still easily resolvable. In intense 
regions only the 2?° scale pulses are resolvable. Due to the 
approximately one microsecond resolving time of the de- 
tector and its associated circuitry, it is unlikely that the 
counting rates wil become high enough to make it impossi- 
ble to resolve the 2?° scale pulses. The output of the 2° 
scaler is transmitted as digital output. 

2) Ionization Chamber—Geiger Counter (University of 
Minnesota): The University of Minnesota experiment 
uses two detectors, an ionization chamber and a Geiger 
counter. The ion chamber is of the integrating type de- 
veloped by Neher and Millikan of California Institute of 
Technology, and recently used extensively on University 
of Minnesota balloon flights. It consists of a 3-inch sphere 
filled with argon gas to an absolute pressure of approxi- 
mately 6 atmospheres. The ions are collected by a quartz 
rod, the upper end of which is coated with Aquadag. The 
conducting rod is charged to a positive potential of 270 
volts by an 8-micron gold-coated quartz fiber mounted 
on a side arm. When the fiber is connected to the positive 
potential, it is attracted to the rod, charging the system to 
the full potential of 270 volts. The fiber is then repelled 
by. the electrostatic charge of the system. As the central 
rod collects electrons formed by ionizing radiation in the 
argon gas of the chamber, its potential drops, and the 
fiber slowly moves back until it is close enough to be 
violently attracted by the image charge. The fiber then re- 
charges the rod and the charging pulse actuates the exter- 
nal electronic circuit. Thus each pulse from the chamber 
represents the collection of approximately 10 coulomb 
of electricity. Data are recorded in both analog and digital 
form; the subcarrier frequency on the analog output is devi- 


134 


ated for every 32 ion chamber pulses, and, in addition, the 
digital output for very low rate measures the time between 
two pulses of the ion chamber by counting the impulses 
from a standard time source and storing them in the digital 
accumulator by means of electronic gates. For very high 
pulsing rates, for example in the earth’s corona, the time 
between pulses is too short to store clock signals so every 
ion chamber pulse is fed directly to the digital accumulator. 
Ambiguities which may arise in intermediate regions can 
be resolved by comparison with the analog record. 

The second detector used in this experiment is a small 
Halogen Geiger counter (Anton type 302). This detector 
has been previously used on the Pioneer III and IV space 
probes. The counting rate is recorded in digital and analog 
form. The analog data are recorded with scale factors of 
2° and 2" and transmitted continuously on the same chan- 
nel as the ion chamber data. The three signals are dis- 
tinguished by their amplitudes. The digital data are handled 
in the same manner as the digital output of the ion cham- 
ber except that a scale factor of 21° is used. The ambiguity 
in the counting rate in high flux regions introduced by 
deadtime effects is resolved by the ion chamber since in 
high flux regions the ionization is high. 

Neither the ion chamber nor the Geiger counter has been 
deliberately shielded. The shields encountered by the in- 
cident radiation are the walls of the vehicle, an aluminum 
can surrounding the complete experiment, the walls of the 
detectors, and the other apparatus present in the vehicle. 
The minimum energy particles which can enter the de- 
tectors are about 1.25 mev electrons or 25 mevy protons. 
Due to the different shielding in different directions, the 
geomagnetic factors of the detectors are functions of en- 
ergy, measured during calibration experiments. 

Preliminary results of the Pioneer III and IV space 
probes indicate a large variation in the outer radiation belt, 
believed to be the result of several solar storms which im- 
mediately preceded the paunch of Pioneer IV. The Able-3 
vehicle will be in its orbit long enough so that it will be 
possible to correlate changes in the intensity of the radia- 
tion belts with solar activity. These changes could be due 
to changes in the amount of trapped radiation or to streams 
of charged particles from the sun passing in the vicinity of 
the vehicle but not necessarily being trapped in the ter- 
restrial magnetic field. 

The intensity of trapped radiation can be expected to 
vary with the magnetic field, and therefore the results of 
the radiation experiments will be correlated with the mag- 
netic field measurements. 

3) Proportional Counter Telescope (University of Chi- 
cago): About a year ago, it was decided to build an instru- 
ment which would measure high energy charged particles 
in the presence of strong bremstrahlung. The apparatus is 
composed of a triple coincidence counter system sur- 
rounded by 5 gm/cm? of lead. By making the counters 
operate in the semiproportional region, the deadtime of the 
counters was reduced to the order of one microsecond, and 
threefold coincidences due to charged cosmic ray particles 
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can be readily measured even though the detector is situ- 
ated in the low energy trapped radiation where bremstrah- 
lung might produce a high counting rate in the individual 
counters. To avoid the problem of directionality while 
passing through the earth’s field, the telescope was made 
essentially omnidirectional. The apparatus consists of seven 
cylindrical proportional counters, six in a concentric ring 
around the seventh, with the outer counters connected in 
two adjacent groups of three. The pulses from each group 
of counters are first amplified before entering the triple 
coincidence circuit. The pulse output of the triple co- 
incidence circuit is recorded with binary scalers having 
a storage of 2"! to 21° depending upon the particular flight 
requirements. The central counter, which forms the second 
counter of the triple coincidence, is also used independently 
to determine the single counter count rate for charged par- 
ticles inside the lead shielding. Thus the single counter 
provides a measure of the bremstrahlung and is used as a 
monitor of this radiation, in addition to providing the 
single count rate needed for correction of accidentals at 
very high levels of flux. The output of the single counter 
is likewise fed to a scaling circuit. Both the single-counts 
scaler and the triple coincidence counts scaler feed their 
outputs into binary storage units provided by STL, which 
may then be read out upon command. 

The minimum energy required to penetrate the triple 
coincidence counter with shielding is 70 mev for protons 
and 12 mev for electrons. Thus, only relativistic electrons 
will be measured as triple coincidence events and only pro- 
tons in excess of 70 mev will be registered. 

This detector system was used on the attempted lunar 
shot of November 8, 1958, in which the detector entered 
the lower edge of the Van Allen region near the equator. 
From the count rates in the triple coincidence detector it 
has been determined that near the lower edge of this 
zone the flux of energetic particles was approximately 
400/cm?-second and is composed of either protons above 
70 mev or electrons about 12 mev. The evidence favors a 
proton component. 


C. Supporting Experiments 


In addition to the satellite-borne radiation experiments 
described above, an extensive program of supporting ex- 
periments is planned by the university groups. Both uni- 
versities plan to continue their programs of balloon flights 
to monitor the cosmic ray intensity at the top of the atmos- 
phere. The University of Chicago maintains a network of 
ground stations running northward from the equator which 
will continuously monitor the cosmic ray intensity through- 
out the period of satellite data acquisition. All participants 
have planned an extensive alert network to insure maximum 
data acquisition from the satellite during times of intense 
solar activity. 


D. Objectives of the Experiments 


Although each of these experiments is independent and 
complete in itself, the results cannot be completely evaluated 
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separately. Each measures the integrated flux of particles 
above its energy cut-off. However, due to the different 
shielding used, each has a different cut-off, as shown in 
Table I. Thus it should be possible to obtain a three-point 
intensity vs energy spectrum. 


TABLE I 
Detection Threshold (mev) 
Experiment = 
Electrons Protons 
Scintillation Counter On225 2.0 
Geiger Counter—lonization Chamber ez 25 
Proportional Counter—Telescope 12 70 


Since the University of Minnesota experiment measures 
both ionization and counting rate with similarly shielded de- 
tectors, it is possible to compute the mean ionization per 
particle and thus obtain some indication of the composition 
and energy of the radiation. 

The principal objectives of the University of Chicago 
and University of Minnesota experiments require meas- 
urements from space probes as proposed under the Able-4 
program. The use of the same instrumentation on the Able-3 
satellite serves three purposes: first, it is identical with the 
equipment to be carried in the Able-4 vehicles and therefore 
represents a checkout of the system for these shots. Sec- 
ondly, there exists the possibility of simultaneous measure- 
ments from the Abe-3 and Able-4 instruments, which would 
provide information not possible from a single vehicle. 
Thirdly, the University of Chicago experiment provides a 
means for investigating the high-energy portion of the 
charged radiation trapped in the earth’s magnetic field. 

The experiments from the two universities cover two 
areas of physics: first, they are concerned with the solar 
production of charged particles of cosmic ray and lower 
energies and secondly, the University of Chicago experi- 
ment is designed to study the modulation of the extra 
solar cosmic ray particles by transient magnetic fields in 
the interplanetary medium. These measurements have a 
bearing upon the description of the solar-generated plasmas 
in the interplanetary medium. 

1) Solar Production of Charged Particles: The second 
problem is the observation of transient intensity increases 
of nonrelativistic particles which may occur around the time 
of the solar activity maximum. Already several events have 
been reported which suggest that solar flares may produce 
low-energy particles in great numbers and that these parti- 
cles are dominantly protons and alpha particles. These 
bursts from the sun are of considerable importance, not only 
for verifying the characteristics of the acceleration mecha- 
nism, but also for their use as probes to explore the mag- 
netic fields in the intervening interplanetary space. 

2) Magnetic Fields in Interplanetary Space—The Mod- 
ulation of Cosmic Ray Intensity: It is now well established 
that over the period of a solar cycle the flux of relativistic 
particles that reaches the earth may be reduced by more 
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than a factor of two at times of maximum solar activity. 
Aside from this 11-year cycle variation in which the in- 
tensity is reduced at the times of solar maximum, there 
exist sharp decreases of intensity that take place in the 
order of hours and usually represent somewhat smaller 
reductions of total intensity. There is good evidence indi- 
cating that these two phenomena represent the modulation 
of pre-existent cosmic rays—cosmic rays which come from 
the galaxy. We currently believe that the modulation takes 
place through disordered magnetic fields in interplanetary 
space. The questions now open relate to the location of 
these fields, their scale size, and whether the same fields 
account for both types of modulation effects (11-year and 
F-decreases). There is some evidence that the two mecha- 
nisms, although physically similar, involve phenomena of 
different sizes, namely, that the 11-year cycle requires the 
distribution of magnetic fields over the entire inner solar 
system, whereas the mechanism producing the Forbush- 
decreases may be more localized in the vicinity of the earth. 
A space probe in the interplanetary medium which shows 
that the cosmic radiation remains at the intensity observed 
near the earth would favor a large-scale modulation mech~ 
anism, whereas an observation showing that the intensity 
rises by factor two or more beyond the environment of 
earth would prove that the 11-year cycle effect is of local 
origin. 


III. Wuistter-Mope EXPERIMENT 
(STANFORD UNIVERSITY) 


A form of VLF radio noise, usually having a pitch 
decreasing with time over a period of about one or two 
seconds, has been observed for decades. It has been estab- 
lished that such a noise signal, called a “whistler,” is a 
pulse of VLF energy produced by a lighting stroke and 
propagated along the magnetic field lines from the source 
up through the ionosphere, reaching heights of several 
earth radii from the surface, and then guided by the mag- 
netic field down through the ionosphere to reach the sur- 
face at a point (the conjugate point) lying on the same 
magnetic line of force as the source. The pulse may then 
be reflected from the ground to make a number of trips 
along the field lines between the source and the conjugate 
point. One essential requirement for the existence of this 
propagation mode, called the “whistler mode,” is the pres- 
ence of electrons having a density of several hundred per 
cubic centimeter at all points along the path. Thus the 
whistler mode is connected with both magnetic intensity 
and electron density in regions well above the main portion 
of the ionosphere. 

A direct observation of naturally occurring whistlers 
from a satellite is subject to the disadvantage that the lo- 
cation of the source pulse is uncertain, and the identifica- 
tion of the ground-propagated pulse, if observed, is a 
matter of chance. Much of the information concerning 
natural whistlers as observed at ground stations is obtained 
from a study of their amplitude and frequency changes 
with time. To carry out similar observations from a satel- 


136 IRE TRANSACTIONS ON 
lite requires a wide-band receiver and either wide-band 
telemetry (several kilocycles of bandwidth) or airborne 
processing equipment. The whistler-mode experiment for 
the Able-3 program is based on the observation of the 
whistler-mode propagation of coded signals transmitted 
from existing VLF transmitters. The transmission of such 
“artificial whistlers” between a 15.5-kc transmitter near 
Washington, D.C., and a receiver located at the southern 
tip of South America has already been demonstrated. 

The airborne receiver for the Able-3 program is a tuned 
radio-frequency (TRF) receiver tuned to a frequency 
of 15.5-+0.05 ke. The receiver output amplitude modu- 
lates the appropriate telemetry channel. The transmitter is 
an existing unit transmitting keyed CW at a carrier fre- 
quency of 15.5 ke. A receiver bandwidth of 100 cps is 
used to give time resolution of approximately 10 milli- 
seconds for measurement of group delay. The telemetered 
information is compared with the outputs of 15.5-ke re- 
ceivers located at the STL tracking stations. Observations 
of the whistler-mode are also made from a Stanford Uni- 
versity ground station at the point conjugate to the trans- 
mitter location, i.e., from the southern tip of South 
America. Both the VLF airborne equipment and the VLF 
ground station equipment at the tracking stations have been 
provided by Stanford University. 

A comparison of the time at which a code group is 
transmitted with the times at which it is received in the 
satellite and at the ground conjugate point makes it possible 
to determine the transmission times over two parts of the 
whistler-mode path. Such information is directly related to 
the average group velocity over the path, which in turn is 
dependent upon the electron-density distribution along the 
path. Also, in the course of the experiment, the amplitude 
of the received waveform is included in the telemetered in- 
formation so that the experiment provides information 
concerning the attenuation of signals propagated via the 
whistler mode. The antenna used with the airborne receiver 
is the vehicle telemetry antenna to which the VLF receiver 
is multiplexed. 

When the digital telemetry is operating, the payload be- 
comes a phase-preserving receiver that permits measure- 
ment of phase velocity. Since the calibration includes 
measurement of antenna capacitance, some information 
pertaining to the dielectric constant of the medium may be 
possible. 

This experiment on the propagation via the whistler 
mode is directly applicable to the interpretation of radiation 
and magnetic field data obtained by the other experiments. 
Any large-scale disturbance in the magnetic field structure 
of the earth due to magnetic storms is likely to exert an in- 
fluence on the whistler-mode propagation. The experiment 
should also provide measurement of the physical size of 
propagation ducts for VLF signals, the coupling loss in 
propagating through the ionosphere, the magneto-ionic duct 
propagation cutoff in the region where the gyromagnetic 
frequency equals the signal frequency, and an over-all 
survey of VLF noise above the ionosphere. 
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The measurement of noise intensity is an important part 
of the VLF experiment. When analyzed in conjunction 
with the simultaneous data received from the corpuscular 
radiation, electron density, and magnetic field vector meas- 
urements, it may be possible to identify the mechanism 
which generates VLF energy. 


TV. Erectron Density ABOVE THE IONOSPHERE 


A. Introduction 


The determination of the contents of solar space has 
proved to be a puzzling and controversial problem, pri- 
marily for lack of good experimental data. The data on 
which to base an estimate of solar space electron densities 
are especially meager. 

One could, for example, begin with data from the solar 
corona. Using polarization measurements and surface 
brightness of the K-corona, an interpolation formula for 
the electron density, N(7), has been determined: 


N(r) = 108(37-18 are at 87-3!) 


where N(r) is electrons/cm® and 7 is in solar radii. This 
equation holds only in the plane of the solar equator and 
only to the vicinity of the orbit of Mercury, since presum- 
ably the solar corona should terminate there. However, 
by extrapolating r to the earth’s orbit even though one 
should not expect continuity between the corona proper 
and the interplanetary gas, a value for the electron density 
of N = 500 cm is obtained. This procedure, of course, 
is speculative since we are extrapolating over a range of r 
of nearly a hundredfold, and there is no theoretical reason 
for believing that N should be a function of r at this dis- 
tance. Nevertheless, this result is approximately the same 
as results obtained by other means. 

Another effect which gives a measure of the electron 
density above the ionosphere is atmospheric whistlers 
(Section III). Theoretically, to support this type of propa- 
gation requires a minimum electron density of about 500 
cm=* at the crest of the propagation path. Since whistlers 
have been observed where propagation paths apparently 
extend out many earth radii, electron densities of at least 
500 cm? presumably exist in interplanetary space close 
to the earth’s orbit. 

Using Doppler measurements with very high altitude 
missiles and satellites, observers have obtained measure- 
ments of electron densities several hundred kilometers 
above the maximum of the ionosphere. The data indicate 
that the electron density decreases much more slowly be- 
yond the maximum than it increases below the maximum. 
Using satellite data and extrapolating the model, N = 
NimaxC-™", Russian experimeners? find that the number of 
electrons above the ionosphere maximum is about 3.6 times 
that below the maximum (a result consistent with Faraday 


*I. L. Alpert, F. F. Dobriakova, E. F, Chudesenko, and B. S. 
Shapiro, “Electrons concentrated in the upper ionosphere as deter- 
mined from radio signals of the first earth satellite,’ Doklady Akad. 
Nauk. SSSR, vol. 120, pp. 743-746; 1958. 
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rotation experiments described below) and, that at alti- 
tudes of about 2000 to 3000 km, the electron density is on 
the order of 200 to 300 cm. They also express the sup- 
position that this is the altitude where the earth’s atmos- 
phere is continuous with interplanetary gas. 

Faraday polarization rotation measurements have been 
made with radio echoes from the moon.? These results are 
in agreement with the missile and satellite measurements 
mentioned above; i.e., there is apparently a very gradual 
tapering off of the earth’s ionosphere; they show that 
about 75 per cent of the electron content between the earth 
and the moon lies beyond the ionospheric maximum. 

Finally, the observed polarization of the zodiacal light 
has been interpreted as due to scattering of sunlight by 
electrons and dust particles in the neighborhood of the 
earth’s orbit. Since light scattered by both electrons and 
dust particles will be polarized, it is extremely difficult to 
separate the two. Van de Hulst® estimates that if the polari- 
zation due to dust scattering were negligible, the electron 
density near the earth’s orbit would be about 600 cm”*. 
Taking into account polarization due to dust he obtains 
values of 200 cm™® and 800 cm under two extreme as- 
sumptions. Russian experimenters claim that, based on 
their measurements of zodiacal light over the Libyan des- 
ert, there is essentially no polarization due to interplane- 
tary electrons. Their supposition is that the electrons in 
space near the earth’s orbit are due to corpuscular radia~ 
tion from the sun. Zodiacal light measurements should 
then be symmetrical about the plane of the solar equator 
and should vary in intensity in accordance with the sun’s 
activity. Since they did not observe these effects, they con- 
clude that sun-earth space is composed primarily of dust 
particles and has a negligible electron component. Black- 
well argues, on the other hand, that the above results would 
tend to prove that there is a stationary interplanetary gas 
(electrons and protons) independent of the corpuscular 
radiation from the sun. 


B. Doppler Difference Frequency Measurement 


The method used to determine the electron density above 
the ionosphere for Able-3 is similar to that employed by 
Seddon‘ for electron density measurements in the iono- 
sphere. It is based upon measuring the difference between 
two Doppler frequencies transmitted from the satellite. A 
CW signal transmitted from the ground at a frequency of 
fo = 401.8 mc is multiplied in the satellite by two factors 
p’ and q’ and these two coherent signals, f: and fo, are then 
retransmitted to the ground. The two multiplication factors 
are p’ = p’ =16/17 and q’ = pq, where q =2/7. Hence, 
the two frequencies transmitted from the satellite are 


fi = fop = 378.2 mc and fe =fopq = 108.06 me. 


2J. V. Evans, “The electron content of the ionosphere,” J. At- 
mos. & Terrest. Phys., vol. 11, pp. 259-271; 1957. 

°HC. Van de Hulst, “Vistas in Astronomy,” London, Eng., 
vol. 2; 1956. , , 

+7. C. Seddon, “Propagation measurements 1n the ionosphere with 
the aid of rockets,” J. Geophys. Res., vol. 58, pp. 323-335; Septem- 
ber, 1953. 
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The frequencies of the two received signals at the 
ground tracking station are then essentially 


(1) 
(2) 


Ri = (fo — Afo)p — Af, 
Re = (fo — Afo)pq — Afe 


where Afo, Af:, and Af, are the three Doppler frequencies. 
The Doppler frequency Af, is 


geet 


fo . Ono 
ndSo = ZY us “fe [ ae ass | (3) 
whi So C Sy OH 
where the integration is taken along the ray path, So, from 
the ground station to the satellite. S, is the component of 
the satellite velocity along the tangent to the ray path at 
the satellite, m, is the index of refraction along the ray path, 
and mx is the index of refraction at the satellite. At these 
high frequencies the index of refraction may be approxi- 
mated by 


, KN 
1 4 
n 2 (4) 
where NV is the electron density in m™*, K = 40, and f, is 
given in cps. Using this approximation, the Doppler fre- 
quency Af, becomes 


SS Ks alien 
Afo = So j | Naso + ie zy ass. (5) 
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Similar expressions hold essentially for Af; and Af. How- 
ever, the ray paths So, Si, and S» and the tangential com- 
ponents of the satellite velocity So, Si, S2 differ because 
the index or refraction is slightly different for the three 
frequencies transmitted. At the high frequencies being 
considered here, the ray paths are sufficiently close that we 
may to a very good approximation use S$, = S, = S, and 
So = S, = S:. With this approximation, the two received 
signals [(1) and (2)] become (after R» is multiplied by 


1/q): 
Ri = fop — Afop — Afi 


2S K 1 
= fo — sk P+ lus Hp e|(1 +—) (6) 
C foc p? 
R A 
ee ome le! 
q q 
2S K 1 
= pfs — nee +—2[ vus + El(1 + ) (7) 
foc 9’ 
where 
peed (aaa (8) 
may i 


Hence, we have a system of two equations with two un- 
knowns. S and (Ny + E/S), where E/S is to be considered 
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an error term. Solving (6) and (7) for the electron density, 
Ny, yields 


fo’ W E 
Nu aes 


K 1 1 1 ©) 
w(: es) of ee =) 
qg 
where the two measured quantities are the difference be- 
tween the Doppler frequencies. 


R A 
detent. Seat ye (10) 
and the Doppler frequency at 378.2 me, 
U = bfo oa Ry = pAfo aon Afi. (11) 


It is not necessary to measure the Doppler frequency U 
since it can be computed from trajectory data. However, 
measuring U is slightly more accurate and makes the sys- 
tem completely self-contained. The error term (E/S) 1878 
fictitious electron content along the ray path. It is shown in 
the next section that at certain times this term can intro- 
duce a significant error in the electron density measure- 
ment. To help us correct for this, we are making a Faraday 
polarization rotation measurement simultaneously with the 
Doppler measurement. 

The values of the parameters entering into (9) depend 
upon the frequencies being used, the electron densities ex- 
pected, and the orbit of the satellite. From nominal trajec- 
tory data, estimates can be obtained on the minimum and 
maximum range rates ($) which will prevail while making 
electron density measurements. These rates will vary from 
a minimum of about 10° feet per second to a maximum of 
about 3 X 104 feet per second. As noted above, the best 
available evidence indicates that the electron density will 
vary from a maximum of about 10° cm™ near perigee (200 
miles) to densities of about 10° cm™ at altitudes greater 
than a few thousand miles, although the minimum density 
is extremely uncertain since the available evidence is 
meager. The Doppler frequency shifts corresponding to 
the range rates being considered are, for a 400-mc signal, 
about 8 X 10? cps at minimum range rates and about 
2 X 10* cps at maximum range rates. For this experiment, 
the most critical parameter is the difference between the 
two Doppler frequency shifts, W (10). For minimum 
range rates and an electron density of 10® cm”, the mini- 
mum Doppler difference, W, is about 10-* cps. It is the 
value for W under conditions of solar space where electron 
densities are low, range rates are low, but where measure- 
ment times can be long. 

The ground equipment used to measure the Doppler 
frequency difference is quite simple. It is designed to work 
with two Able-3 Microlock receivers at 378 and 108 mc. 
The inputs to the equipment are the VCO outputs which 
contain the Doppler frequencies from the receivers. After 
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suitable multiplication and mixing, the difference in phase 
between these two signals is measured and recorded on a 
Sanborn chart. This gives W. A separate measurement of 
the Doppler frequency U at 378 mc is made and recorded. 
(The range rate, S', is also known from trajectory data.) 
With these two quantities known, the electron density can 
be computed from (9). 


C. Faraday Rotation Measurement 


The error term (E/S') entering (9) is due to the time 
change of the electron content along the radio path from 
the satellite to the tracking station. The change in electron 
content is caused by changes in the sun’s ionizing radiation 
during every 24-hour period due to the earth’s rotation, 
and to ionospheric scintillation caused by the drift of elec- 
tron patches across the propagation path. Since the satellite 
is many thousands of miles from the tracking station dur- 
ing electron density measurements, its angular position 
with respect to the earth changes by only a few degrees 
during the time of measurement, and the propagation path 
through the ionosphere is essentially constant. The effect 
of the growth and decay of the electron content can be esti- 
mated from Faraday polarization rotation measurements 
of moon echoes. These show that during major portions of 
every 24-hour period, a fictitious electron density on the 
order of 10? cm™ can be introduced. This is approximately 
the same as the expected electron density in solar space. 
During other parts of the day, especially just before dawn 
when the ionosphere is relatively unchanging, the error 
term is negligible. 

An indication of the error introduced by ionospheric 
scintillation can be deduced from radio star® scintillation 
measurements. These data are still quite sketchy, but indica- 
tions are again that ionospheric scintillation can introduce 
a fictitious electron density comparable with the expected 
actual density. However, the scintillation phase fluctuations 
presumably have a mean of zero; 1.e., sometimes they add 
and other times subtract a fictitious electron density. Con- 
sequently, over a long enough measurement time, this error 
can be “averaged” out. 

To obtain further quantitative data on this error term and 
thus correct the electron density measurements, the time 
rate of change of the Faraday polarization rotation is 
made at 108 mc simultaneously with Doppler difference 
measurements. At this high frequency, the Faraday polari- 
zation rotation angle for quasi-longitudinal propagation is 
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where H is the earth’s magnetic field in gauss along the 
path and 9 is the angle between the direction of propaga- 
tion and the magnetic field. At the great distances that 
electron densities are measured, the time rate of change of 
the Faraday rotation is approximately (the term due to 
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satellite motion is negligible at great distances) 
oN 
(H cos @) i} — ds 
s Ob 


where (H cos 0) is a suitable average of this quantity 
through the ionosphere. This gives a measure of the error 
term (E/S), contained in (9). The rotation, ¢, is measured 
both manually and by recording the phase difference of 
the two circularly polarized components of the 108-mc sig- 
nal from the satellite. 


2.3 X 104 
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V. SEARCH-COIL AND FLUX-GATE MAGNETOMETERS 


The magnetic field of the earth is usually represented by 
a model consisting of an eccentric dipole plus higher order 
multipoles. Up to the present time, th: data used in the 
establishment of this model have been obtained from sur- 
face measurements and a few isolated rocket measure- 
ments. More recently, the extraterrestrial magnetic field 
has been sampled by the Pioneer I space probe, the telem- 
etered data from which have indicated departures from 
the rate of variation to be expected from a spherical har- 
monics expansion of the surface measurements. These de- 
partures originate in extraterrestrial phenomena, such as 
the controversial ring current, or the magnetic field of the 
solar corona. 

A satellite is inherently well adapted to map the earth’s 
magnetic field at high altitudes. A measurement of the 
vector magnetic field from a high-altitude satellite has the 
advantage of being unaffected by local surface anomalies 
on the earth. In addition to the mapping of the steady-state 
magnetic field, the field observations from the Able-3 ve- 
hicle permit the direct investigation of the ring current, of 
transient phenomena such as auroral current, and of vari- 
ous hydromagnetic effects. 

One such hydromagnetic effect is the hypothesized prop- 
agation of hydromagnetic waves along magnetic field lines. 
Turbulence at the transition between the coronal wind and 
the geomagnetic field, or turbulence in the coronal wind 
itself, may provide a mechanism for the generation of such 
waves in the geomagnetic field. Such effects may have been 
observed in the Able-1 program. Although the effects of 
interest occurred at a distance of about 10 earth radii from 
the earth’s center on the Pioneer I flight, magnetic condi- 
tions were exceptionally quiet on that day. On a magneti- 
cally disturbed day, the transition between the geomag- 
netic and solar magnetic fields is expected to occur at about 
5 earth radii, thus making the effect observable from the 
Able-3 orbit, which has an expected apogee of about 7 
radii. For several reasons, therefore, a continuous moni- 
toring of the extraterrestrial geomagnetic field, during a 
period of time extended enough to include a variety of 
forms of solar activity, is of great fundamental value. 

The two magnetometers carried in the vehicle operate on 
different principles. A search coil mounted near the pe- 
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riphery of the vehicle measures the component of the mag- 
netic field perpendicular to the vehicle spin axis, and a flux 
gate also near the periphery, oriented parallel to the spin 
axis, measures the component of the magnetic field along 
this axis. By using the output of the television sensor to 
indicate the vehicle rotational position at the time of maxi- 
mum search-coil output, the direction and magnitude of the 
component of the magnetic field normal to the vehicle spin 
axis is determined. This information, combined with the 
output of the flux-gate magnetometer, determines the mag- 
nitude and direction of the vector magnetic field provided 
that the orientation of the vehicle spin axis is known. 

The search-coil magnetometer is similar to that used in 
the Able-1 program. Its sensing element consists of a mu- 
metal core wound with 5000 turns of number 40 copper 
wire. The output of the coil is coupled to a transistorized 
amplifier tuned to the nominal spin rate of the vehicle. In 
a steady magnetic field, the amplifier output is a sinusoid 
at a frequency equal to the spin rate. This output is applied 
to a subcarrier oscillator to modulate a telemetry channel. 
In the Able-1 program, the greatest possible sensitivity and 
dynamic range were desired. To attain this objective, a 
very large amount of AGC feedback was employed in the 
amplifier to extend the dynamic range. Such extension was 
obtained at the expense of the high-frequency response of 
the system. In the Able-3 program, on the other hand, the 
amplifier is modified to reduce the dynamic range to a value 
of approximately 300 to 400. This reduced dynamic range 
provides precise measurements over a range from about 1 
milligauss, the smallest expected field for the Able-3 orbit, 
up to 300 to 500 milligauss, the surface value of the geo- 
magnetic field. The reduction in dynamic range results in 
an improvement in the transient response of the search- 
coil magnetometer. 

The flux-gate magnetometer, purchased from the Schon- 
stedt Engineering Company, has a high-frequency cutoff 
of about 10 cps and, as a consequence, relatively rapid 
changes in the magnetic field amplitude can be followed. 
The flux-gate magnetometer is, however, subject to dc 
drifts, and is temperature sensitive. The search-coil mag- 
netometer, on the other hand, being an ac device, is essen- 
tially independent of temperature and other causes of drift. 

In addition to providing a complete description of the 
vector magnetic field, the two devices complement one an- 
other in other respects. In particular, the flux-gate mag- 
netometer, and thus the flux-gate output, is more useful in 
the study of rapid transient variations in the magnetic 
field. On the other hand, the insensitivity of the search- 
coil magnetometer to drifts makes it the better device for 
the observation of possible secular variations. 

Although neither the flux-gate nor the search-coil mag- 
netometer possesses the accuracy offered by the proton 
precession magnetometer, the search-coil and flux-gate 
magnetometers each measure the magnetic field continu- 
ously, while the proton precession magnetometer samples 
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at intervals of about 10 seconds. Furthermore, the power 
requirements for the proton magnetometer are severe. 
Since hydromagnetic phenomena constitute a very impor- 
tant class of magnetic phenomena, a continuously-reading 
magnetometer is of the greatest utility. 


VI. Rapio SCINTILLATION EXPERIMENT? 


The object of the radio scintillation experiment is to 
gain information about the size, shape, orientation, mo- 
tions, and stability of ionospheric irregularities which 
cause amplitude and phase scintillations in radio waves. 
Almost all information has been derived from measure- 
ments in the field of radio astronomy using radio stars as 
sources, which have input signal-to-noise ratios of ap- 
proximately 13 db. Therefore, the expected signal-to-noise 
ratios of order 25 to 30 db offered an excellent opportunity 
for studying the important physical parameters of the iono- 
spheric irregularities. 


A. Equipment 


The two 40-foot diameter radio telescopes on the present 
1500-foot east-west base line at Table Mesa, Colo., are 
used as the two elements of the phase-lock, phase-sweep 
interferometer system. Two STL receivers of the phase- 
lock type are used. Each receiver consists of four super- 
heterodyne stages. The first, second, and fourth stage local 
oscillators of the two receivers are locked together, with the 
first receiver phase-locked to the incoming signal. The 
third stage local oscillators will be maintainend 1000 cps 
apart by means of a 1000-cps local oscillator, so that the 
RF phase difference of the incoming signals at the two an- 
tennas show up as phase modulation on a 1000-cps carrier. 
It is for this reason that the term “phase-sweep” is applied 
to the detection scheme. The bandwidth of this 1000-cps 
signal is 100 cps. This 1000-cps signal is phase-compared 
with the 1000-cps oscillator to obtain the desired angular 
scintillations out to 100 cps. The amplitude of the 1000-cps 
signal is also continuously recorded to obtain a continuous 
record of the signal power received by the system. 

As explained below, two such phase-lock, phase-sweep 
interferometers permit simultaneous observations over 
various baseline lengths and orientations. The second unit 
utilizes one of the 40-foot paraboloids and a portable Yagi 
antenna, but otherwise is identical to that described above. 


B. Expected Signal-to-Noise Ratio 


The signal-to-noise ratio to be expected at the input of 
a single paraboloid receiver system of appropriate polari- 
zation is given by 


(2) = cs) ol 
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® This experiment is part of a continuing BMD-STL weapon 
systems effort. 
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where P is the effective transmitted power, the second term 
is the space transmission loss factor, the third is the effec- 
tive area of the paraboloid of radius, r, and the last term 
is the effective noise power (receiver plus sky noise) in 
the input bandwidth, Bin. 

At the output of the receiver, the signal-to-noise ratio of 


such a system is 
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If two such antenna-receiver systems are coupled to 
form a phase-sweep interferometer whose input and out- 
put bandwidths are both equal to B, the signal-to-noise 
ratio at the output is 
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In this experiment, the radiated signal is linearly polar- 
ized. To avoid periods of loss of signal due to Faraday 
rotation, circularly polarized antennas will be used, with a 
resultant loss of receiving efficiency of a factor of two. 


(5) = Gr Tl 


The signal-to-noise ratio on the proposed phase-sweep 
interferometer system is therefore a factor s(=5 db) 
worse than that obtained by a single-antenna, total-power 
system of correct linear polarization. This loss of signal- 
to-noise ratio is, however, more than compensated for by 
the improved stability of the phase-sweep system, whose 
output is relatively insensitive to variations in sky noise, 
receiver noise, or receiver gain. 

Substituting P = 500 mw, R = 36 X 10° meters (apo- 
gee), Nx R? = 60m’, T. = 1000°K, B = 100 cps in (14), 
we obtain an output signal-to-noise ratio of approximately 
26 db. This compares with an output signal-to-noise ratio 
of about 14 db for the Cygnus source when using the same 
antennas in a phase-sweep interferometer system at an 
output bandwidth of 0.5 cps. 


(14) 


C. Use of Portable Yagi Antenna 


In the case of where the effective spacing, D, of the 
antennas of an interferometer is small compared with the 
0.5 correlation distance in the incident diffraction pattern 
of radio energy, one can assume that the phase gradient is 
uniform between the two antennas. Under such circum- 
stances, the rms angular scintillation observed on a single 
antenna is given by 


r 
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where ®,ms is the root mean square value of the irregular 
deviations in relative phase of the signal at the two an- 
tennas. 

In the case where the amplitude and phase scintillation 
are not well correlated at the two antennas, the actual ir- 
regular angular deviation of the incident wave front would 
be significantly greater than that predicted by (15). 

While it is probable that (15) can be used on the present 
radio-star data without great error, any search for high- 
speed scintillations, implying structure that is small com- 
pared with 1500 feet, means that (15) cannot be used. To 
overcome this difficulty (by actually determining the scale 
of the structure in the diffraction pattern), measurements 
are made at various antenna separations, using a simple 
portable Yagi antenna and one of the 40-foot paraboloids 
as the two elements of an interferometer. Since there is 
high probability that the scale of the diffraction pattern 
across the ground is a function of direction, observations 
are being made with variable orientation of baseline, as 
well as with variable baseline length. (For example, such 
observations could, under certain circumstances, permit the 
simultaneous measurement of the angular scintillations in 
azimuth and in elevation. ) 

The signal-to-noise ratio of the Yagi-paraboloid combi- 
nation is some 4 db worse than for the two-paraboloid in- 
terferometer ; if necessary, this deterioration in signal-to- 
noise ratio can be overcome by using a narrower tracking 
bandwidth (— 25 cps). 


D. Expected Capabilities of System 


1) Measurements of Relative Phase: 

a) Effect of finite signal-to-noise ratio: The expected 
signal-to-noise ratio of about 26 db (400:1) at apogee re- 
sults in an rms phase jitter of approximately 1/400 radian 
or 0.2° or relative phase, for an input-output bandwith of 
100 cps. 

b) Phase stability of electronics: Brief tests with the 
electronics of the existing ratio-star interferometers sug- 
gest that the minute-to-minute relative phase stability of 
the electronics is significantly better than 0.2°, and there- 
fore is nota limiting factor. 

c) Stability of electrical centers of paraboloids: The 
phase-centers of the dishes vary irregularly by one or two 
orders of magnitude more than 0.2° under extreme wind 
conditions. Only during calm conditions is it likely that 
phase stabilities of 0.2° are possible; a median value is 
probably about 1.0°. 

d) Tropospheric limitations: Discussions with Dr. M. C. 
Thompson of CRPL suggest that at low angles of eleva- 
tion (< 10°) the median relative rms phase shift due to 
the troposphere is of the order 0.1°. At higher elevation 
angles the effect should be less. 

e) Frequency stability of system: In general, the satel- 
lite signals at the point where they mix will have traveled 
several wavelengths further via one antenna than via the 
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effects upon the relative phase of the signals ; however, for 
a 10-wavelength difference, the relative phase would 
change by less than 0.1° provided the stability of trans- 
mitter was better than 1 part in 3.6 X 10* (say 3 in 10°). 
Minute-to-minute frequency stability of this order is avail- 
able. 

f) Conclusions: While under typical conditions the rela- 
tive phase measurements are limited by mechanical distor- 
tion of the paraboloids due to fluctuating wind loading, 
under very favorable conditions it should be possible to 
reach a sensitivity of 0.2° rms phase shift. For large-scale 
irregularities, this corresponds to a rms angular glint of 
about 0.4 second of arc. 

2) Measurement of Amplitude Scintillations: 

a) Effect of finite signal-to-noise ratio: The expected 
signal-to-noise ratio at apogee is 26 db for a 100-cps input- 
output bandwidth. This implies that the limiting sensitivity 
to amplitude scintillations (without additional integration ) 
is 1 part in 400 (= 0.01 db). 

b) Effect of fluctuations in receiver gain: While minute- 
to-minute rms stability of the order 0.01 db can be obtained 
on sensitive radio astronomy equipment, it is by no means 
certain that this can be achieved using the proposed phase- 
lock, phase-sweep interferometer system. However, it is 
possible that a figure of 0.03 db rms (0.6 per cent) can be 
obtained. 

c) Effect of satellite rotation: Satellite rotation can in- 
troduce changes in received signal by varying the polariza- 
tion and/or the power radiated towards the receivers. The 
former effect is minimized by the use of circularly polar- 
ized feeds; the effect of the latter depends upon the axis 
of rotation of the satellite relative to the axis of the dipole 
antenna. Any effect of satellite rotation should be periodic 
in nature, and can therefore (to a first order) be allowed 
for. 

d) Conclusion: Given adequate stability of the effective 
transmitted power, it should be possible to measure the 
moment-to-moment fluctuations in received signal strength 
to a limiting sensitivity of about 0.03 db. 


VII. Imace TRANSMISSION SYSTEM 


The present knowledge of the solar spectrum, together 
with the photo response characteristics of an image scan- 
ning device, allows for the determination of reflection 
characteristics of planetary or lunar surface. These data, 
accumulated over a broad band ranging from the infrared 
through the vacuum ultraviolet, can reveal information 
concerning the composition of planetary atmospheres, 
cloud cover, and general physical characteristics of the sur- 
face. There an image transmission system capable of op- 
erating under the weight, power, and bandwidth limitations 
of the Able-3 vehicle has been included in these payloads. 

The optical unit is a Newtonian telescope incorporating 
a General Transistor Corporation type 2N469 phototran- 
sistor mounted behind a stop with a pinhole aperture. The 
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plane of the stop is located at the focus of a two-inch aper- 
ture spherical mirror which looks out through a hole in the 
payload shell. The output voltage of the phototransistor is 
amplified, stretched in time, and telemetered in analog form 
to the ground stations. The transmitted output is a quan- 
tized function of the phototransistor output giving eight 
possible brightness levels. 

The 2N469 phototransistor is sensitive to wavelengths 
between 2000 and 17,000 A. Its response peak is in the 
8000 and 16,000 A region and consequently is primarily a 
near-infrared detector. 

The angular field of view of the instrument is 0.5°. Un- 
less a portion of a planet lies within the field of view, the 
output signal of the phototransistor is zero. 

The spin axis of the vehicle is roughly in the plane of 
its orbit. Therefore, it is possible to make use of the orbital 
motion and the vehicle spin to obtain a two-dimensional scan 
of the brightness of a planetary disk. When a planet lies 
within the line of sight of the system, a strip of the plane- 
tary disk of 0.5° angular width is scanned by the system. 
The brightness of this disk is measured by the voltage out- 
put of the phototransistor as it scans the disk. The system 
is designed to telemeter this voltage output pattern to a 
ground station, where the original intensity pattern can 
then be reproduced. 

Due to limitations introduced by the narrow informa- 
tion bandwidth available for telemetered data, it is impos- 
sible to transmit the intensity pattern from a complete scan 
during a single spin cycle of the vehicle. Therefore, it is nec- 
essary to rescan the same area many times and measure 
and telemeter one portion of the information during each 
spin cycle. 

In practice each scan line is divided into 62 spots of 0.5° 
angular diameter. The scan line starts at the limb of the 
planet and covers 31° of the spin motion of the vehicle. No 
data are transmitted until the system first sees the limb of 
the planet. Then the brightness of the first spot of the scan 
line is measured and transmitted during the remainder of 
the spin cycle. Each time the limb of the planet comes into 
view, a timer is started which selects the portion of the 
scan line whose intensity is to be measured. For example, 
on the tenth spin revolution after the planet has come into 
view the timer instructs the system to wait until the vehicle 
rotates through 9 X 0.5° or 4.5° after seeing the limb of 
the planet before measuring and transmitting the bright- 
ness. The timer is then advanced and the process repeated 
for 62 spin revolutions. Then the timer is reset to zero and 
for two spin revolutions, a sync signal is transmitted to 
show that this has been done. Thus, it takes a total of 64 
spin revolutions to scan a single scan line. 

During the time required to scan one line, the orbital 
motion of the vehicle changes the portion of the planet in 
the field of view of the image transmission system. Thus, 
when the scan of a single line is completed the next strip of 
the disk scanned is shifted slightly from the original strip. 
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In addition, the first strip is inclined at an angle since all 
the intensity points are measured at different times. Thus, 
it is seen that the reproduced image consists of a series of 
strips. each of which covers a portion of the disk of the 
planet. Since each strip has the same angular length, they 
may include only a portion of the planetary disk or extend 
beyond it. 

Under ideal conditions these strips are exactly adjacent, 
but for a satellite moving in an elliptic orbit about a planet, 
it is impossible to obtain a scan of the planetary surface in 
which all scan lines are exactly adjacent. In the case of the 
Able-3 vehicle, the image transmission system is oriented 
in such a way that the gaps and overlaps are never greater 
than half a scan line width. The scan lines will always 
cover the complete disk of the planet. The overlap of gap 
between lines can be calculated from the orbital parameters 
of the actual satellite orbit and may be taken into account 
when reducing the data. 

The nominal spin rate of the vehicles is 2.8 rps. At this 
spin rate, there is no overlap between 62 spots of a single 
scan line. However, if the spin rate is slower or faster than 
this nominal value, there will be overlap or gaps, respec- 
tively. If the vehicle precesses, the scan lines will be greatly 
distorted. This would lead to a major data reduction prob- 
lem in attempting to reconstruct the image of the planetary 
disk. A mercury damper ring is installed in the payload to 
minimize precessional motion. 


VIII. MicroMeteorITrE MoMENTUM SPECTROMETER 


Visual and radio observations of sporadic meteor trails 
in the earth’s atmosphere have established the presence in 
space of small grains of matter ranging in size from a few 
grams down to the order of 10° gram. The results of many 
measurements have indicated that the density of these 
grains in space is approximately inversely proportional to 
their mass ; that is, in any volume of space the average total 
mass of particles of a given size is the same for all sizes of 
grains. It is not known, however, whether this empirical 
law can be extended to dust particles whose masses are 
much smaller than 10° gram since such minute particles 
produce insufficient ionization in the atmosphere to be 
studied by radio or optical means. 

It is of considerable interest to determine the properties 
of the micrometeoritic material, not only to extrapolate the 
mass distribution to the millimicrogram region but to pro- 
vide a more substantial basis of fact for theoretical inter- 
pretation of other, related, astrophysical phenomena. There 
are considerable data available as a result of many precise 
astronomical observations which depend for their inter- 
pretation on assumptions regarding the nature of the inter- 
planetary and interstellar medium, assumptions which are 
further reflected in the form of theories of galactic struc- 
ture and evolution. Clearly, direct measurements of the 
physical properties of the interplanetary dust particles will 
be scientifically invaluable. 
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Although a complete set of measurements (charge, mass, 
velocity, composition, magnetic moment, etc.) is precluded 
by present payload restrictions, an initial step in such a 
study has already been made with the successful operation 
of the micrometeoritic spectrometer in the Able-1 program, 
and will be continued with the Able-3 flight. The relatively 
low density of the micrometeorite population requires a 
long sampling time in order to make statistically reliable 
conclusions concerning the interplanetary medium. In this 
respect, Able-3 should provide considerably more valuable 
data than did Able-1. Furthermore, because of the long 
observation time available, one can reasonably expect to 
obtain information on the spatial and temporal distribution 
of micrometeorites of both high and low momentum, indi- 
cating whether the particles in each momentum range are 
uniformly distributed or occur predominantly in localized 
regions of space and whether or not fluctuations in the 
density of particles are correlated with the occurrence of 
meteor showers as observed on the earth. Moreover, if an 
enhanced flux can be correlated with sporadic E layer ioni- 
zation, this will furnish valuable information concerning 
cosmic dust as a mechanism for ion production in the E 
region. 

This experiment, performed under the direction of the 
Air Force Cambridge Research Center, separates the mi- 
crometeorite flux into high and low momentum components 
in a simple way by having two levels of sensitivity of the 
detector. A particle striking the diaphragm of the spec- 
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trometer will transfer its entire momentum inelastically to 
the diaphragm, and the 100-kc component of the resulting 
impulse is picked up by a microphone containing a piezo- 
electric crystal which rings at this frequency ; the pulse is 
then amplified and the envelope detected. A two-way chan- 
nel (corresponding to two momentum thresholds) pulse 
height analyzer classifies the pulse according to whether it 
was produced by a low momentum (> —5 X 10% gm- 
cm/second) or a high momentum (> —5 X 10° gm-cm/ 
second) particle. The spectrometer has been calibrated 
with both low-velocity and high-velocity grains ; however, 
since the highest velocity of the calibration particles obtain- 
able in the laboratory (4 km/second) is considerably 
smaller than the probable velocities (11 to 70 km/ second ) 
of micrometeorites impacting on the satellite, the response 
of the detector to the momentum of actual micrometeorites 
is somewhat uncertain. The necessity of extrapolating the 
momentum calibration to high velocities together with the 
fact that the detector does not respond uniformly to par- 
ticles striking at different points of the diaphragm results 
in an estimated uncertainty of at least a factor of three in 
the momentum measurement. Any inferences made about 
the masses of the micrometeorites will be subject to even 
larger errors since the velocities of the particles can only 
be estimated. In spite of the relatively large uncertainties 
expected in the measurements, the experiment remains 
sufficiently accurate to contribute significantly to our 
knowledge of the properties of micrometeorites. 


Able-3 Payload Design* 
JOHN E. TABER} 


Summary—This paper discusses the Able-3 payload instru- 
mentation. The design parameters and functional characteristics 
of this instrumentation are discussed in terms of the two goals 
of this program: to build a payload which gathers as much 
scientific information as possible concerning the space environ- 
ment about the earth, and to develop and check out communica- 
tions, power supply, and temperature control systems suitable 
for subsequent interplanetary space probes. The basic electronics 
of the satellite and their relationship with the experimental equip- 
ment as well as with structural and other aspects of the satellite 
are discussed in some detail. Particular emphasis is placed on the 
communications system, and, in particular, on the new digital 
telemetry system and its value to deep space exploration. 


* Manuscript received by the PGMIL, July 21, 1959. . 
* Space Technology Labs., P.O. Box 95001, Los Angeles 45, Calif. 


INTRODUCTION 


Y THE TIME this paper appears in print, the 
B Able-3 launch will have taken place, and other re- 

lated vehicles will follow shortly. At the date of 
writing, the payload has been sufficiently well defined to 
warrant its description nearly two months before scheduled 
launch. 

The Able-3 payload, an earth satellite in a highly ellipti- 
cal orbit, has two purposes—that of studying the space 
environment about the earth, and that of providing a test 
vehicle to evaluate components and techniques developed 
for use in interplanetary space probes. The four-stage 
Able-3 vehicle consists of a Thor first stage, an Aerojet 
10-101 second stage, an ABL 248 solid propellant spin- 
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stabilized third stage, and a spin-stabilized fourth stage or 
payload weighing about 140 pounds. The first three stages 
are used to place the payload in a large geocentric ellipti- 
cal orbit whose apogee is about 20,000 nautical miles and 
whose perigee is about 150 nautical miles. A five-pound 
solid propellant injection rocket in the payload can be used, 
if necessary, to increase perigee and thereby increase the 
payload’s lifetime. 

Had the design of a Thor-Able launched payload been 
concerned solely with studying the space environment about 
the earth, it could have proceeded in a rather straightfor- 
ward manner and, with the exception of the number and 
complexity of the experiments carried, would have incor- 
porated many features found in payloads already orbiting 
the earth. An already developed FM-PM telemetry system 
would have been used, power supply design would have 
been conventional, and heat control could have been handled 
be proven techniques. The additional requirements that the 
payload serve as a “flying laboratory” for the evaluation 
of components and techniques useful in interplanetary 
space probes, however, has supplemented the usual design 
requirements with a new and unusual group of problems. 
Most of the new problems which arise can be associated 
either directly or indirectly with the need to transmit intel- 
ligence over multimillion miles, or with the severity and 
duration of both changing and static environmental con- 
ditions. 

The range equation, 

GrGr 
ees 

(4r°-R)?L 
expresses the functional relationship between the power 
received on the ground and the power transmitted in the 
payload, the payload and ground antenna gains, the wave- 
length used, range, and losses due to polarization, diplexers 
and cables. A second equation, 

Py = k[T. + To(n — 1) |B, (2) 
expresses the noise that will be superimposed upon the re- 
ceived signal as a function of Boltzmann’s constant, the 
temperature of outer space and at the receiver, the re- 
ceiver’s noise figure, and the bandwidth of the receiving 
system. The final relationships which complete the specifi- 
cation of the communication problem are the required 
signal-to-noise ratio at the receiving site as determined by 
the quality of data desired and the bandwidth needed as 
a function of the desired rate of information transmission. 
Moreover, this latter functional relationship varies with 
the degree of sophistication employed in choosing a modu- 
lation method. 

Any attempt to substitute realistic values of the param- 
eters into this system of equations quickly demonstrates 
that a conflict of interests exists. For instance, by choosing 
the largest ground antennas available (1.e., the 250-foot 
radio telescope of the University of Manchester), the most 
sensitive receivers in existence today, a payload antenna 


Pr= (1) 
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gain consistent with the present state of the art, and a 
frequency which minimizes extraneous losses, it is apparent 
that only a few bits of information per second can be con- 
veyed over interplanetary distances using a payload trans- 
mitter whose output power is as small as 100 watts. From 
this example two things become clear: first, the system 
designer must accept the need for a payload transmitter of 
relatively large power, together with its attendant power 
supply problems; and, second, any improvement in the 
modulation efficiency or in elimination of pretransmission 
redundancy will provide a better utilization of the limited 
transmission rate capability. 

The high transmitted powers required, coupled with the 
duration of interplanetary flights, precludes the develop- 
ment of a payload whose power is derived solely from 
conventional batteries. Since the state of the development 
of atomic power supplies does not permit their serious con- 
sideration at this time, the payload designer must turn to 
solar cells as a necessary choice for the source of primary 
power. The solar cell configuration discussed later in this 
paper approaches the greatest solar cell surface area which 
could be accommodated. It is capable of delivering between 
15 and 20 watts continuously to the payload. 

Early in the initial system design a value of 150 watts 
of payload transmiter power was specified for use in inter- 
planetary probes because on the one hand the largest trans- 
mitter and power supply which could be carried was limited 
to this value ; and at the same time using existing techniques 
and components it was sufficient to provide reliable com- 
munications over distances up to 55 million miles. After 
taking into account the efficiency of the power amplifier 
and of the static converters needed to develop the high 
voltages required, it was found that a primary power sup- 
ply capable of delivering peak powers of 500 watts would 
be required. 

It is thus immediately clear that the transmitter must 
operate at a reduced duty cycle. This, in turn, implies that 
a power storage system capable of accumulating power 
over extended periods of time and discharging it at very 
high rates must be carried. The payload weight which is 
available for a storage battery determines the maximum 
storage capacity and therefore the maximum ON time for 
the transmitter. 

The problem of how to obtain maximum information 
from the scientific experiments is a consideration which 
arises from the transmitter duty cycle limitations. Fortu- 
nately many of the data obtained are most meaningful in the 
form of events-per-unit-time. Micrometeorite impacts and 
radiation particle impingements are of this type. While it 
is desirable to keep the unit of time at some reasonably 
small value, information is obtained, though with greater 
smoothing, even when the unit of time is the several hours 
between successive transmissions. Hence there is a need 
for accumulating the number of events as they occur and 
later transmitting the total to the earth. 
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The information content explicit in analog data, how- 
ever, is so great as to preclude its storage in the payload 
for later transmission. Therefore experimenter must con- 
tent himself with periodic spot checks of such quantities as 
magnetic field strength and direction and payload tem- 
perature. 

Another factor implicit in the use of a low duty cycle 
transmitter is the need for an automatic controller or a 
command system capable of turning the transmitter on and 
off. Inasmuch as a Doppler transponder was to be carried 
aboard the payload for tracking purposes, it was logical to 
decide to use the greater flexibility inherent in a command 
system. The command system is also used to fire the injec- 
tion rocket, to change the rate of transmission of tele- 
metered information, to turn on a simplified television sys- 
tem, and other functions. 

Using the vehicle capabilities existing today, it can be 
anticipated that interplanetary distances will require 
months to traverse. Thus, equipment designed for inter- 
planetary space probes must have the reliability to endure 
unattended through the launching and for several months, 
or even years, in an airless environment while maintaining 
reasonable equilibrium temperature. 

The Able-3 payload incorporates in its design solutions 
which are felt to be adequate for all of these problems. 
While these solutions will not be evaluated under actual 
interplanetary travel conditions, it is believed that enough 
information will be obtained to increase greatly our con- 
fidence in the success of later launchings of truly inter- 
planetary probes. Two photographs of the interior of the 
payload are given in Fig. 1. 


COMMUNICATIONS AND TRACKING 


Transmitting and receiving equipment within the pay- 
load permits two-way communication. Payload transmitters 
convey telemetry information to the earth and a payload 
receiver permits the reception of earth-transmitted com- 
mands. When interconnected coherently, the payload re- 
ceiver and one of the transmitters form a transponder cap- 
able of providing both range and range rate information 
as an aid in tracking the payload. 

For communication to the earth, three transmitters are 
employed. One transmitter accepts an RF signal from the 
payload receiver, multiplies it eight times in frequency, and 
amplifies it to a 5-watt level. In the process the resultant 
signal, whose frequency is 378 mc, has been modulated 
with a 1024-cps subcarrier containing the time-multiplexed 
pulse-code-modulated output of a digital telemetry system. 
Biphase modulation is employed to impress the telemetry 
output on the subcarrier. This system, when supplemented 
with a 150-watt power amplifier, whose development was 
not included for the Able-3 program, forms the basic 
vehicle-to-earth portion of the interplanetary communica- 
tion system. The transmitter, about 3 by 6 by 1.6 inches, 
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Fig. 1—(a) Able-3 payload (top view). 
(b) Able-3 payload (bottom view). 


weighs about a pound, and has an efficiency of about 15 per 
cent. 

Each of two transmitters conveys analog data by means 
of a six-subcarrier FM-PM telemetry system operating at a 
carrier power of 100 milliwatts and a frequency near 108 
me. The transmitters together weigh one pound and con- 
sume about 1.5 watts of power. 

Except for some vehicle performance inputs, both tele- 
meters convey the same information, and thus the older 
and more proven FM-PM system can be used to monitor 
the performance of the new digital telemetry system. 
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The payload command receiver is a transistorized 
double-conversion, phase-lock-loop receiver which pro- 
duces a coherent output at 2/17 of the received frequency. 
It can be operated with either a 250-cps or a 40-cps band- 
width. The receiver operates continuously and, since its 
band-width is considerably less than the frequency uncer- 
tainty of the received signal, it repeatedly sweeps over a 
range of 30 kc searching for a carrier. The sweep period is 
10 seconds for the wide band and three minutes for the nar- 
row band. When the receiver acquires and locks on a signal 
from the earth the sweeping stops and the receiver can then 
accept any of 30 possible commands. A list of the 14 com- 
mands employed in the Able-3 payload is given in Table I. 


TABLE I 
COMMANDS TO PAYLOAD 


1) Digital Telemetry, 64 pps 
2) Digital Telemetry, 8 pps 
3) Television On 
4) 378-mc Transmitter Off 
5) 5-watt Transmitter On 
6) Accelerometer On 
7) 108-mc Transmitter Off 
8) 108-mc Transmitter On 
9) Over-Voltage 
10) Receiver Wide Band 
11) Receiver Narrow Band 
12) Clear Rocket 
13) Arm Injection Rocket 
14) Fire Rocket 


Note: When command 4 is given, command 6 is negated and the 
digital telemetry is restored to 1 pps. 


The digital command demodulator was developed after 
an extensive examination of various time and frequency 
multiplexed command systems. Investigation indicated that 
excessive weight and/or unreliability characterized each of 
the multiple-tone command systems capable of meeting the 
required bandwidth and frequency stability. The develop- 
ment of the digital decoder, however, required only one 
narrow-band filter and made use of many of the com- 
ponents already existing in the digital telemetry system. 

In the narrow-band position, the receiver exhibits a 
sensitivity (based on a 12-db noise figure) of —140 dbm. 
At the present state of development, commands can be 
given down to a sensitivity of only —133 dbm. 

Signals from the earth to the payload are transmitted by 
using a high-power carrier, phase-modulated with a 512-cps 
subcarrier. Amplitude modulation of the subcarrier using 
a train of 13 pulses permits the transmission of a syn- 
chronizing pulse, a simple address code, and a 5-bit com- 
mand together with its complement. 

The receiver weighs four pounds, occupies 1300 cubic 
inches, and draws 1.5 watts at 16 volts. The command de- 
coder used with the Doppler-command receiver is shown in 
Fig. 2. 

In its role as a transponder the payload receiver accepts 
a CW signal from the ground and after suitable process- 
ing delivers it to the 5-watt transmitter for retransmission 
to the ground at 16/17 of the received frequency. 
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Fig. 2—Command decoder. 


Range rate is measured on the ground to accuracies of 
one foot per second by extracting the Doppler frequency 
shift between the transmitted and received signal after 
correcting for the frequency offset introduced by the ve- 
hicle transponder. Range is measured by frequency modu- 
lating the ground transmitter by low-frequency sinusoids. A 
modulating frequency of 2 cps provide a range resolution of 
better than 500 miles when adequate smoothing exists. A 
second modulating frequency of 0.25 cps is used to resolve 
the phase ambiguities in the 2-cps signal. Further ambigui- 
ties are resolved by measuring the time delay between the 
start of modulation on the ground and the beginning of 
modulation in the return signal. Tracking in angle is per- 
formed by nodding the ground antenna alternately in ele- 
vation and azimuth. Angular accuracies of about 0.2° are 
possible. 

TELEMETRY 


The analog telemetry system carried aboard the Able-3 
payload is essentially a double-scale version of the system 
carried aboard the Pioneer I and II payloads. Eleven input 
signals are applied to two groups of subcarrier oscillators 
occupying the lowest five and the lowest six IRIG chan- 
nels respectively. A twelfth signal, a VLF propagation ex- 
periment occurs directly as a subcarrier. Two multiplexers 
sum the two groups of subcarrier oscillator outputs and 
the VLF experiment and deliver the resultant composite 
signals to two low-power VHF transmitters. One trans- 
mitter delivers a signal whose frequency is 108.06 mc, 
while the other delivers a signal whose frequency is 108.09 
me, 

The digital telemetry system. is carried aboard the Able-3 
payload principally to evaluate its usefulness in interplane- 
tary probes. However, the automatic data handling tech- 
niques which are available for digital information give it a 
definite purpose for use in an earth satellite. It is antici- 
pated that digital techniques will prove to be better suited 
to the transmission of information over interplanetary 
distances than are analog techniques because of the nature 


1959 


of the thresholding characteristics of the analog system.? 

The digital telemetry unit developed by STL, the Telebit 
system, is shown in block diagram in Fig. 3. A photograph 
of the unit is shown in Fig. 4. This system accepts both 
analog and digital inputs from various experiments. The 
converted information at its output appears as a binary 
coded subcarrier (1024 cps) which then phase modulates 
the signal of the 5-watt transmitter. 

The binary output of the Telebit system occurs at a 
synchronous rate and is composed of repeating sets of 
frames of words. For Able-3, 11 words per frame are 
used. One word of each frame is used as a frame sync 
and is read out as all zeros, while the balance of the words 
are coded with the digital representation of the input in- 
formation. Each word contains 12 pulses. The first two 
pulses (for information words) are always coded the same 
(zero, one) and define the start of a word; i.e., these two 
pulses provide a word sync. The other 10 pulses may take 
on any combination of binary values to represent a number 
from zero to 1023. 

Another advantage of the digital telemetry system over 
the analog system is that the rate of transmission of infor- 
mation can be adjusted as a function of range. Provision 
is made for the transmission of information at the rate of 
one, eight, or 64 pulses per second. This flexibility of in- 
formation rate will be found to have its greatest usefulness 
during the transit of interplanetary distances where large 
changes in range are encountered. At 55 million miles, 
using 150 watts of power, information can be conveyed 
only at one pulse per second; while at less than 500,000 
miles, five watts are sufficient to transmit at 64 pulses per 
second. 

A 12-bit combination binary counter and shift register, 
referred to as a shifting accumulator, is provided for each 
word. Pulses from a digital-type experiment are applied 
directly to the counting input of a shifting accumulator 
while an analog input is applied to an analog-to-digital con- 
verter whose output is then applied to a shifting accumula- 
tor. An electronic commutator running synchronously at 
the word rate gates 12 shift pulses to each shifting ac- 
cumulator during one word interval each frame. 

These shift pulses cause the information in the shift- 
ing accumulator to be delivered to the biphase modulator 
and at the same time the output of the digital shifting ac- 
cumulators is returned to the input so that after 12 shift 
pulses the state of the shifting accumulators is exactly as it 
began. The outputs of all of the shifting accumulators are 
connected together but since only one is shifting at a time 
no interference results. 

The conversion of analog to digital information is ac- 
complished with the aid of a digital ramp and a voltage 
comparison circuit. In essence conversion results from 
counting the number of steps in the ramp below the level 
of the analog input. The counting is actually done in a 


1G. E, Mueller and J. E. Taber, “An interplanetary communi- 
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Fig. 3—Block diagram of telebit system. 
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Fig. 4-STL telebit unit. 


shifting accumulator just as for digital experiments. The 
Able-3 payload has six digital and four analog shifting 
accumulators. 

The biphase modulator accepts the pulses emerging 
sequentially from the shifting accumulators and produces a 
subcarrier whose phase shifts by 180° each time there 
is a “one” to be transmitted and retains its phase each 
time a “zero” is to be transmitted. This biphase modu- 
lated subcarrier is then delivered to the transmitter for 
phase modulation upon the carrier. 

The pulses which cause the electronic commutator to 
step and the shift registers to shift originate in the pro- 
grammer, and the programmer in turn receives its excita- 
tion from the master clock. The programmer is equipped 
such that application of an outside signal derived from the 
digital command decoder causes the pulse rate of the digital 
telemetry system to change. 
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ANTENNAS 


The Able-3 payload uses both VHF and UHF for data 
transmission and UHF for guidance functions. The VHF 
transmission circuit from the two transmitters consists of 
a coaxial ring hybrid for coupling both transmitters to the 
single antenna system, a shunt tuning stub to match the 
antenna system to the hybrid, a diplexer for diverting sig- 
nals for the VLF whistler mode experiment, and two mono- 
pole antennas. 

The VHF diplexer consists of a resonant circuit tuned 
to 15 kc for the whistler experiment and is shunted by a 
band-pass filter employing a helical transmission line. The 
VHF monopole antennas are conventional in design and 
operate against the payload shell as the ground plane. Al- 
though the resultant radiation pattern is perturbed by the 
nonplanar ground plane and the solar cell paddles, pattern 
coverage and gain are adequate for the mission require- 
ments. The UHF antennas, similar in design to the VHF 
radiators but appropriately scaled in length, are two mono- 
poles operating against the payload shell located on the op- 
posite end of the payload from the VHF antennas. Both the 
antenna systems are symmetrically located with respect to 
the spin axis of the payload. 


PoweER SUPPLY 


Placing solar cells directly on the surface of the payload 
presented two difficulties: 1) the surface area of the pay- 
load was too small to permit adequate solar energy inter- 
ception; and 2) inasmuch as the solar cells also absorb 
large quantities of heat, the payload temperature would be 
intolerably high. An early decision was made to locate the 
solar cells on paddles extending essentially radially from 
the payload. These paddles fold within the nose fairing 
during launch but when fairing is jettisoned they spring 
out and latch into place. The paddles extend radially 22.5° 
from the payload equatorial plane as shown in Fig. 5. Ap- 
proximately 8000 silicon solar cells are carried on the 
Able-3 satellite. Because of attitude and spin considera- 
tions, however, only about 2000 of the solar cells are re- 
ceiving solar energy at one time. A prototype paddle with 
solar cells mounted on it is shown in Fig. 6. 

The payload storage battery possesses a 50-watt-hour 
capacity at 16.8 volts. A type of battery which would have 
been lighter for the same capacity, such as silver zinc, was 
not chosen because of doubtful reliability. Part of the 
power developed in the solar cells is immediately consumed 
in experiments and the receiver, and the charging rate of 
the battery is thus lessened by this amount. It is anticipated 
that between six and ten hours will be required to com- 
pletely recharge the battery. 

An over-voltage control operated on command from the 
ground can be used to reduce the charging rate if it appears 
that the battery is being consistently overcharged. An 
under-voltage control is provided which automatically re- 
moves the 5-watt transmitter from the load on the batteries 
in case battery discharge goes so far as to disable the re- 
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Fig. 6—Solar cells mounted on paddle. 


ceiver and thus prevent commanding off the transmitter 
from the ground. 

The batteries operate at 16.8 volts. Since this voltage is 
inadequate for all of the electronics, a series of static con- 
verters is employed to provide a variety of voltage levels. 


TEMPERATURE CONTROL 


Temperature control of the payload is achieved by ar- 
ranging the thermal balance between input solar radiation 
and infrared emission so that it occurs at a temperature 
within the desired operating range. This is done by select- 


1959 


ing a surface coating that will give a value of the ratio of 
solar absorptivity to long wavelength emissivity (#/¢) in 
the range of about 1.0 to 1.3. 

A similar problem exists in the solar cell paddles but 
there is less freedom of solution. Furthermore, the solar 
cells are quite temperature-sensitive and their conversion 
efficiency drops off about 0.6 per cent per degree Centigrade 
above a nominal value of 25°C. 

Because the absorptivity-to-emissivity ratio of a silicon 
cell is about 3, the cells will run excessively hot unless their 
surfaces are coated to alter the long wavelength emissivity. 
Thus glass plates, 0.0003-inch thick, cemented individually 
to each cell, are used to increase the long wavelength emis- 
sivity without appreciably altering the cells’ conversion effi- 
ciency. To protect the cement from ultraviolet radiation 
and to reduce the cell temperature even more by eliminat- 
ing a part of the solar spectrum to which the solar cells are 
functionally insensitive, an ultraviolet reflective coating is 
applied to the glass before it is cemented to the cell. For 
the glass plates to reduce the temperature successfully, it 
is essential that the glass be bonded to the cell over the 
whole surface so that there is no gap to cause a greenhouse 
effect. 

The large thermal capacity of the payload prevents large 
drops in its internal temperature for eclipses expected 
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during the first 1000 days of flight. Temperature drops of 
15 to 20°F are the most severe that are expected during 
this time. After about 1000 days, however, long eclipses 
(about 2%4 hours duration) will occur; the payload tem- 
perature drop by the end of such eclipse will be the order of 
DOeu 

Temperatures of the solar cells, on the other hand, drop 
rapidly during an eclipse. While there is, of course, no 
operating requirement for the solar cells during an eclipse, 
they must be able to survive the low temperatures. For a 
one-half hour eclipse, the lowest temperature on the solar 
cell paddle should be about —50°F, and about —175 to 
—200°F for a 24-hour eclipse. 


CoNCLUSION 


Considerable effort has gone into the solutions briefly 
outlined here. The component and system designs repre- 
sented by the payload are felt to be sound and laboratory 
environmental testing has indicated their reliability. How- 
ever, long-term operation in space can provide the only 
adequate test and the results of this test are bound to 
indicate areas requiring further development. In any event, 
{future space probe vehicles necessitating considerably more 
sophisticated instrumentation at considerably further 
ranges will be able to profit from this experience. 


Three-Dimensional Interplanetary Trajectories* 
A. B. MICKELWAIT?, E. H. TOMPKINS, JR.t, anp R. A. PARK} 


Summary—Developments in guided missiles and related sci- 
entific areas have reached a state where ballistic flights to the 
planets Venus and Mars are feasible and imminent with existing 
hardware. This paper discusses considerations involved in making 
actual ballistic flight to Venus. Conventional two-dimensional 
analysis is compared to the actual three-dimensional case. Burn- 
out conditions required to minimize the problems of guidance, 
payload, range safety, booster aerodynamics, midcourse guidance, 
and capture are established. Problems relating to launch times, 
launch latitude, and interplanetary communications are also 
discussed. 


INTRODUCTION 
\ LL significant design criteria for interplanetary ve- 


hicles come directly from the interaction between 
the gross properties of interplanetary ballistic 
trajectories and their sensitivities to small initial dis- 
turbances. In other words, the dynamic and geometric 
problems must be considered simultaneously with the 
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accuracy requirements. In fact, the two subjects are so 
interrelated that propulsion and guidance requirements, 
the engineering counterparts of orbital analysis, must 
be regarded as inseparable in interplanetary vehicle de- 
sign in order to adequately resolve the often conflicting 
demands of payload needs and high accuracy. 

It is commonplace to consider the effects of the rela- 
tive orientation of the earth’s rotational axis, the earth’s 
orbital plane, and the orbital plane of the other planets 
as niceties that can be relegated to the final high-preci- 
sion orbit determination. For earth-bound satellites this 
is largely true. Indeed, although for lunar missions 
three-dimensional considerations cannot be neglected, 
they do not usually contribute decisive vehicle design 
factors. However, even a preliminary analysis of inter- 
planetary flights indicates that they can miss their goals 
by an order of magnitude, from the point of view of 
both payload and accuracy, if an oversimplified two- 
dimensional analysis is relied upon. 

Because Venus has an almost optimum orientation in 
its orbital plane with respect to the earth’s orbital plane 
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in 1959, the trajectory design problem for a flight to 
Venus in 1959 will be analyzed in the most detail. The 
characteristics of this trajectory will be contrasted with 
the problems of flights to Mars in 1960 and Venus in 
early 1961. To keep the paper a reasonable length, pri- 
mary attention will be given to purely ballistic flights, 
that is, flights with no orbit corrections after leaving 
the earth’s vicinity. Furthermore, the paper will be 
limited to the range of propulsion and guidance capabil- 
ities available at the present time. However, it will be- 
come apparent that no interplanetary mission even with 
the crudest requirements would have a reasonable 
probability of success unless some midcourse or terminal 
guidance is carried out; therefore some trajectory char- 
acteristics directly pertinent to these problems will be 
introduced. 


Two-DIMENSIONAL KINEMATICS 


First, we will assume that the earth and target planet 
move in circular orbits in the same plane and that the 
earth’s equator is in this plane. Let us assume further 
that the vehicle is launched and flies in this common 
equatorial plane. To handle the analytically intractable 
multiforce problem posed by the attraction of the sun 
and all the planets simultaneously, we will assume that 
the vehicle’s orbit consists of three phases (see Fig. 1): 
first, escape hyperbola centered at and determined 
solely by the earth’s gravitational field; second, the sun- 
centered ellipse controlled only by sun’s gravity and 
using injection conditions taken from the escape hyper- 
bola leaving the earth; and third, a reverse of the first 
phase, the approach hyperbola centered at the target 
planet and controlled entirely by its gravity. 

When leaving the earth, the vehicle is launched from 
a platform (the earth) moving about 100,000 feet per 
second around the sun. For most practical trajectories 
to Venus or Mars, the vehicle leaves the earth with a 
velocity of about 36,000 to 40,000 feet per second rela- 
tive to the earth. By the time it has left essentially all of 
the earth’s gravity along the escape hyperbola relative 
to the earth, it either adds (going to Mars) or subtracts 
(going to Venus) only about 10,000 feet per second to 
the earth’s much larger orbital velocity. Although, in 
principle, it is possible that a significant spatial error at 
injection into the sun’s field could occur (that is, in two 
dimensions, an error in the radial distance from the sun) 
other constraints upon launch and burnout conditions, 
which will be discussed later, make this possibility 
negligible. Therefore, this source of error will not be 
discussed. The guidance problem is then simply that of 
achieving those conditions at burnout which will give 
the vehicle the correct velocity when it enters the sun’s 
field. If the mission requirements are made sufficiently 
stringent—say, for example, that one wishes to impact 
at a particular point on Venus—position variables at in- 
jection into the sun’s field would need to be considered, 
but midcourse corrective maneuvers could be used to 
reduce this source of error. 
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Fig. 1—Two-dimensional coordinate systems for the three flight 
phases. (a) Escape hyperbola relative to the earth; (b) solar ellipse; 
(c) approach hyperbola relative to Venus. 


At this point it should be noted that our lack of 
knowledge of certain physical constants, especially the 
relevant dimensions of the solar system, will introduce 
their own uncertainties in position and velocity. How- 
ever, the most recent estimates of the error contributed 
by our lack of exact knowledge of the astronomical unit 
indicates that it can be safely ignored when compared 
with other sources of error, even if the present estimate 
of experimental uncertainty is doubled. 

Fig. 1(a)-1(c) shows the coordinate systems used and 
the relative orientation of earth and Venus at take-off. 
It is apparent from the geometrical relationships shown 
in this figure that the conditions in the sun’s frame of 
reference for a successful flight are given by 


6r(V «, Bs) = WyT( van Bs) ae Or(to) (1) 


1959 


when 


3 
T 
= 


where 
6;=the angle in the sun’s frame of reference turned 
by the vehicle between injection and impact 

®,=mean angular velocity of Venus 

T=total flight time 

6,=initial lead angle of earth relative to Venus 
r=vehicle distance from the sun 

R,= Venus’ mean orbital radius. 


Eq. (1) says that the properties of the trajectory after 
injection are determined entirely by the magnitude of 
the velocity at injection into the sun’s frame, | V,|, and 
the angle between V, and the radius trom the sun, Bs. 
It follows that if the launch day is fixed, that is, if @z is 
fixed, then there exists a functional relation between 
V,| and @,;. Referring to Fig. 1(b), we see that this 
relation may not be one to one; that is, there may be 
two sets of V,| and 8, allowable on a given launch day, 
one set (Class I orbits) corresponding to an impact 
when r=R, as the vehicle crosses the Venus orbit the 
first time and the other set (Class II orbits) to that 
trajectory when the vehicle impacts as it crosses the 
Venus orbit the second time. 

The classification of orbits into classes is useful for 
several reasons. One important system consideration, 
communication transmission distance, depends strongly 
on the orbit class. Class II orbits generally involve 
longer times of flight than the minimum energy orbits 
(~147 days to Venus and 260 days to Mars). Since the 
target planet is moving away from the earth at the end 
of a minimum energy trajectory at the rate of about 10 
nautical miles per second, relatively small changes in 
flight time can radically increase or decrease the trans- 
mission distance. In fact, only several hundred feet per 
second change in take-off velocity can vary the trans- 
mission distance by as much as 30 million miles. 

It is clear that the most efficient use of V., is made 
when it is collinear with the velocity vector of the earth, 
V., and when 6,=90°. If the magnitude of V,, is fixed 
and its direction changed, V, rotates so that 8, does not 
equal 90° and a larger | V,| results (Fig. 2). A larger 
| V. implies, of course, a shallower penetration into 
Venus orbit. When 6,=90°, reducing | Vin | increases 
| V. until the two points where the vehicle crosses 
Venus’ orbit coalesce. The vehicle’s trajectory is then 
tangent to Venus’ orbit at r= R,, and the vehicle travels 
exactly 180° from launch to impact as measured relative 
to the sun (the “Hohmann” transfer ellipse). A smaller 
| Vo | and therefore a smaller take-off energy will give 
an orbit around the sun that falls inside of the earth’s 
orbit but does not touch Venus’ orbit. 

Recalling that (1) related 6, | V.| , and B,, we see that 
if B, and | V,| are adjusted for a minimum energy orbit, 
then a unique launch day is determined. For other 
values of 6, (varying roughly one-half degree per day for 
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Fig. 2—Relation of the earth’s velocity, Ve, and vehicle’s velocity 
with respect to the earth after escape, V., to vehicle velocity in 
the sun’s frame, Vs. 


the earth and Venus) and for the most efficient 6B, =90°. 
smaller values of | V,| will be required since the two points 
of intersection must move away from 6; = 180°, implying: 
a deeper penetration inside of Venus’ orbit than the 
Hohmann ellipse and, therefore, greater reduction of 
the original | Ve 

Fig. 3 shows the increase in velocity requirements 
over and above that required for the two-dimensional 
Hohmann ellipse as the angle 9; is increased or de- 
creased away from 54°. The system implication of the 
velocity increase is that there is a definite period of days 
during which launching to another planet is feasible 
with limited propulsion capabilities. On the other hand, 
since the angular velocities of Venus, earth, and Mars 
are quite different,! the required relative angular orien- 
tation will reoccur periodically with the so-called 
“synodic” period; for Venus every 584 days, for Mars 
every 780 days. As we shall see, however, when the in- 
clination of the various orbit planes are taken into ac- 
count, the time between reoccurrence of identical 
launching conditions becomes much larger. The key 
dates to be kept in mind for planetary probes are June 
7, 1959 to Venus, October 11, 1960 to Mars, and Janu- 
ary 7, 1961 to Venus again. On these days, at least with 
a two-dimensional analysis, the minimum energy 6, 
occurs. 

In addition to conserving energy, the minimum energy 
trajectory is also desirable for another reason. Most mis- 
sions attempt either to maximize the length of time 
spent in the target planet’s vicinity or enter the target 
planet’s atmosphere as slowly as possible to diminish 
deceleration and heating. In any case, the slowest pos- 
sible relative approach velocity is desirable. From Fig. 
1(c), and recalling that the inertial velocities of target 
and vehicle are of the order of 100,000 feet per second, 


. 


1 Earth, w7=0.986°/mean solar day, 
Venus, wy = 1.600°/mean solar day, 
Mars, wy =0.524°/mean solar day. 
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Fig. 3—Burnout velocity in excess of minimum as a 
function of lead angle, @z. 


we see that even a small inclination of the vehicle’s 
velocity, Vven, to that of the target planet, V», can 
greatly increase the relative velocity, Vre1. Therefore, 
the tangential approach afforded by the two-dimen- 
sional Hohmann ellipse provides the desired minimum 
approach velocity. 

Considerations concerning the approach phase of the 
trajectory may be made more quantitative by relating 
the injection conditions into the approach hyperbola to 
the distance of closest approach. By utilizing one-body 
formulas expressing conservation of energy and angular 
momentum relative to the target planet, we arrive at 
a relationship between the impact parameter, 0}, ap- 
proach velocity, Vre, and distance of closest approach, 


Dp 
eae ie a 2GM, (2) 
4 Dee 


The impact parameter is simply the distance the vehicle 
would miss the target in the absence of any attraction 
by the target planet. The quantity (GM,) is the uni- 
versal gravitational constant times the mass of the tar- 
get planet. For an impact, D would, of course, be set 
equal to the target radius.’ For the lowest possible ap- 
proach velocity to Venus, ~9000 feet per second the 
maximum allowable impact parameter is about twice 
D; that is, the gravity field of Venus has effectively 
doubled the geometrical size. Clearly, then, from this 
standpoint alone it is desirable to keep Vre as low as 
possible. 


Two-DIMENSIONAL GUIDANCE ANALYSIS 


The first approach to studying sensitivity (of the 
trajectory to errors at burnout) is to look at the sun- 
frame ellipse by itself. From (1) we see the expressions 


2 Radius of Venus 3410 naut mi, 
radius of Mars=1820 naut mi. 
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give the distance from the target planet when the ve- 
hicle crosses the target orbit due to injection errors AV, 
and AB,. However, this is not the actual miss distance 
due to injection dispersions since the vehicle is usually 
closer to the target either before or after crossing the 
target orbit. Either the variation of distance of closest 
approach, D, or variation of impact parameter, b, are 
more suitable measures of orbit sensitivity. The impact 
parameter is often used extensively since its sensitivity 
is easier to calculate. If the approach velocity is kept 
within a small range, the maximum allowable dispersion 
in b for the success of a given mission may be estab- 


lished. Then 
2 0b Dea yl 
AV; AB; 3 
y+ ae) | © 


0b 
nai 
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determines the allowable dispersions in 8; and V, over 
this range of approach velocities. Eq. (3) assumes im- 
plicitly that AV, and A@, are random, independently 
distributed quantities. In Figs. 4 and 5, we see the 
variation of the impact parameter miss coefficients 
0b/dV; and 06/08, over a range of injection conditions. 
The important behavior to note here is the vanishing 
of B; sensitivity at minimum energy and the vanishing 
of V, sensitivity at velocities slightly above minimum 
energy. Although not shown in Figs. 4 and 5 (the curves 
are for Class I orbits only), Class II orbits do not pos- 
sess the zero in velocity sensitivity. 

The two-dimensional analysis presented so far is de- 
ficient for two important reasons. First, the effect of 
escape phase has not been included in the analysis. Re- 
ferring to Fig. 1, we see that a change in B» at burnout 
varies the orientation of V,, although not its magnitude. 
But changing the orientation of V., changes both 8; and 
| V,|. Likewise, a change in | Vol will change | V..| and 
therefore change both 8, and | V,| simultaneously. Only 
when V., becomes collinear with V, do By and Vo interact 
directly with their counterparts in the sun’s frame of 
reference. However, as will be shown in the next section, 
when the complete three-dimensional situation is taken 
into account, it is very difficult to have V,, parallel to 
V. except for a unique orientation of the earth relative 
to Venus. Furthermore, when the third degree of free- 
dom is introduced, the relations between burnout condi- 
tions and injection conditions become more complex. 
The net result is that the strong mimina in sensitivity 
are flattened out and in fact will vanish altogether ex- 
cept for isolated situations. In the light of the failure of 
the two-dimensional analysis to carry over into three 
dimensions, it is not fruitful here to carry out the com- 
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Fig. 4—Two-dimensional velocity miss coefficient, 0b/dVs, as a func- 
tion of the vehicle velocity in the sun’s frame, V,,!for three flight 
path angles in the sun’s frame, Bs. 
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Fig. 5—Two-dimensional angular miss coefficient, 0b/d85, as a func- 
tion of vehicle velocity in the sun’s frame, Vs, for three flight path 
angles in the sun’s frame, Bs. 


plete analytical treatment of connecting burnout dis- 
persions to solar injection dispersions and then to im- 
pact dispersions; rather we will present in the next sec- 
tion some results of the exact computer analysis of orbit 
sensitivities.® 


THREE-DIMENSIONAL CONSIDERATIONS 


Although the orbital planes of the earth and Venus 
are inclined to each other by only a small angle (~3.4°), 
the three-dimensional aspects of their relationship has 
a strong influence on the energy, guidance, and target 


3 A complete two-dimensional analysis is presented in T. A. Mag- 
ness, J. B. McGuire, and O. K. Smith, “Accuracy Requirements for 
Interplanetary Ballistic Trajectories,” Space Technology Labs., 
Rept. TR-59-0000-00682, May 18, 1959; presented at the Ninth 
Annual Congress of the International Astronautical Federation, 
Amsterdam, The Netherlands, August 25-30, 1958. 
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Fig. 6—Three-dimensional coordinate system in the sun’s frame. 
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Fig. 7—Three-dimensional aspects of the orbits of Venus, 
earth, and the probe. 


approach of the free-flight trajectory. We shall consider 
first the energy relations for the three-dimensional 
situation. 


A. Energy Considerations 


Let a probe be launched with a burnout velocity, Vo, 
in a coordinate frame, S,, fixed at the center of the earth 
(but not rotating with the earth). As the probe leaves 
the vicinity of the earth it travels along a hyperbolic 
escape orbit with its velocity vector gradually approach- 
ing the asymptote of the hyperbola, and its magnitude 
gradually approaching a constant value, V.,. As before, 
the velocity of the earth with respect to the sun is V.. 
The vector sum of V., and V, is again V;, which is the 
velocity of the probe in the sun’s reference. To a close 
approximation, the probe’s orbit will lie in a plane 
through V, and the center of the sun. Finally, the angle 
between the ecliptic (the plane of the earth’s orbit) and 
the probe’s orbit plane is 7, and the angle between V; 
and a radius vector from the probe to the sun is B,, as 
shown in Fig. 6. The quantities V;, 7;, and 8, are the 
three-dimensional “injection conditions” for the probe 
orbit about the sun. 

Fig. 7 shows the inclined orbital planes of the earth 
and Venus intersecting in the line AB (the line of 
“nodes”). Imagine first that the inclination is set equal 
to zero so that both orbits lie in the same plane. This case 
is identical to the two-dimensional case described in the 
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Fig. 8—Burnout velocity, Vo, as a function of the in-plane angle, 6, 
for various flight path angles in the sun’s frame, 8, (Venus flight’ 
1961 era). 


previous sections, and the Hohmann transfer ellipse to 
Venus will occur when 6z=54°. For the real solar sys- 
tem, the inclination of Venus’ plane is not zero, and in 
each synodic period the two-dimensional minimum 
energy lead angle of 0; =54° will occur at varying places 
on the earth’s orbit. Suppose 6,=54° when the earth is 
in the position shown in Fig. 7. If a probe is launched 
then with 7,=0°, that is, in the ecliptic, it can reach 
Venus only at a point such as the ascending node, Q, 
where Venus crosses the ecliptic plane. However, for 
this trajectory 6; is much less than 180°, implying a 
quicker trip and expenditure of more energy (by chang- 
ing V, or B,) than in the two-dimensional case. One 
might take a longer trip and reach Venus at the de- 
scending node, @, but this would also require more 
energy. Consider now the situation when 7, is changed 
from zero. In the two-dimensional case, the minimum 
energy flight traveled 180° around the sun. Study of 
Fig. 7 will show that a value of 6; = 180° in three dimen- 
sions is possible only for the case when the trajectory 
travels almost over the ecliptic pole (7,~90°), a maneu- 
ver which would require an extremely large burnout 
velocity (> 100,000 feet/sec). Let us look at some of the 
details on realistic energy requirements. 

Fig. 8 shows the burnout velocity, Vo, on earth as a 
function of the in-plane angle, 07, for constant values of 
B, for flights to Venus in the 1960-1961 era. Each point 
on the curves corresponds to a certain | V,, i,} set and 
thus to a certain launch date and flight time. In this 
figure, the following orbit classification is used: 

Type I trajectory—has a flight angle 67<180°, 

Type II trajectory—has a flight angle 6;> 180°, 
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Fig. 9—Burnout velocity, Vo, as a function of launch dates for various 
flight path angles in the sun’s frame, @; (Venus flight, 1961 era) 
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Fig. 10—Detailed analysis of one curve, flight path angle, 8, =88°, of 
burnout velocity, Vo, as a function of launch date (from Fig. 9). 
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reaches Venus, D, (in millions of nautical miles) are given at 
various points on the curves. 


Class I trajectory—arrives at Venus on first crossing 
Venus’ orbit, 
Class II trajectory—arrives at Venus on the second 
crossing of Venus’ orbit. 

Fig. 9 shows Vo asa function of launch date for constant 
values of 8;. From Figs. 8 and 9, it can be seen that the 
three-dimensional minimum energy flight to Venus in 
the 1960-1961 era occurs about January 15, 1961 (over 
a month after the earth crosses Venus’ line of nodes) 
with 6,=90° and 6;=150°. For this flight, Vo= 36,920 
feet per second (about 300 feet per second above the 
two-dimensional absolute minimum velocity), 7;=1.1°, 
V,=90,000 feet per second, and flight time t; = 126 days. 
The probe orbit makes an angle of ~4.3° with the 
Venus plane, and, on the launch date, the earth leads 
Venus by ~36° rather than the 54° predicted by two- 
dimensional analysis. Figs. 8 and 9 also show, as was 
concluded from Fig. 7, that considerable additional 
velocity in excess of the minimum must be expended to 
obtain an in-plane angle near 180°. Fig. 10 shows in de- 
tail the parameters corresponding to several points on 
the B, = 88° curve of Fig. 9. Considering such points, one 
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Fig. 11—Burnout velocity, Vo, as a function of launch dates fot 
oa flight path angles in the sun’s frame, 8, (Venus flight, 1959 
era). 


may begin a systems analysis of a Venus flight by ex- 
amining such factors as the following: 

1) Each operating point on the curves determines a 
particular out-of-plane angle, 77, and in-plane 
angle (through 67) at the crossing of Venus’ orbit. 
These angles have a strong effect on the guidance 
sensitivity of the trajectory and hence on the 
probability of success of the flight. They also affect 
the relative approach velocity in Venus’ frame and 
hence the capture, impact, or landing require- 
ments. 

2) Similarly, each point implies a certain transmission 
distance at the time of arrival of the probe. The 
fact that this transmission distance may vary from 
23 to 138 million nautical miles will have a strong 
influence on communication capability in the early 
years of interplanetary flight and will, in large 
measure, determine the actual trajectory flown. 

3) Each point implies a particular set of burnout 
angles for a given launch site. As will be shown in 
a later section, these angles may be quite unac- 
ceptable either because they affect the efficient 
operation of the powered flight part of the trajec- 
tory or because they may result in the booster 
stages impacting in inhabited areas. 

In contrast to Fig. 9, Fig. 11 shows Vo as a function of 
launch date with a constant 8, for Venus flights in the 
1959 era. It should be noticed that the 1959 flights pre- 
senta very special set of characteristics resulting from 
the fact that the three-dimensional minimum energy 
trajectory occurs very close to the time of the earth’s 
crossing of Venus’ line of nodes; which is to say that the 
three-dimensional trajectory is almost identical with the 
two-dimensional trajectory. It was therefore possible to 
fly with near absolute minimum velocities on essentially 
two-dimensional trajectories which entered Venus’ 
plane at the descending node-crossing and arrived nearly 
tangential to Venus’ orbit at the ascending node. This 
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Fig. 12—Coordinate system at burnout in the earth’s frame. 


coincidence will occur every eight years, but there will 
be a gradual secular change away from the exact two- 
dimensional conditions. Once the conditions become a 
few weeks apart, the two-dimensional analogy is no 
longer relevant; therefore, the 1959 case should be re- 
garded as essentially illustrative. 

One should also note that by increasing burnout 
velocity only a relatively small amount above the 
minimum, flights to Venus in each era are possible over 
a period of several months and that a wide variety of 
types of trajectories is possible during each period. 
These problems related to energy considerations are 
only one factor to be balanced against others in the 
complete systems analysis. 


B. The Three-Dimensional Escape Trajectory 


We shall now define the variables used to specify 
completely the burnout conditions of the powered flight 
on earth and show how these variables are constrained 
by the conditions V,, 8., and 7, of the sun-frame transfer 
orbit. 

It is convenient to employ the following set of initial 
conditions in studying the powered flight, specified at a 
burnout time fo (to =date-+time of day): 

ao =right ascension‘ of burnout point 

5)=declination® of burnout point 
vo = geocentric radius of burnout point 

Vo=magnitude of inertial burnout velocity 

8)=angle between inertial velocity vector, Vo, and 

geocentric radial line (“velocity path angle”) 

Ay=azimuth of Vp from North, measured in a plane 

perpendicular to the geocentric radius at the 
burnout point. 

These definitions are illustrated in Fig. 12. In reality, 
the escape hyperbola will blend into the elliptical orbit 
around the sun and no clear-cut place will exist which 
one can speak of as the injection point into a central 
force orbit about the sun. Thus, one must ultimately 


4 Right ascension is the astronomical measure of longitude from 
the vernal equinox, which is the ascending node of the sun in its 
apparent motion around the earth. 

5 Declination is the astronomical measure of latitude from the 
equatorial plane. 
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ous burnout velocities, Vo (Venus flight, June 6, 1959, Cape 
Canaveral launch, burnout altitude 22,000,000 feet from center of 
the earth). 


determine a, 50, 70, Vo, Bo, and Ao at ¢o as initial condi- 
tions of a system of differential equations of the probe’s 
motion such that the probe reaches a specified point on 
or near Venus. These equations must include the effects 
of the earth’s oblateness and the attraction of the moon, 
Venus, etc. We shall now follow this approach. The sun- 
frame conditions V;, 8;, and z,, will remain useful guides, 
however. 

Let us hold 79 fixed at a reference altitude of 22,000,000 
feet (~177 nautical miles above the earth’s surface). 
The powered flight will normally lie approximately in a 
plane through the launch site at lift-off and the burnout 
point at to. To conserve energy, it is desirable that the 
free-flight continue in this same plane. Therefore, 69 and 
the longitude of the burnout point are determined once 
the range of the powered flight arc over the earth, the 
location of the launch site, and the azimuth, A,, are 
chosen. From the longitude of the burnout point and 
the burnout time, fo, the right ascension, a@, may be 
computed. Three parameters, azimuth, Ao, flight path 
angle, Bo, and burnout velocity, Vo, remain to be chosen. 
Fig. 13 shows Ao as a function of 8» with curves of con- 
stant burnout velocity, Vo, fora flight to Venus launched 
from Cape Canaveral, Fla. on June 6, 1959. The curves 
were obtained by integration of the full equations of 
motion in a computer program which iterated the initial 
conditions to obtain a probe impact at Venus. Different 
points on the constant Vo curve correspond to different 
a values, and hence to different burnout (and launch) 
times.® The curves of constant Vo correspond roughly to 
curves of constant 1%,. 


6 The right ascension of a fixed point on the earth’s surface varies 
as the earth rotates since it is referred to the vernal equinox which is 
a fixed line in inertial space. 
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For a given mission and launch site, the range of 
“possible” flight azimuths is tied closely to the sun- 
frame orbit and the planet geometry. The curves in 
Fig. 13 permit powered flight azimuths well within the 
bounds of the Cape Canaveral test range. However, the 
“possible” azimuths for any given mission may be at 
odds with range safety problems. But using “safe” azi- 
muths may have adverse effects upon the other system 
factors such as energy, guidance, communication dis- 
tance, etc. This problem can be clearly illustrated utiliz- 
ing the reasonable launch time to Mars of September 23, 
1960 when the earth is at the autumnal equinox, as 
shown in Fig. 12. 

A flight to Mars on this date can be made a simple 
geometric problem for illustration since 7, is equal to 0° 
and B,=90° at a reasonable velocity. To achieve these 
values for 7; and B,, we must make V,, and the asymp- 
tote of the escape hyperbola be parallel to the earth’s 
velocity vector. V, must also add to V., carrying the 
probe out towards Mars’ orbit. One way in which this 
may be accomplished is to fly due north (azimuth, A., 
equals 0°) from the launch site at the instant when Cape 
Canaveral passes through the plane containing the 
earth’s axis and V,. The escape geometry then lies in a 
plane as shown in Fig. 14. If the powered flight is as- 
sumed to cover a 15° arc over the earth’s surface, the 
escape hyperbola must turn through an in-plane angle, 
6; = 113°, to bring V,, parallel to V.. Analysis shows that 
a Class II flight to Mars with a Vo of ~37,600 feet per 
second at 7)9=22,000,000 feet can achieve 7,~0° and 
Bs-~90° (approximately the minimum three-dimensional 
Vo for a 1960 Mars flight). Given Vo, 79, and Bo, the 
angle 6, may be computed from simple two-body 
formulas. Conversely, if Vo, 70, and 6, are known, the 
burnout velocity path angle, Bo, is determined. Fig. 
15 shows 6, vs Vo for various values of 8) with 7 
= 22,000,000 feet. From this figure, we see that our 
Mars flight requires a By of ~72°. 

However, suppose it is necessary to increase azimuth 
A from 0° to, say, between 100° to 125° for reasons of 
range safety. This may be done by delaying launch until 
the launch site has rotated farther east which inclines 
the escape hyperbola plane (Fig. 16). However, this 
will decrease 6, and, as is seen in Fig. 15 fora fixed 
Vo, decrease Bo. This generates curves of By vs Ao 
similar to those of Fig. 13. In fact, however, for an 
efficient powered flight, Bis limited to a practical range 
of 65° to 90°. Thus, a practical limit on the azimuth 
which can be followed is reached when By goes below 
65°. For the case shown, it appears that the hypothetical 
range safety considerations make the flight impossible 
on the assumed launch date. However, there are several 
alternatives: 


1) Change the initial conditions at injection into the 
sun’s frame, and/or the launch date so that V, 
will have a more satisfactory direction. This will 
usually involve the expenditure of additional 
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Fig. 14——Geometry of escape hyperbola which lies in a plane 
including the earth’s velocity vector. 
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Fig. 15—In-plane angle of escape hyperbola, 6s, as a function of burn- 
out velocity, 8o, for various flight path angles at burnout, Ao. 
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Fig. 16—Rotation of escape hyperbola with change in launch time. 
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energy to increase Vo, and also this may produce 
trajectories which have undesirable guidance or 
communication properties. 

2) Launch at an earlier time when the trajectory 
plane through the launch site is inclined in a north- 
west-southeast orientation, launching in the south- 
east direction. This would also be incompatible 
with a reasonable Bo since 6, would be too large. 
Now, however, the probe may be placed in a satel- 
lite orbit and allowed to coast approximately half- 
way around the earth until the proper @, can be 
obtained after a final stage has been fired to inject 
into the escape hyperbola. 

3) Relocate the launch site. 


Some or all of these alternatives may be undesirable. 
This example merely emphasizes that launch azimuth 
considerations may have a very strong effect on the 
systems analysis of interplanetary flights in the early 
years of such flights.” 


C. Three-Dimensional Guidance Requirements 


After nominal interplanetary trajectories have been 
determined, the next step is to study the effect of dis- 
persions of the nominal produced by variations in the 
initial conditions (burnout conditions of the powered 
flight). Such a study will establish the requirements 
which must be met by a guidance system acting during 
the powered flight parts of the trajectory. But, as men- 
tioned earlier, the guidance requirements themselves 
will have a strong bearing on the choice of a nominal 
trajectory since specific “good guidance” trajectories 
exist which are less sensitive to errors by a factor of 10 
or more than neighboring trajectories. 

The terminal (or impact) portion of a three-dimen- 
sional interplanetary trajectory may be described in a 
manner analogous to the two-dimensional discussion. 
The approach hyperbola centered at the target planet 
may be characterized by the velocity relative to the 
target planet at a great distance from it, Vrei 7, and the 
impact parameter, 6. The magnitude of the impact 
parameter, now a vector quantity, is still defined as the 
distance of closest approach which would occur if the 
target has no gravity. (More precisely, 6 is the vector 
from the center of target planet that is perpendicular to 
the asymptote of the approach hyperbola.) 6 then shows 
the direction and is a measure of the distance of closest 
approach without gravity (not the actual miss with 
gravity) of the probe at the target planet. The sensitiv- 
ity of the miss to changes in the initial conditions (apo, 
50, Bo, Ao, Vo, 70) back on earth may be measured by de- 
fining the vector miss coefficients: 


0b 0b 
Cy Ve 


Bo dVo 


7 The problem will change from era to era; for example, azimuths 
required for the 1960-1961 Venus flights are much more in a northerly 
direction than for the 1959 Venus flights. 
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APPROACH ORBIT 


Fig. 17—Geometry of approach phase to target planet. 


where 6,, is the miss coefficient of the flight path angle, 
A,, of the azimuth and V,, of the velocity. Suppose, for 
example, a nominal trajectory has been obtained whose 
asymptote passes through zero (Fig. 17) providing a 
“dead-center” impact on the planet. A change of 1° in 
the 8) of the nominal will move the trajectory away 
from zero, say to the trajectory 73. 6, will be a direc- 
tional measure of the resulting miss in nautical miles per 
degree. 

Let us now return to the Bo vs Ao curves of Fig. 13 
and examine how the miss coefficients Vin, Bm, and An 
vary with the operating point on these curves and with 
launch date. Figs. 18 and 19 show curves of Vm, Bm, and 
Am as a function of Vo for two launch dates for Venus 
probes in June, 1959.° These velocity profiles correspond 
to Ay and {> values along a constant ap (~launch time) 
curve on the 8» vs Ao plots. Machine computations 
show that ordinarily the miss coefficients along a con- 
stant Vo curve vary only slightly as compared to the 
variations which occur with changes in Vo. Therefore, 
since points on a constant Vo curve correspond to 
changes in launch time, it follows that the miss coeffi- 
cients with constant Vo vary only slightly if launch time 
is varied by an hour or so. 

Figs. 18 and 19 indicate the guidance accuracy re- 
quired for interplanetary ballistic flights—angular error 
sensitivities of 40,000 to 500,000 nautical miles miss per 
degree and velocity sensitivities of 1000 to 20,000 nau- 
tical miles per feet per second being typical. The sensi- 
tivities vary strongly with the class of orbit, the launch 
velocity (~7,), and the relation of the launch date to the 
earth’s crossing of the target planet node. Furthermore, 
the “good guidance” points for Vo, Bo, and Ay will ordi- 
narily not occur at the same operating point on the By 
vs Ap curves so that optimization procedure must 
include a knowledge of the magnitudes of each of the 
errors at burnout that the guidance system cannot elim- 
inate. 


§ These miss coefficient curves were obtained with a computer 
program which solved a system of variational equations along with 
the normal equations of motion of the probe. 
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Fig. 18—Velocity, flight path angle, and azimuth miss coefficients 
Vin, Bm, and Am, as a function of burnout velocity, Vo, for Venus 
flight, June 6, 1959 (before the earth crosses line of nodes). 
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Fig. 19—Velocity, flight path angle, and azimuth miss coefficients, 
Vin, Bm, and Am, as a function of burnout velocity, Vo, for Venus 
flight, June 9, 1959 (after the earth crosses line of nodes). 
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Fig. 20—Residual miss components, azimuth perpendicular, A\, and 
flight path angle perpendicular, 8,, and their sum for 4° errors in 
Agand Bo. The parallel components have been removed at vernier 
correction (Venus flight, June 1959). 


These extreme burnout tolerances required for inter- 
planetary ballistic flights may be met if a very high de- 
gree of accuracy is attained at burnout or if vernier cor- 
rections are made in the first few hours of flight. For 
vernier correction, tracking data obtained in the first 
hour or so of the free-flight may be used to compute the 
actual burnout conditions accurately. The required cor- 
rection can then be determined and sent to the probe by 
radio command. 

Study of the angular orientation of the miss coefficient 
vectors has shown that for certain take-off days and 
burnout conditions, these vectors are almost collinear. 
The angular errors then appear equivalent to velocity 
errors and may be corrected with vernier velocity cor- 
rections on the probe. The closeness with which col- 
linearity can be achieved will vary sharply with launch 
date. Consider that the 8,, vector is resolved into com- 
ponents 6); and #1, parallel and perpendicular to the 
vector V,,, and let 8), be eliminated by vernier velocity 
correction. Fig. 20 shows the variation of the residual 
miss A. and #1 after early vernier correction (per 3° Ao 
and 8» error) asa function of launch date for 1959 Venus 
flights. The points shown correspond to the most col- 
linear operating region on each day. 

In addition to an early vernier correction, midcourse 
guidance corrections using long-term tracking data and 
commanded by radio can also be easily done. Suppose 
that the probe axis is stabilized in its flight on the sun- 
frame orbit at the attitude the axis had near burnout, 
which would be the case if the vehicle is spin stabilized. 
Then the vector Ab, representing the change in vector 
impact parameter from the nominal caused by a mid- 
course velocity impulse or “kick” will rotate in the miss 
vector plane as a function of the time of kick (Fig. 21). 
We may in effect resolve effects of burnout errors into 
components parallel and perpendicular to V,,, remove 
parallel components with an early vernier and perpen- 
dicular components with midcourse vernier rocket firing 
at the proper time. The midcourse corrections may also 
remove residual parallel components, taking advantage 
of the much more precise knowledge of the orbit gained 
by tracking in the early months of the flight. Such cor- 
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Fig. 21—Change in magnitude and direction of vector impact param- 
eter resulting from midcourse correction as a function of time of 
correction (Venus flight, June 8, 1959, V> =37,050 feet/sec). y* axis 
approximately parallel to the ecliptic plane; z* axis is perpendic- 
ular to the asymptotic velocity vector. 


rection schemes have the advantage of requiring mini- 
mal equipment in the payload since computation and 
decision functions will be carried out on earth using 
computers and communications equipment already re- 
quired for tracking and telemetry. 


CONCLUSIONS 


This brief review of interplanetary ballistic trajec- 
tories shows their extreme complexity, but it also in- 
dicates that a satisfactory method of analysis has been 
developed. The three-dimensional analysis of the prob- 
lems associated with these trajectories shows that con- 
ventional two-dimensional analysis is inadequate to 
understand the basic problems. If Figs. 4 and 5 are ex- 
trapolated back to burnout conditions in the earth’s 
field and then compared with Figs. 18 and 19, it is ap- 
parent that the two-dimensional miss coefficients are 
inaccurate by a factor of 5 to 10. Moreover, of course, 
two-dimensional analysis is wholly misleading in that it 
cannot consider the many practical problems, such as 
powered flight considerations, range safety, target ap- 
proach velocity, etc., which are essential to a successful 
mission. 

What this paper attempts to show is that if a proper 
trajectory analysis is carried out and ingeniously ex- 
ploited, existing propulsion systems, unsophisticated 
stabilization techniques such as spin, and relatively 
crude guidance systems are satisfactory for sending ve- 
hicles on interplanetary missions near each minimum 
energy time. However, if target approach velocity is not 
a central consideration, as in an impact case, or if com- 
munication transmission distance is not considered, then 
six to eight months each synodic period are available 
for launch to Venus and similar periods for launches to 
Mars. But such trajectories require an even more 
sophisticated combination of trajectory analysis and 
guidance techniques at burnout, post-burnout, and mid- 
course. The point is that whatever the accuracy guid- 
ance system adopted, or whatever the capability of the 
propulsion system, interplanetary flight can only be suc- 
cessfully achieved if the realistic dynamics of the trajec- 
tory problem are fully examined and understood. 
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Detection and Evaluation of Space Radiations” 
pEPAUL T. CORKHILL} anp R. A. HOFFMAN; 


Summary—In the very near future man will escape from his 
terrestrial environment and venture into space. The first flights 
will be limited space voyages circling the earth only one or sev- 
eral times at altitudes from 100 to 500 miles. As technological 
advances permit, these space journeys will become longer in 
duration and farther from earth until travel between planets in 
our solar system is accomplished. Ultimately of course, travel 
outside our solar system is envisioned. 

The hazards of space flight are numerous. An environment 
paralleling that of the earth’s must be simulated as closely as 
possible inside the space capsule. In providing such an environ- 
ment, those problems of air regeneration, temperature control, 
pressure control, etc., are overcome. One potential hazard from 
which man in space cannot be effectively protected is radiation. 
Therefore, the intensity of space radiations and their biological 
effects must be evaluated prior to prolonged space journeys if 
such journeys are to be accomplished without undue risk. 

The aeromedical field laboratory is now studying a space 
radiation detection system which would be capable of detecting 
and identifying the several forms of space radiation and the 
intensity of each form. Once the types and intensities of radia- 
tion prevalent in space are ascertained, it will be possible to 
study, in earth laboratories, the biological effects of most of 
these radiations. Radiobiological studies should lead to realistic 
RBE (relative biological effectiveness) values for these radia- 
tions, and thus biological dose rates (rem/hr) can be established. 

A radiation survey system such as described will be a very 
intricate electronic assembly. The biological studies to determine 
RBE values will require careful study. The unit must be com- 
pact, rugged, and lightweight if it is to be adaptable to satellite 
flights. 

Presently, a feasibility study is underway to determine the 
most practical and efficient electronic assembly. Then follows 
the fabrication of the complete detection system. Orbital flights 
of this system must then be made to evaluate the types and in- 
tensities of space radiations. Biological studies of these radiations 
must then be made so that biological dose rates may be estab- 
lished. Once the biological studies have been completed, the RBE 
factors can be incorporated in the radiation detection system, 
i.e., a system which will be capable of giving the space passenger 
a biological dose rate in rem/hr on a dial type instrument as well 
as his total accumulated radiation dosage during the flight. Ad- 
mittedly, such a system may not be a prerequisite for the first 
toddling steps of man into space, but it must certainly form an 
integral part of future space vehicles when bolder journeys farther 
from earth are made. 


AN is presently at the threshold of escaping his 
M terrestrial confines and venturing into space. 
Many problems have evolved and have been 

solved in approaching this new era of space flight; many 
are yet to be solved and undoubtedly others will be pre- 
sented after man has made his first orbital flights. Yet we 


remain undaunted, and certainly interplanetary flights 


* Manuscript received by the PGMIL, July 21, 1959. 
+ Operations Sec., Aeromedical Field Lab. AF Missile Dev. 
Center, Holloman AFB, N.Mex. 


within our solar system will become a reality in the not too 
distant future with space voyages to other solar systems 
following as technological advances permit. 

Technical advances in rocket design, propellants, elec- 
tronics, etc., have made possible the launching of artificial 
earth satellites, moon rockets, and tiny artificial planets 
circling our sun. The hostile environment of space presents 
additional problems to biological systems. These problems 
are now under study and must be solved prior to man’s 
initial venture into space. On earth, man is surrounded by 
an atmosphere to which he is adapted. Proper concentra- 
tions of oxygen, nitrogen, and carbon dioxide are present, 
and atmospheric pressure is such that gaseous exchanges 
can take place within our biological system. Also, our 
bodies are adapted to this pressure and any marked devia- 
tion will present serious consequences. The thin blanket of 
atmosphere surrounding the earth provides us with protec- 
tion against harmful components of solar radiation, cosmic 
radiation, and meteorites, besides serving as a buffer 
against sudden severe temperature fluctuations. Earth’s 
gravity provides us with a sense of weight which allows us 
to orient ourselves with our surroundings. These comprise 
only a few of the comforts which we enjoy on earth, and 
which are lacking in space. Design of a space vehicle in 
which the earth’s environment can be closely simulated is 
the ultimate goal. This is not impossible, but for at least 
two factors—one of providing weight, and second, of pro- 
viding protection against space radiations—it is extremely 
difficult. This paper is concerned with space radiations and 
will be limited primarily to this subject. 

The types of radiation prevalent in space are not 
thought to differ appreciably from terrestrial forms ex- 
cept for one type—heavy primary cosmic particles 
(Z > 6). These interact with our atmosphere and are 
broken into smaller secondary fragments before reaching 
earth. The primary cosmic radiation is composed of nuclei 
of elements ranging from hydrogen through iron in ap- 
proximately the following proportions—80-85 per cent 
protons, 14-17 per cent alpha particles and approximately 
one per cent heavier elements. When these primary par- 
ticles interact with the space atoms and molecules or are 
influenced by strong magnetic fields, they produce the 
other types of radiations found in space such as gamma 
rays, beta particles, positrons, 7 and w mesons, secondary 
protons, neutrons, neutrinos and secondary alpha particles. 
It is obvious, then, that the types of space radiations are 
similar to those found on earth except for these heavier 
primary cosmic particles. The energies of space radiation 
is quite another matter. While the energy of terrestrial 
produced radiations is limited in most instances to MEV 
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or a few BEV, the space radiations may possess energies 
of many BEV. The maximum energy levels are as yet un- 
known. It is presently thought that these high energy levels 
will appreciably affect the biological damage caused by 
these radiations. This will be discussed in more detail in 
another section of this paper. 

Detection of space radiation and determination of the 
intensities of the various components present still another 
problem. Certainly detection instruments are now available 
which are capable of detecting most of these radiations. 
However, it is desirable to have one detection system, light- 
weight, rugged, and compact, which, when exposed to a 
field of mixed radiation, such as in space, is capable of de- 
tecting the various types, separating each type into sub- 
groups according to energy, and finally determining the 
intensity of each subgroup. Obviously, this differs greatly 
from laboratory conditions where the type of radiation and 
even the energy distribution may usually be known. A 
study of a detection system with the capabilities listed, 
which would be adaptable to satellite use, is now in prog- 
ress. 

The biological damage caused by ionizing radiation is 
generally attributed to the formation of ion pairs as the 
ionizing particle or photon traverses matter. One possible 
exception to this mode of action would be the damage pro- 
duced by the heavy primary cosmic particles. These may 
produce minute damaged areas per particle as these par- 
ticles terminate in tissue. More study must be done in this 
area of radiobiology to determine if such an effect exists 
rather than the general toxic effect of free radicals as is 
the case with other radiations. Further, the damaging char- 
acteristics of one radiation relative to another (relative 
biological effectiveness, RBE) is thought to be related to 
the ionization density or linear energy transfer (LET) 
as the particle passes through matter. The greater the LET, 
the greater the biological damage which is due to a greater 
concentration of ion pairs per unit volume of tissue. 

The LET of ionizing radiation is influenced by the energy 
of the radiation and even more by the charge of the par- 
ticle. For example, an electron of 500 ky energy has a 
greater LET than one of 5 MEV and an alpha particle 
(Z = 2) has a much greater LET than an electron. As was 
mentioned earlier, the types of space radiation may not 
differ appreciably from terrestrial forms except for the 
heavy primary cosmic particles, but it is expected that 
much more energetic radiation will be encountered. Since 
the energy of a radiation is related to LET and this, in 
turn, to biological damage, it is obvious that the damage 
produced by these very energetic space radiations may 
differ considerably from the same type of radiation, pos- 
sessing less energy, which is found on earth. It is therefore 
necessary to ascertain the energy and intensity of the vari- 
ous space radiations and increase our studies of the bio- 
logical effects of the most energetic radiations available in 
laboratories so that realistic RBE values may be formed 
for the very high energy space radiations. These RBE 
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values would then be incorporated into the detection sys- 
tem allowing a biological dose rate in rem to be made. 
This is, of course, the ultimate goal—to have a radiation 
survey system which integrates the dosage attributable to 
each form of radiation into an over-all dose rate and total 
accumulated dose which could then be read on a dial similar 
to the speedometer of an automobile. 

A study of a space radiation detection system is now 
underway by the Lockheed Missiles and Space Division 
under a contract with the Aeromedical Field Laboratory of 
the Air Force Missile Development Center, Holloman Air 
Force Base, N.M. Two approaches to a detection system 
were considered. One approach would provide a dose 
rate and accumulated dosage but would not identify the 
various radiations. The dose rate on such a system would 
automatically be based on the LET of the radiations 
as measured by the system. The second approach would be 
a system of detectors capable of analyzing the radiation 
spectrum into its various components and energy levels. 
Subsequent incorporation of RBE values would allow this 
type of system to provide a biological dose rate and total 
accumulated dosage. The first approach would probably 
prove simpler in design, and, after evaluation of ambient 
space radiation with respect to types and energy levels, 
this may prove the system of choice for future routine in- 
corporation into space craft. It is felt that first the space 
radiations should be identified and the intensity and energy 
levels determined, and it was because of these factors that 
the second approach was chosen for the present study. 

The system chosen to survey the forms and energy of 
ambient space radiations will undoubtedly prove to be an 
extremely complex electronic apparatus. As was stated 
before, this unit should be able to survey all the space 
radiations encountered, subgroup some of the forms ac- 
cording to energy, determine the intensity of each group, 
and telemeter and record these data. Such a system must 
ingeniously combine almost all of the presently known 
detection instruments plus some possible new ones. The 
equipment must be miniaturized to the maximum extent 
still maintaining excellent performance. It must be light- 
weight and yet sturdy enough to withstand the rigors re- 
quired of space flight. 

We have limited our discussion to the detection and 
evaluation of space radiation and the possible biological 
hazards posed by these radiations. This is by far the pri- 
mary objective of the proposed radiation survey system ; 
however, any radiation in the space cabin produced by 
nuclear propulsion methods may also be monitored. The 
proposed unit will be capable of including such radiation 
in its survey since these radiations will not be unique 
among those encountered in space. Therefore, all ionizing 
radiations present in the space cabin, and the dose rate and 
dosage from these radiations will be adequately monitored. 

In summary, the primary objective of the radiation de- 
tection system now under study is to detect and analyze 
ambient space radiations. Information concerning types of 
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radiation encountered and energy levels obtained from 
initial flights should be used as the basis for additional 
laboratory studies on very energetic forms of radiation. 
Establishment of realistic RBE values from these studies 
and incorporation of these values into the detection system 
should provide a manner in which the dose rate and total 
accumulated dosage will be available to the space passenger 
at all times. 
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A radiation survey or monitoring system such as the one 
proposed will probably not be necessary for the initial 
manned space voyages at altitudes of 100 to 500 miles. As 
technological advances permit, man will journey farther 
into space—to other planets of our solar system and, 
eventually, beyond. Certainly the radiation hazard of such 
journeys must be evaluated, and constant monitoring dur- 
ing the flight accomplished, if undue risk is to be avoided. 


High-Accuracy Electronic Tracking of Space Vehicles* 
P. F. von HANDEL} anp F. HOEHNDORFF 


Summary—It is the purpose of this paper to show that optical 
accuracies can be reached with electronic tracking if the re- 
fractive index is known at the site of the radar, as well as at 
the site of the target. This is the case in tracking space vehicles 
with microwaves. The index can be measured at the tracking 
site and it is unity beyond the denser atmosphere. Ionospheric in- 
fluences can be neglected in the higher kilomegacycle range. It 
is also shown that the uncertainty of the course of any particular 
profile between the surface of the earth and the target has only 
minute effects on range and angular errors. Accuracy limits are 
determined and presented numerically. The inherent high pre- 
cision of modern electronic tracking systems can be fully utilized 
under these conditions; this is not the case in tracking vehicles 
flying in the atmosphere. 


HE inherent angular precision of modern electronic 

guidance and tracking systems is basically given by 

the ratio of the baseline between system antennas to 
the wavelength. Using frequencies in the higher kilomega- 
cycle range and long baselines, almost any desired precision 
can be reached. 

The inherent range precision is given by the measure- 
ment of an electric phase at a certain frequency. The phase 
angle at this frequency can be expressed by a correspond- 
ing time interval, t, which is the time an RF wave needs 
for a round trip from the illuminating radar to the target 
and back. The range to be measured is obtained by multi- 


@ 
plying = by the mean wave propagation velocity. The 


phase angle in conventional pulse radars, for instance, is 
measured at the pulse-repetition frequency (PRF). It also 
can be measured at an IF or at the RF. Assuming a given 
accuracy in phase measurement, the precision in measur- 
ing ranges is proportional to the frequency used. Going up 
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from measuring phases in the PRF to an IF, or to the RF 
again, almost any desired inherent precision can be reached. 
With modern equipment an angular precision of a few 
seconds of arc and a range precision of a few feet can be 
achieved. 

But the accuracy to be reached in actual tracking targets 
flying in the atmosphere is limited by wave propagation 
effects, that is, by refraction in the denser atmosphere. 
Errors in angle and range caused by the uncertainty of the 
refractive index profile exceeded the inherent precision of 
a system by up to two orders of magnitude. The system’s 
precision, therefore, could never be fully used. In fact, it is 
unreasonable to build a system to the highest possible in- 
herent precision, at considerable expenditure, when the ac- 
curacy of the final results is much inferior. 

The purpose of this paper is to show that in tracking 
targets beyond the atmosphere this situation no longer 
exists. The refractive index profile, of course, is still un- 
certain. The index can be measured at the sites of the 
electronic system on the surface of the earth, as usual. But 
now that the index at the target site is also known, it can 
be shown that it approaches unity to less than 10-* at 
altitudes beyond 200,000 feet. This is not the case for tar- 
gets flying in the denser atmosphere. In that case, the re- 
fractive index at the target site normally has a high degree 
of uncertainty and this is the reason for the uncertainty in 
angular and range errors. : 

The refractive effects of the ionosphere will not be 
treated in this paper because it can be shown that these 
effects are negligible if frequencies in the kilomegacycle 
range are used. We shall confine our considerations to 
these higher frequencies because modern high-accuracy 
tracking systems will, for other reasons, use these fre- 
quencies exclusively, 

The only other assumption we make is that of a spheri- 
cally stratified atmosphere concentric with the earth, ne- 
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glecting horizontal gradients of the refractive index. This 
assumption must be made, because such gradients along the 
path of a wave cannot be determined. The azimuthal error 
limits caused by such gradients will have to be determined 
by experiment. They are the same as those common to 
optical measurements in astronomy and are not likely to 
exceed a few seconds of arc. 

The new situation, then, consists of the fact that we now 
know the refractive index on both ends of the wavepath, at 
its origin on the earth’s surface and at its end on the site 
of the target. It will be shown that we now can define a 
family of standard refractive-index profiles, each member 
of the family starting at a certain index, mo, on the earth’s 
surface and approaching unity to less than 10-* at 200,000 
feet altitude. The errors in angle and range produced by 
these standard profiles will be computed. Proper correction 
terms, therefore, can be applied to the corresponding num- 
bers as measured by the tracking radar; consequently we 
shall call these terms correctable errors or corrections. 
They do not impair the accuracy of the system; the correc- 
tions are based upon measurable data and precise calcula- 
tions. 

The crucial question now arises concerning the magni- 
tude of the noncorrectable errors which originate in the un- 
certainty of the individual profiles between the known data 
at the beginning and at the end of the line. In Fig. 1 we see 
the site of the radar and three positions of a satellite. Of 
course, the scale of this picture is strongly distorted in 
order to make clear what we wish to show. It is easy to see 
that the refractive-index profiles cannot be measured from 
the radar at point O to all the points Si, S2,... 5; along 
a satellite path. The dashed lines indicate three of these 
propagation ways. 

Meteorological measurements are made along a more or 
less vertical line and each measurement needs a finite time. 
But the index profile is subject to changes in time and 
space and therefore it is impossible to cover the entire 
profile-plane at the proper instants. 

We then have to accept the uncertainty of individual 
profiles and herewith the neglecting of horizontal index 
gradients, and we have to find out the resulting uncorrect- 
able errors. 

In Fig. 2(a), two refractive index profiles have been 
drawn schematically, representing two extreme weather 
conditions. The index on the surface of the earth is the 
same in both cases and the index always approaches unity 
at 200,000 feet altitude to less than 10-°. At profile 1, the 
air is almost saturated all the way up to higher altitudes, 
whereas at profile 2 the air is dry from high altitudes down 
to the lower atmosphere and the humidity increases rapidly 
from there to the surface. A similar profile will occur with 
a temperature inversion in the lower atmosphere. In Fig. 
2(b) we see the two schematic wave paths corresponding 
to the two profiles. Path 1 has a slight curvature correspond- 
ing to an over-all low gradient of index vs altitude (profile 
1), whereas path 2 has a greater curvature in the lower part 
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Fig. 2—(a) Refractive index profiles; (b) wave paths correspond- 
ing to profiles. 


of the atmosphere corresponding to a high index gradient 
and a very low curvature in the upper atmosphere where the 
gradient of profile 2 is low. But note that both paths leave 
the atmosphere beyond 200,000 feet as parallel straight lines. 
This is the vital point in our investigation. No matter how 
the distribution of a stratified atmosphere might be, in other 
words, no matter how an individual profile might be, the 
wavepaths will always come out parallel beyond the atmos- 
phere. We designate a the zenith-angle measured by the 
radar. Then Aa, and Ag, at 200,000 feet are the correctable 
errors which can be precisely computed from information 
received only on the refractive index measured on the sur- 
face (mo). Both wavepaths are parallel (Ae, = Aa), which 
means that the uncorrectable errors are zero for a target 
at infinity. A real target will not be at infinity but relatively 
far from an altitude of 200,000 feet. The computation, 
therefore, will show that the uncorrectable errors decrease 
with increasing range. It will be seen later that they do not 
exceed a few seconds of arc in angle and a few feet in 
range for all practical purposes. 

The situation is quite different for targets inside the 
denser atmosphere (Fig. 2). Assume a target, say, at 
50,000 feet altitude. We see that the error angle Aa, is quite 
different from Ag, for the two profile cases. But the point 
is that neither Aa, nor Ae, can be computed because the 


164 


profiles are not known. In the case of a target inside the 
atmosphere, therefore, Aa, and Aa, are not corrections, but 
uncorrectable errors of considerable magnitude. Modern 
electronic tracking and guidance systems reach an inherent 
system precision which goes far beyond these uncorrectable 
errors. The system’s capacity, therefore, could not be fully 
utilized with targets in the atmosphere. But it can be used 
on targets beyond the denser atmosphere and there is no 
reason why it should not reach the known accuracy of opti- 
cal instruments without being biased by their restrictions 
on clear weather conditions. 


REFRACTIVE INDEX PROFILES 


The index of refraction, n, for microwaves penetrating 
the atmosphere, is a function of temperature, air pressure, 
and partial pressure of water vapor, according to: 


ay, 1) x10 = 2 (p+ =) (1) 
v= (n" = T 2, re 4 
where 

c=79 [°K/mb] 

T=temp [°K] 

p=air pressure [mb ] 

b=4800 [°K] 


e=partial pressure of water vapor [mb | 
n=refractive index. 


Figs. 3-5 show profiles of refractive index vs altitude 
(MSL) for medium, very low, and very high indexes at 
the surface. The very low case (Fig. 4) can only occur at 
the site of a high mountain. A meteorological sounding 
from Holloman Air Force Base has been used and the meas- 
ured data were transferred to the Sacramento Peak (10,148 
feet MSL), about ten miles east of Holloman AFB, N.M. 
The medium case (Fig. 3) was taken at Valparaiso, Fla. 

A sounding for the very high case (Fig. 5) was not 
available. Therefore, its construction has been based upon 
meteorologically reasonable assumptions of the temperature 
gradient and the supposition of a very high humidity 
(75 -- 90 per cent) in all altitudes. These conditions may 
occur some time before the outburst of a tropical thunder- 
storm. 

From measured data of pressure, temperature, water va- 
por content, and altitude, the profiles were calculated ac- 
cording to (1) and plotted in curves a, Figs. 3-5. In 
all three cases the relative humidity at the surface was 75 
per cent, as measured in Valparaiso. In order to obtain pro- 
files b and c, corresponding to somewhat extreme weather 
conditions at the same surface index, Vo = (mo-1) X 10%, 
we assume that at curves b the air is saturated at altitudes 
beyond 1000 feet above the surface while pressure and tem- 
perature distribution vs altitude remains the same. In 
curves ¢, to the contrary, the air is assumed to be extremely 
dry from high altitudes down to about 1000 feet and to 
become rapidly humid in a small transition layer from 
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Fig. 3—Refractive index profiles vs altitude: the medium case. 


there down, reaching 75 per cent relative humidity at the 
surface. The same conditions will also occur in case of a 
temperature inversion between zero and 1000 feet altitude, 
as shown by (1). 

Eq. (1) also shows that at altitudes higher than 50,000 
feet, for example, all the profiles m (or v) vs h will almost 
coincide independently of time and location because p is 
small, e is practically zero, and T does not vary noticeably. 

A y profile, as shown in Fig. 5 curve c, with such a dras- 
tic reduction of humidity in a small transition layer of 
only 1000 feet near the surface, seems to be meteorologi- 
cally unrealistic. However, it has been chosen intentionally 
in order to learn the uncorrectable errors in angle and 
range caused by two radically different profiles originating 
at the same surface index. 

Curves a and b in Figs. 3-5 were approximated by ex- 
ponential functions. Curves c were approximated by second- 
order parabolas in the lower part and by exponential func- 
tions in the upper part. The accuracy of these approxi- 
mations can be seen in the figures. 
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Fig. 4—Refractive index profiles vs altitude: the low case. 


At this point, we have to anticipate the results of the 
computations of angle and range errors made for the eight 
profiles shown in Figs. 3-5. The theory of these computa- 
tions will be treated in the next section. Table I shows the 
values in question (see Fig. 6). 


Qq=apparent zenith-angle (81.75°) 
s=geometrical range on the curved path 
vo = (mo — 1) 108 
Aa=angle correction 
d=real range to the target (4.864 X 108 feet) 
Sq = apparent range as measured by the radar on the 
curved path with apparent velocity v,=299,773.1 
[km sec? | 
(propagation velocity in vacuum c=299,792.9 
[km sec~?]) 
As, =d—S,=range correction 
A’a =Ac(a-curvey -AQ@. or ceurvey = uncorrectable angle 
error 
A?’s, =ASa(a-curve) —ASa(@. or c-curvey =uncorrectable range 
error. 
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Fig. 5—Refractive index profiles vs altitude: the high case. 


Table I shows the uncorrectable angle and range errors 
(A?a@ and A’sg) which originate from the fact that since 
only vp can be measured, it is not known which of the pro- 
files, a, b, or c, may actually prevail. They are two orders of 
magnitude smaller than the correctable profiles Ae and Aso. 

This fact clearly proves that there is no point in trying 
to fit meticulously the mathematical expressions, exponen- 
tial functions and parabolas, to the measured profiles. In 
other words, the fitting accuracy seen in Figs. 3-5 is 
of very minor importance. This becomes clear if we note 
that the uncorrectable errors are very small even though 
they result from such extremely different profiles as a, D, 
and c, and are calculated for a zenith angle of almost 82°, 
1.e., 8° above the horizon. 

It is not the uncertainty of the individual profile which 
is important but the value of vo. This means that the way in 
which intermediate layers are stratified does not substan- 
tially affect deviation in angle and range. The important 
factor is the surface index of refraction, vo, and this value 
must and can be measured accurately at the very instant of 
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TABLE I 
Curve a Curve } Curve ¢ Curves a, b, ¢ 
Aa AAS a5 4’ 42.6” 4’ 44.1” s-d 3.8 feet 
yo=250 Ma +0.9” —0.6” 
Fig. 4 Asa 282 feet 281.1 feet 282.8 feet 
A2s, +1.5 feet —0.2 feet 
Aw 8’ 7.4" ee Shes 8/10 .2” s-d 3.1 feet 
yp = 368 Ara +3.6” —2.8" 
Fig. 3 ASa 265.4 feet 259.8 feet 268.2 feet 
A’sq +5.6 feet —2.8 feet 
Aa 10’ 48.8” Same as curve a NO Sif s-d 3.1 feet 
vp =486 A2a —8.6" 
Fig. 5 ASa 384.8 feet 374.2 feet 
A’sq +10.6 feet 
p= poe hho (2) 
p= (vo: etho/#) ‘7 
Xo g WAVE PATH 
Ba v = vg a (3) 
Vo = tA A ee (4) 


Fig. 6. 


the tracking measurements and in the immediate vicinity 
of the radar. 

The vital point of our investigation lies in this fact. It 
entitled us to introduce “standard” profiles. 

Forty-four profiles which were taken at different alti- 
tudes (MSL) and seasons have been calculated according 
to (1). Among these were five individual soundings pro- 
vided by the U. S. Weather Service. Ten profiles were 
taken from “Climatic Data for White Sands Missile 
Area,” October, 1958. Twenty-two were taken from 
“Upper-Air Climatology of the United States,” U. S. De- 
partment of Commerce, Weather Bureau, T.P. No. 32 
(1957). 

All profiles taken on the East Coast and in the Continen- 
tal U.S. and Alaska could easily be approximated by an 
exponential function. All profiles taken from the coast of 
the Gulf of Mexico and the Caribbean islands showed a 
flatter slope of the curve at low altitudes. Nevertheless, in 
all cases there is a function, 


vy = vo:e 4, that fits the v profiles. 


In plotting values of v, the constant, H, for all the 
44 profiles, it was noticed that all points from White Sands 
(altitude 4000-feet MSL) and all points from coastal loca- 
tions (altitude 0 to 20 feet) were lying near straight lines, 
whereas all other points appeared to lie at random. This led 
to the attempt to reduce all values vo to sea level according 
to the transformation: 


where vg, is the refractive index reduced to MSL, 


H=a constant [feet] 
h=any altitude in feet above mean sea level (MSL) 
hy =altitude (MSL) of the observation site. 


The result is shown in Fig. 7. The dots represent values de- 
rived from factual soundings ; the crosses are derived from 
profiles constructed for Figs. 3-5. A smooth curve was 
then drawn through the points. Five checks for the points 
farthest from the curve showed that the approximation with 
the H values given by the curve was better than’the origin- 
ally chosen H indicated by these five farthest points. 

Let us define the H values given by the curve (Fig. 7) as 
“standard.” 

It will be noted that almost all points lie between the 
limits Hstanaara + 1500 feet. Later, therefore, we shall de- 
fine the uncorrectable errors, A?a@ and A2S,, as those which 
result from a deviation of + 1500 from the standard Hs; 
value at the same measured vo. 

From (4) a family of standard curves has been drawn 
in Fig. 8. For each vo measured at a station of altitude hy 
(MSL) the proper H value can be found by interpolating 
between the curves of parameter H. vg, is found in the in- 
tersection of the proper H curve with the axis of abscissa. 


COMPUTATION OF ANGLE AND RANGE CORRECTION 


General Equations 


The deviation of microwave propagation from a straight 
line in passing through the atmosphere can be calculated by 
Snell’s law of refraction. In its original form it holds only 
for the special case of plane surfaces: 


ee (5) 
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Fig. 7. 


For our purposes it has to be generalized to apply to 
media with curved surfaces. This derivation has been made 
by Knothe,? and independently by Dixon,? who stated from 
Fermat’s law and arrived at the same results: 


sin a 


R+ ho no 


= — - 6 
SIN ao R+h (0) 


(See Fig. 6.) 


It can be seen immediately that the general law (6) re- 
duces to (5) if the radii of curvature R +h) and R+h 
approach infinity, their differences remaining finite, which 
means the plane case. 

Any electronic range measurement between points O 
and S (Fig. 6) basically consists of the measurement of a 
phase; it may be an RF phase, an IF phase, or a PRF 
phase. In conventional pulse-radars, for instance, the elec- 
tronic phase angle is the PRF, which occurs between the 
illuminating transmitter pulses and the pulses reflected 
from a target is measured. At a known frequency the 


1P. yon Handel and H. Knothe, “The Curvature of Microwave 
Propagation Through the Troposphere and the Ionospheric Layers,” 
Holloman Air Dey. Center, Tech. Rept. No. 3; January, 1956. 

2H. M. Dixon, “Study of Refraction Errors in Radar Propaga- 
tion,’ White Sands Signal Corps Agency, Tech. Rept. 18A; April, 


1956. 
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Fig. 8—Graph for the determination of Vs: and H for standardized 
curves if ho and Yo are given. 
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phase angle can be expressed by a corresponding time in- 
terval, t, which is the time for a round trip of a wavefront 
from the radar to the target and back. The one-way trip, 


7 
therefore, needs Pe The distance to be measured is ob- 


tained by multiplying m by the mean signal propagation 


velocity, v. The propagation velocity v is not known ex- 
actly ; itis a function of the refractive index along the path 
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WAVE PATH 
K x 
ex 
Nees 
Pas eae 
ds 0G | 
c dh 
ZF (Rha 


and is subject to our calculations. Each electronic system, 
therefore, has a certain apparent velocity, va, which is cali- 
brated to the instrument as a constant. The measured value 
of range, indicated by the radar system, is an apparent 
range, Sa. It is a fictitious range which a wavefront would 
transverse if propagated with constant velocity, va, along a 
path s (Fig. 6). 


Seer 7 
5 ee (7) 


With sa measured and vg known (vq = 299,773.1 km), we 


T 
have —. 
2, 


Next, we have to determine the range s along the curved 
path. We note the well-known expression for the propaga- 
tion velocity in a medium with refractive index, n: 


C C 
v — = = —— (8) 
ve 
It follows from Fig. 9 that 
dh 
ds = ? (9) 
cosa 
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and we obtain from (6): 


R h 2 2 1/2 
cosa = E — sin? ap ( zn ‘ (“) | A(tO) 
R+h nN 


Substituting (10) in (9), we obtain: 


ike n(R + h) 
ae ho [n?(R +h)? — m0°(R + ho)? sin? aug |? 


.dh, (11) 


s is the geometrical range on the curved path (see Fig. 6). 
The real range we are seeking is d, the path on a straight 
line. We shall calculate it later. 
Meanwhile, we have to find a way to compute the alti- 
tude of the satellite, 4,, which is still unknown in (11). 
According to (8), we have: 


C 
ds = 2-di = —di, (12) 
n 
and combining (12) and (11), we find: 
n(R + h) 1 
c[n?(R + h)? — m02(R + ho)? sin? ao}? 

and 

1 hs 2(R h 
—-= if ae : - dh. (13) 
2 C ho [n?(R 43 h)? ao no(R + ho)? sin? arg |*/2 


From (13) we find hs by methods of numerical integra- 
tion on an electronic computer.® All the other values are 


t 
known: es is measured [see (7) ] and ~ can be found from 


(2) and (3): 


n= 1 (i — 1) ee: 


where H is determined by Fig. 8 for a measured m and a 
known fo. With hs now known, we find s from (11). 

Next, we have to determine the geometrical distance, d, 
(Fig. 6) that is the real range to the target. 


Applying the law of cosines of plane trigonometry, we 
find: 


d = [(R + ho)? + (R+ hy)? 


— 2(R + o)(R + hy) cos $s]? (14) 


In this equation @, is still unknown. In order to find ¢ as a 
function of any altitude, h, we proceed as follows. It can 
be seen from Fig. 9 that 

tan a-dh 


do = 15 
3 R+h ty) 


Eliminating « from (15) and (6), we obtain: 


3 All integrals have been calculated on a Remington Rand 1103A 
Univac Scientific. 
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6 oe le no(R + ho) sin ao 
} np (R +h) n2(R +h)? — m02%(R + ho)? sin? avo |? 


d can now be calculated from (14). 

Finally, we want to determine the correction in range, 
Asa, which has to be added to the apparent range, sa, meas- 
ured on the radar in order to obtain the real range d. From 
(7) and (14), we find: 


[(R + ho)? + (R + hy)? 


AS. =a — Ss, = 
= 2(R + Ii)(R + hn) cos by)! — = a. (17) 


The correction in angle, Ag, can readily be found in apply- 
ing the law of sines (see Fig. 6): 


8 


R 
sin (ao + Aa) = 


sin ¢s; 


and 


: Roar 
Aa = sin ( 3 sin c) — a. (18) 


Computations 


In all practical cases of tracking space vehicles, the fol- 
lowing values are known from measurements: 


Va=the propagation velocity calibrated with the ra- 
dar 

ho=elevation of observation site [feet, MSL] 

R=20, 891, 110 [feet] radius of the earth (sphere of 
equal volume) 

no=1-+7)X10-*=refractive index at the site of the 
radar antenna 

Oy) =apparent zenith angle as measured by radar 

Sa =apparent range as measured by radar 
H and vgz from Fig. 9; v from (3) 

¢ = propagation velocity in vacuum. 


According to (7), we obtain 


2, Va 


and solving (13), we find hs by methods of numerical 
integration. Next, we solve (16) obtaining g, and then 
(14) finding the true distance d. Finally, we determine 
the correction for range by (17) and the correction for 
angle by (18). 

In fact, a somewhat different course of computation has 
been used. We assumed the values of 


C 
no = (<); ho, se Q0, IF ip big 1? vgip el#, 


Va 


and s to be given. From these values, hs was calculated 
7 

from (11) and oo: =o from (13), gs from (16), d from 

(14) and, finally, Av and As, from (18) and (17). 
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(16) 


? 


This procedure enabled us to make some checks on the 
computational results, such as the respective values of d, s, 
and ¢, the magnitude of which must lie in this sequence 
whereas the difference s-d must be smaller than o-s. 
Furthermore, for any set of ¥o and a, 9-s is constant and 
s-d must increase with increasing s. The latter values must 
decrease with decreasing %. 

Eq. (11) has been approximated by step-by-step inte- 
gration according to Weddel’s rule with steps of width 
Ah = 2; 4; 8; 16; 32; and 64 feet, and 1620 steps in each 
case. This led to an altitude h, = 204,120 feet and a cor- 
responding s,. From here on » = 1 and the rest of the 
integral can be written explicitly ; hs can then be calculated 
from: 


h, = [(s — s,)? + 2(s — s)VA + (R+ hy)? [12 — R, 


where 
A =(R+h,)? — m0?(R + ho)? sin? a ¢ 


Eq. (13) was rewritten: 


Ts hs 
¢—=c= . 
Z h 


0 


This integral was calculated by the same method as (11) 
and in the same sequence and width of steps. Analogously, 
we define o, to be the value of ¢ at altitude 4,, and we ob- 
tain o=o,+ (s—s,) because if »=1, the integrals 
in (13) and (11) are identical. 

Eq. (16) was also integrated with the same procedure 
up to h, = 204,120 feet and the corresponding value 9, at 
that altitude. The explicit integration gives: 


$ ay $ che | cos no(R oe ho) sin ao 
's ir R a he 
no(R + ho) sin =] 
COS! at 7 ae 
Kicks 


Eqs. (14), (17), and (18) were calculated straightfor- 
wardly. No random errors could be found to an accuracy of 
0.001 foot and 0.001 second of arc. However, systematical 
errors, originating in the computer capacity, may occur. 
Therefore, the values of Asa and Aa listed in Table II 
should be considered accurate to two units in the last deci- 
mal. These values were calculated for the following param- 
eters: 

a = 87.5°; 85°; 80°; 70°; 50° and 10° 

s X 10-*[feet] = 5; 4.5; 4;3.5; 3.0; 2.5; 2.0; 1.5; 1.0 

vo = 400; 370; 340; 310; 280 
H{feet] = He, — 1500; Hot; Hee + 1500 
ho[feet | = 0. 
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TABLE II(a) 
2 ap=87.5° ao= 85° ay= 80° ao = 70° ao=50° ao= 10° 
vo Seay : ne Na nh Ara Aa Aa ; Aa Ara ? Aa La Aa zo 
wee (min) (sec) (sec) |(min) (sec) (sec) |(min) (sec) (sec) |(min) (sec) (sec) |(min) (sec) (sec) | (sec) (sec) 
400 5 23 OR OEE2ZO 8 TOMAS aS eeeiinh 7 2 Selene line 3 42.2 0.4 1 Sis Uni 14.45 0.01 
4.5 De SNA PHEAY) ey Lies wo) is 24,0. les 3 41.3 0.4 1 Sine) OR 14.45 0.01 
4 A See. UN. SRS = the 7 224 19 3 41.5 0.4 1 SMO Osa 14.44 0.01 
Sige) Die Day SRY) 15° 1S 2 1.5 7 20.4 159 3 40.7 0.4 1 Sibel Oe 14.42 0.01 
3 22 Se 2eZ SHO, 1S 244 Se A889 7 Whi oe eel 3 ANIA ORS) 1 30. OmOed 14.42 0.01 
Dos DN, 39) <6) UO) Tote ele OR BORG: 7 1450 2253 3 39 OLO 1 36.008 OFZ 14.39 0.01 
2 — — os 1! 6.83, IDES 7 Ses 27, 3 Siar WAG) 1 94 Hail 14.36 0.02 
5, — —_ — —_ — —_ 6 SS, geo 3 SATO me OKS 1 Sb UE 14°30" 70.02 
1 ; o— — 3 ak - AEG 1 Ser Was) 14.22 0.03 
370 5 200 3252220 We 2ifeee Ones 6 49.7 1:5 3 PS Go Mes} 1 SOeOs (eu 13.36 0.01 
4.5 DOR ASS eek 1 BS 10 6 A325 neo 3 24°58 0.3 1 2050 RO at 135360 10,02 
4 2OR SOUS oS WR DRS A 7 0d 6 AND ihe 3 DA 4 003 1 2928 a ORL nao. OO! 
3n0 20 ees Oe 38 WE PODS HS) 6 45.0) 138 3 23.9 OL4 1 29 irae Oak 16333) ORO 
S (Ors SH Le DAES IES OUR 7a) 6 4253, el. 9 3 23.2) 20.4 1 2935). 20.51 13:32) 0201 
Ds) O70) ee or6 12 LEBER S25 6 RG Peal 3 PRN OW 1 DOR Seen (eel 132308 OZ0L 
2 —_ _- — Limos 7a 5 9)25 6 SPR 9255) 3 AVIS OLS 1 DoS One 13226eeOn0n 
AS — — —_ — — _- 6 23545 553).0 3 is EY 1 28) 2 Oe 13,21) F001 
1 —- — — —_ — 3 WA Sl (0) il 26.9 0.2 yet ORO 
340 5 USB e422 eel 920 TD SPO sh SIGS) 6 | De 3 So leeORS 1 22 50m et 12271 0201 
4.5 Tea Aly} ibe PRATT 6 ilyacie ale: 3 Hoes WES} 1 22 50m On WAAK  OLOil 
4 1320) tel 9eS TSS Ore 6 LIF SSeS 3 7.4 0.4 1 PRS Oza 12226 30201 
Sis) 18 4.4 20.1 the “ieee Ged 6 Sy SIE) 3 6.8 0.4 1 224 Ot 12.24 0.01 
3 i Ae OES 11 RROD AREAL 6 Tel 1.8 3 6: be 2074 1 PPS West 12523 OLOt 
25 WH” ANB ae DAL ik) Sake oS 6 Qa Bw oS Mall =U).8) 1 21D Ort 1232S OxOn 
2 — — — 1O MSs) o.3 5 Se Sule 3 3.4 026 il PHNasy UE I 2507. 202 01 
15 —_ _- —_ — —- — 5 Jia ya PANT 3 OF Se ORG 1 Pers. War? 125127, 40302 
aL _ — — — 2 55105 2089 1 LOSS SO? 12.04 0.02 
310 5 Geno 224 eons IO Neh By 5 Sele ala! 2 Sle OOo 1 155 Site OF il 11513.) 20700 
4.5 Kel ASB MGA KORE BRIBE Shas! 5 Sad) 1k! 2 SOe7 One 1 1532S Ont ib Srivee OaCont 
4 16 Ae Gr lGeS 10 870s 255 5 SOR eso 2 Keine Es 1 15) Ot 1PSt7 OLOH 
30 1S OROmel On 10 1S so 5 34 5eale5 2 49.8 0.3 1 15.0 0EH 1.15 0201 
3 is) Sieg) Ie GA O) Bek Heil 5 He EO 2B 49.1 0.4 1 Wah (Obeak 1d ta OROT 
DES 15 Sey 7-9 OF 41.97 1635 5 Sr) ene Lie 2 48.0 0.4 1 14e5 oe Ont 112s 0 Ou 
2 —_— — — ONE A Aa gic 5 PIS Peet 2 46.4 0.5 1 LAAs Ost 112039 4001 
155 — _ — — _ — 5) 13 2bG cane: 2 AG ONE 1 13,45 SuOet 11) 025 0202 
1 a — — — _ Dy Soe E ORS 1 12> et OR 10.95 0.02 
280 5 14 26.0 12.8 9 6.5 4.4 5 3) Cel ieal 2 34.0— 0.3 1 Oe Or 10.09 0.01 
4.5 Aneel Roiael SAO 9 2.9 4.5 5 Dai Seal pp 39,0) 025 1 he Oe Opell 10.08 0.01 
4 14 6.4 13.2 Se BOSso) 456 5 ee As he 2 3320.3 Tt (paih <2 0ait 10.08 0.01 
he) LOE FOZ RZ S26 Sao 2930 43.9 4 92 ee lien 2 Ole Od 1 Oman Oe 10.06 0.01 
3 TS Meo lees Acul Saya one 4 SOnOi al eA: 2 SO Wee 1 ede ORL 10.05 0.01 
DS. 13 6.3 14.6 San ose Ole on 4 522 Oo 126 2 30.9 0.4 1 Telaeee ON 10:03 -0.01 
2 — — _— Carel 520 8) (6m 4 AA eel Be 293: .Oc4 1 One Onl 9.99 0.01 
iS — — — 4 SS Omens 2 26.7 0.6 if OOS 0) al 9.93 0.01 
1 oe ee — — — — — —_— —_— 2 21.4 0.8 1 AST “02 9.86 0.02 


In all the calculations, it was required that hs = 210,000 
feet in order to be sure that the target is outside the denser 
atmosphere. 

The differences of the values Aw and As, calculated for 
the three different values of H are denominated A?a and 
A’sq, respectively, and represent the uncorrectable errors 
caused by the uncertainty of an individual refractive index 
profile. These values are of considerably higher accuracy 
than Aa and Asg because the possible systematical errors 
are vanishing as well as the random ones. The values A?a 
and A’sq, listed in Table II, always represent the larger 
error in case the differences for Hs; + 1500 feet should 
not be equal. 


Accuracy Limits oF ELectronic TRACKING 


From Table II(a) it can be noted that considerable 
angular corrections, A%, must be applied specifically at 


large zenith angles, large distances, and large refractive 
indices. These corrections reach values up to more than 
twenty minutes of arc at an % = 87°. At zenith angles of 
% = 85° they decrease to about ten minutes of arc, and 
around % = 40° they diminish to less than one minute of 
arc. But even at @& = 10° from zenith they still amount to 
about 10 to 15 seconds of arc, which certainly is not 
negligible if the inherent accuracy of the system is better 
than two seconds. 

The uncorrectable errors, A?a, on the other hand, in- 
crease with decreasing range. The influence of & and Vo 
on A®a, however, is the same as on Aw. These uncorrectable 
errors are of somewhat alarming magnitude (up to % 
minute of arc at @ = 87.59, vo = 400, and s=2.5 X 
10° feet) for targets just appearing above the horizon, 
But A’a decreases rapidly from there and is down to about 
two or three seconds of arc at % = 80°, down to less 


=< 
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TABLE II(b) 

we sX10- ao=87.5° ap=85° ao= 80° ao=70° ag= 50° ay= 10° 
feet Asa feet A?sq feet~| Asa feet Asg feet | Asa feet As feet | Asg feet Asa feet | Asa feet As, feet | Asa feet A2sq feet 

400 5 168.6 8.7 25585 5.8 278.3 ore 302.0 17. 314.7 0.9 SD it 0.6 
4.5 135.8 8.6 202.5 See) 245.4 O80 269.2 1.8 281.8 0.9 288.1 0.6 
4 103.1 8.6 169.6 5.8 212.6 3.3 236.3 1.8 249.1 1.0 25 5—e 0.6 
3.5 70.2 8.7 136.8 5.8 179.8 3.3 203.5 eZ 216.2 1.0 DIED 0.7 
3 37.4 Sled 104.1 5.8 147.0 Be) 170.7 ihe 7 183.3 1.0 189.1 0.6 
pes 4.7 8.7 es 5.8 114.3 ne 138.1 1.8 150.4 0.9 156.1 0.6 
2 — — 38.8 5.8 81.7 ae) 105.3 1.8 15 iO oral 0.6 
RS —- — — — 49.2 Bis) Wt 1.8 84.7 1.0 90.2 0.6 
il — — — — —— -- 40.3 lead by) 1.0 SES 0.7 

370 5 LT es 235.0 5.3 278.3 3.0) 301.9 1.6 314.6 0.9 SpAle il 0.6 
4.5 138.9 7.8 203.0 Se: 245.3 3.0 269.1 1.6 281.7 0.9 288.1 0.6 
4 106.1 7.8 170.1 Sree 21255 3.0 236.2 1.6 249.0 0.9 250m: 0.6 
O25 13.3 7.8 137.3 So 179.7 3.0 203 .4 1.6 216.2 0.9 220 0.6 
3 40.4 7.8 104.6 S08} 147.0 3.0 170.6 1.6 183.3 0.9 189.1 0.6 
2a5 7.6 7.8 71.8 Ses 114.3 3.0 138.0 1.6 150.3 0.9 156.1 0.6 
2 -— - 39.3 ORS 81.6 Sell 105.3 126 117.4 0.9 123 0.6 
125 — — — —_ 49.1 3.0 UP Th 1.6 84.6 0.9 90.1 0.6 
1 — ~- — — — —_— 40.2 1.6 2 0.9 ies 0.6 

340 5 LR 6.9 236.9 4.8 278.4 2.8 301.9 135 314.6 0.8 Salil ORS 
a5 142.6 Us) 204.0 4.8 245.5 2.8 269.1 iS) 281.7 0.8 288.1 0.5 
4 109.8 120 Alea 4.8 DH? sh 2.8 236.2 15 249.0 0.8 2551 0.6 
Ono Hh al 7.0 138.3 4.8 179.9 2.8 203.5 15 216.2 0.8 D207 0.5 
3 44.2 6.9 105.4 4.8 147.0 2.8 170.6 fess 183.3 0.8 189.1 0.5 
25 11.4 6.9 72.8 4.8 114.4 2.8 138.0 bess 150.3 0.8 156.1 0.5 
2 — _- 40.3 4.8 81.9 2.8 105.3 15 14 0.8 1230 0.6 
15 _- — — — 49.3 2.8 72.7 125 84.6 0.8 90.1 0.5 
1 - — — — — a 40.2 185 2a 0.8 Silo O75 

310 5 178.8 Ont Rei as 4.3 278.5 BS 301.8 13 314.5 ORY 321.0 OFS 
4.5 145.9 6.2 204.6 Abs! 245.5 De 269.0 1.3 281.7 0.7 288 .0 0.5 
4 11322 6.2 Lids? 4.3 212.6 DS 23651 1.4 249.0 0.8 255.1 O55 
3:45 80.4 6: 1 138.9 4.3 179.8 Mes) 203 .4 1.4 216.1 0.7 222ml ORS 
3 47.6 6.1 106.0 4.3 147.1 BS 170.5 1.4 183.2 0.7 189.1 OFS 
225 14.8 6.1 73.4 4.3 114.3 IS 137.9 123 150.3 0.8 156.0 0.5 
2 — —- 40.9 4.3 81.8 BS 105.2 1133 117.4 0.8 123.0 0.5 
1.5 — —_ — — 49.3 2e5 72.6 1.4 84.5 0.8 90.1 OLS) 
il — — — — — — 40.1 lh 52.0 0.7 57.3 O85 

280 5 182.3 Bs) 238.2 3.8 278.3 2S) 301.7 Woo 31425 0.7 321.0 0.4 
4.5 149.5 BES 205.2 3.8 245.4 Dee, 268.9 122 281.6 0.7 288.0 0.5 
4 116.7 3.3) 17255 Sas 212.6 Dae 235.9 152 248.9 0.7 255.0 0.5 
58) 83.8 5.4 OES SES 179.8 BD 203 .2 12 216.1 0.7 222.0 0.4 
3 51.0 5.4 106.7 3.8 147.0 BS) 170.4 i 183.1 0.7 189.0 0.4 
D225 18.2 ae 74.0 3.8 4S DS 137.8 ge? 150.2 0.6 156.0 0.4 
2 — —- 41.4 3.8 81.7 DD 105.0 12 Lire) Ons 122.9 0.5 
15 — a —_ 49.2 DD) (ee! fe2 84.5 OF, 90.0 0.4 
1 _ — _ — — — 40.0 12 51.9 0.7 Sie 0.4 


than one second at % = 70°, and down to about one hun- 
dredth of a second at % = 10°. This means that we reach 
the inherent precision of high accurate radar systems al- 
ready at angles % < 80°. 

Table II(b) shows the range correction, Asa, and the 
uncorrectable range errors, A?sq. Let us first discuss the 
errors, A?sqa. They depend only on % and vo and are inde- 
pendent of the range, s. At large zenith angles % and large 
refractive indexes v) they reach a magnitude of almost ten 
feet; at % = 0° to 50° they amount to about one foot, and 
diminish to about one half foot at zenith. 

The range corrections, Asa, on the other hand, increase 
with decreasing %. For angles a < 80° the values are 
practically independent of vo. They increase with increasing 
range. At a first glance, it may be surprising that, contrary 
to the family of Aa curves (Fig. 10), the corrections Asa 
(Fig. 11) increase with decreasing % and with increasing 
distance, d. This is because the apparent propagation veloc- 


ity, Va, calibrated to the radar as a constant, was 299,733.1 
km in our case and all calculations in this paper have been 
based upon this figure. If another constant wv had been 
calibrated into a specific radar, the respective number 
would have to be changed in (7) and (17). If, for instance, 
light velocity in vacuum (c = 299,792.9 km) had been 
used for va, then the family of curves on the left side of 
Fig. 11 would show increasing errors Asg with increasing 
% and decreasing d. 

Here again, it can be seen that the range corrections, 
ASq, are in the order of a few hundred feet and must be 
allowed for if accurate data are expected. 

The exact results of our computations can only be seen 
from Table II(a) and (b). In order to enable a general 
survey on trends of Aa, A’a, Asa, and A’5sq, the content of 
Table II(a) and (b) has been plotted in Figs. 10 and 11, 
respectively. Exact values, of course, cannot be determined 
from these figures because they lie below the drawing 
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Fig. 11—Range corrections. 


accuracy. The ordinate in both figures is a. Precisely, the 
curves in the middle of Fig. 10 split up appreciably in their 
lower part (% > 80°) for values of vs; 4 340, as can be 
seen by comparison with Table II(a). These divergences, 
however, have been omitted to keep the figure readable. 
The same is true for the right-hand part of Fig. 10, where 
the curves of A’ also diverge with d in the lower part. The 
same is true for the right-hand part of Fig. 11, where 
A’s, depends on vgz. 


REMARKS ON OPTICAL TRACKING 


From the previous paragraphs we have learned that 
range must be introduced to add further angular correc- 
tions. The same thing, of course, is true for optical track- 
ing of satellities. In Fig. 12 the configuration is shown 
schematically. A satellite, S’, appears in coincidence with 
a star, *’, under a zenith-angle a. The true position of the 
satellite is indicated at S and the true direction to the star 
is shown at ®. It is easy to see that the correction angles 
Aw, and Aa,, respectively, are not equal. They diverge 
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Fig. 12. 


increasingly with decreasing range of the satellite. Aa, is the 
normal astronomical refractive correction. In order to obtain 
a better correction, we have to subtract (Aa, — Aa;) from 
Aa,. The remaining uncertainty caused by the uncertainty 
of individual profiles is characterized by A’a« in Table 
II (a). It is the same as for electronic tracking. 

Table III shows the necessary corrections Aa, — Ads 
dependent on @ and d. 


TABLE III 
OPTICAL CORRECTIONS 


a Ao — Aas in sec of arc 

Se 10° 50° 70° 80° 
d 
1.108 feet 0.33 4.0 etl — 
2.108 0.20 Qe2 9.3 39.9 
32108 0.14 Thsgs) 6.6 30.6 
4.108 0.11 iat 523 25.9 
SoilOy 0.10 1.0 4.6 DSi, 

CoNCLUSIONS 


We have found that angular accuracies of a few seconds 
of arc and range accuracies of less than ten feet can be 
reached in electronic tracking of satellites. This precision 
can be obtained by measuring only the refractive index at 
the radar site without any knowledge of the course of an 
individual index-profile along the way of propagation. The 
accuracy limits caused by wave propagation effects reach 
the inherent precision of high accurate radar systems al- 
ready at 20° above horizon (a < 80°) if the range to the 
target is larger than 3 million feet. Therefore, it will now 
be profitable to put a new emphasis on developing elec- 
tronic systems with high inherent precision. In both the 
electronic and the optical tracking procedures range must 
be known in order to introduce accurate angular correc- 
tions. In this case, the same accuracy limits are reached 
electronically and optically. 
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Some Reliability Aspects of Weapon and Other Space 
Systems Using Ballistic Missile Boosters* 
HENRY TRIWUSH}+ 


Summary—This paper shows that there is a need to differen- 
tiate, at least on occasion, between weapon and nonweapon 
space systems where their reliability is concerned. Three criteria 
are offered as a means of determining whether this differentiation 
need be made in the case of any particular nonweapon space 
system. Where this differentiation should be made, possible 
reliability approaches for the nonweapon space system are sug- 
gested. 


INTRODUCTION 


66 ELIABILITY is everybody’s business.” This 
R well-known platitude is often only the preamble 
to such witticisms as “And what’s everybody’s 
business turns out to be nobody’s business.” But let us con- 
sider it seriously. The fact is, where weapon and other 
space systems are concerned, reliability is everybody’s busi- 
ness because of these systems’ costs, complexities, and na- 
tional importance. “Everybody” includes both technical and 
nontechnical people, within the armed services and outside 
them. The topic has even attracted Congressional interest.* 
With all of this concern, it is essential that certain aspects of 
the reliability of such systems be brought more clearly into 
focus. This involves consideration of the differences be- 
tween the two kinds of systems, as well as of the reliability 
requirements for both. 

Both weapon systems and other (operable in peacetime) 
systems are defined in an initial section of Air Force Regu- 
lation 70-9.2 The sole criterion used to distinguish be- 
tween the two is that of intended use: a weapon system 
“forms an instrument of combat” (to quote the above regu- 
lation) ; a peacetime system—such as a nonweapon space 
system—does not. In other words the ultimate (i.e., 
planned) use of a weapon system involves its being “fired in 
anger.” It has been common usage (as is later found in the 
above regulation) to let the term “weapon system” be ap- 
plied to both kinds of systems, and for many purposes 
this is acceptable. However, when considering the reliability 
requirements of systems involving ballistic missile boosters, 
it will be shown that weapon systems cannot always be 
lumped with other space systems. 

A fundamental concept in considering reliability require- 
ments is that these requirements are useful only insofar as 


* Manuscript received by the PGMIL, July 21, 1959. 

+ Major, U. S. Air Force, Ballistic Missile Division, AF Unit 
P.O., Los Angeles, Calif. 

1See, for example, “Department of Defense Appropriations for 
1960, Pt. 1. Policy Statements,” in “Hearing Before the Subcom- 
mittee of the Committee on Appropriations, House of Representa- 
tives, 86th Congress, First Session. Subcommittee_on Department 
of Defense Appropriations,” U. S. Govt. Printing Office, Washing- 
ton, D. C., pp. 258-295; 1959. 

2“Procurement and Contracting. Airborne Weapon System and 
Support System Procurement,” Dept. of the Air Force, AF Regula- 


tion 70-9; 1956. 


they aid in planning the procurement and/or operation 
of the system to which they are applied. It follows, there- 
fore, that they should be applied only to the extent to 
which they are useful; i.e., there should be no leveling of 
numerical reliability requirements on a given system 
merely for the sake of doing so, or because other systems 
are so treated, or because of other similar reasons. 


BALuistic MisstLE WEAPON SYSTEM RELIABILITY 


The desired output of any weapon system is the delivery, 
when desired, of a certain known amount of damage-creat- 
ing energy upon a target. Regardless of the amount of the 
energy to be delivered (ranging from that of a small arms 
bullet to that of a megaton-yield warhead), the desired 
performance of the weapon system under given conditions 
of range or duration of operations, etc., may be tied to a 
single characteristic, the weapon accuracy. In other words, 
successful performance may be defined as designed energy 
delivery with an accuracy greater than some arbitrary 
limit (usually some simple function of the circular proba- 
ble error). The probability that the weapon system will 
attain this performance under the given conditions of 
warning or preparation time, range or duration of opera- 
tion, etc., is the numerical value of its reliability.® 

For planning purposes, a known minimum numerical 
value of the system reliability must be assigned (with a 
known degree of confidence) to a weapon system involving 
ballistic missiles. This known reliability value is necessary 
for planners either to 1) estimate the most probable amount 
of damage that a given ballistic missile striking force, 
launched at a randomly occurring (uncontrolled) point in 
time, can inflict on an enemy, or 2) determine how large a 
ballistic missile striking force must be procured to be able 
to inflict a desired amount of damage, with a given degree 
of confidence, on an enemy under the same time considera- 
tion asin 1) above. Underlying this consideration is the con- 
cept that this country’s national policy does not permit us 
the initiative in determining the exact time at which these 
weapon systems will be used in their ultimate operational 
environment, i.e., “fired in anger.” (With the present in- 
creasing emphasis on such ballistic missile weapons, maxt- 
mum permissible national reaction times are reduced to 
fractions of an hour, strengthening the validity of this con- 


cept. ) 


’ Note that for ballistic missiles the condition of warning or prepa- 
ration time is at least as important as that of range or duration of 
operation. In fact, ballistic-missile weapon-system reliability is fre- 
quently broken down into at least three factors: probability of being 
ready to launch whenever desired, probability of achieving successful 
launch whenever desired, and probability of achieving successful 
flight. In any event, numerical values of both minimum accuracy and 
reliability are required for planning purposes. 
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RELIABILITY OF A NONWEAPON SPACE SysTEM USING 
Batuistic MissiLE BoosTERS 


At present, budgetary considerations influence the plan- 
ning for such a system’s size (i.e., the number of ballistic 
missile boosters to be assigned to the system program) 
much more than for weapon systems. This is because of 
the fact that while peacetime-operable systems can con- 
ceivably be used for any of the system functions given in 
Air Force Regulation 70-9 (logistics, training, intelligence, 
etc.), the current state-of-the-art limits these programs to 
exploratory, semi-R and D efforts. The present cost and 
scarcity of ballistic missile boosters (including the neces- 
sary manpower, facilities, and supplies for their operation ) 
puts a practical limit to the number of boosters available 
for such a program. As a result of this budgetary rather 
than operational limit, a space-system program of the kind 
being discussed in this paragraph, which will usually have 
many objectives (i.e., desired outputs), will not have exactly 
the same objectives for any two attempted launches. In 
other words, there are no “space launches’; the failure 
(partial or complete) of any launch represents an ir- 
reparable failure to achieve a portion of the program’s 
over-all objective. This implies a reliability requirement of 
100 per cent which is unrealistic.‘ 

Although future systems of this nature may not suffer 
from a budgetary, “no spares’’ limitation, the multiple- 
objective possibilities will always have to be considered. 

Finally, the initiative for the use of such systems, since 
they operate during peacetime, at present generally rests 
with the planners. Therefore, no requirement imposed by 
the need to consider limited preparation or warning time, 
exists for system reliability. This is in direct contrast to the 
case for ballistic missile weapon systems. 


RECOMMENDATIONS AND CONCLUSIONS 


It can be seen that space systems which are operable 
during peacetime and which use ballistic missile boosters, 
merit careful consideration as being distinct from ballistic 
missile weapon systems when their reliability aspects are 
involved. This does not mean that there will always be 
differences, but that sufficient basic differences may exist in 
individual system programs to require a different approach 
to expressing and achieving system reliabilities. Each non- 
weapon space system must be considered individually with 
regard to the following three criteria: 


*A secondary effect implied by this multiple-objective considera- 
tion is that the hardware will vary sufficiently from launch to launch 
so that success or failure in earlier launches of a program can no 
longer be used to predict the reliability (in terms of probability of 
success) of the later launches of the program. This is, of course, in 
direct opposition to the practice for ballistic missile weapon systems. 
This can be put into other terms by saying that the multiple-objective 
nature of such a peacetime-operable program makes the expression 
of its reliability as a numerical probability at least ambiguous, if not 
completely without meaning. 
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1) Does the system have only one all-important objec- 
tive? 

2) Can this objective be expressed in terms of a single 
system output characteristic corresponding to the accuracy 
of a weapon system? 

3) Does the program for the system envisage sufficient 
material, when compared to the demands of the system’s 
over-all objective, to permit redundancy, ie., “spare 
launches” (less than 100 per cent reliability requirement) ? 

Only if all three of these questions can be answered in 
the affirmative can the nonweapon space system be con- 
sidered as if it were a weapon system so far as its reliability 
is concerned, i.e., its reliability can be expressed as a 
simple probability. 

For a system which does not meet the above criteria, the 
approach again is individualistic. Perhaps numerical ex- 
pressions of reliability, if required, may be given in terms 
of maximum permissible losses of system outputs which 
characterize the system’s objectives.** This implies that 
failure to attain these objectives to any extent whatsoever 
is permissible, but this may not be the case. If no such 
failure is permissible for planning purposes, then all that 
can be established as a system reliability requirement is the 
maximum possible reliability commensurate with available 
funds and the state of the art. Imposition of firmer de- 
mands (in the form of numerical values, however ex- 
pressed) merely hampers development of the system by 
misdirected application of the limited funds and manpower 
available for this development. 

Finally, the stupendous costs of the types of system dis- 
cussed here make it extremely unlikely that tests of the 
entire system will ever be used as the sole means to dem- 
onstrate whether any predetermined system reliability re- 
quirements have been met at the end of the system develop- 
ment period. This means that any such demonstrations must 
involve tests of relatively inexpensive components of the 
system. This, in turn, implies the placing of a further re- 
quirement on the space system reliability approach: that it 
be of such a nature that any numerical system reliability 
requirements can be resolved into relatively easily de- 
monstrable, component reliability requirements, preferably 
at the component part level. 

A system reliability approach using the above concepts 
will place no heavier loads than those presently placed on de- 
sign and reliability engineers working on component and 
component part levels; it will insure that these engineers’ 
efforts will contribute more effectively to the desired system 
reliability. 


*A. C. Block, “Methods of relating performance to reliability,” 
Proc. Fifth Natl. Sump. on Reliability and Quality Control in 
Electronics, Philadelphia, Pa., pp. 389-395; January 12-14, 1959. 

° J. L. Burnside, “A Proposal for the Specification of Reliability in 
Military Equipment Contract,” Lockheed Aircraft Corp., Missile 
and Space Div., Sunnyvale, Calif. Rept. No. LMSD-421249; 
February 15, 1959, 
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Radar Beacons for IRBM/ICBM* 


F. J. CLARK+ 


Summary—This paper summarizes problems which confront 
engineers who are responsible for the design of radar beacons 
and beacon antennas used in IRBM and ICBM test vehicles, The 
most important factors which affect compatibility between the 
beacon and tracking radars at the test range are enumerated. 
Cited also are several design areas which should be correlated 
with the test range prior to finalization of beacon design criteria. 


INTRODUCTION 


ESIGN engineers concerned with radar beacons 
and beacon antennas for IRBM or ICBM test 
programs are confronted with many difficult inter- 

related problems. The general scope of those problems 1s 
probably not exceeded in any other area of electronic de- 
sign activity. 

In brief, the items to be considered by the engineer 
include simultaneously the requirements associated with 
uncompromising physical and operational environments. 
Specifically, he will be concerned with environmental ex- 
tremes regarding temperature, pressure, vibration, and ac- 
celeration, His design criteria will be strongly influenced 
also by limitations on size, shape, weight, and power con- 
sumption. With respect to the test range environment, he 
will be fortunate if his initial activities or other commit- 
ments leave time to investigate thoroughly and to correlate 
the beacon design requirements with the capabilities of the 
tracking radars and their specific locations or planned 
utilization, The latter situation is one that will always re- 
quire special attention because of the continual efforts of 
the test range to increase the performance capability of its 
radar facilities. 


COMPATIBILITY CONSIDERATIONS RELATING TO 
CHAIN-RADAR OPERATIONS 


Successful radar beacon utilization in ballistic missiles is 
(usually) based on the necessity for simultaneous tracking 
by two or more widely separated radars (1.e., chain-radar 
operation). This requirement, which presumes a Satis- 
factory degree of compatibility between the beacon and 
the tracking radar, establishes basic criteria for the beacon 
design parameters listed below. Since the operational as- 
pects for each missile test program differ, it is not feasible 
to suggest that either of these parameters is more 1m- 
portant than another. Experience indicates, however, that 
deficiencies or incompatibilities in the following design 
areas cause most tracking failures in chain-radar opera- 
tions: 

1) Beacon antenna radiation pattern and/or polariza- 
tion—this is also allied with changes in aspect angle and 
in-flight stability of the missile. 


* ‘ot received by the PGMIL, July 21, 1959. : 
Se Isienciatea Div., Directorate of Range Dev., Patrick 


AFB, Fla. 


2) Poor frequency stability characteristics in transmit- 
ter and/or receiver. 

3) Insufficient radiated power. 

4) Inadequate receiver sensitivity. 

5) Unfavorable AGC characteristics. 

6) Too stringent a requirement relative to optimum 
spacing of (radar) interrogation code pulses. 

7) Long recovery time. 

8) Poor receiver selectivity. 


Design Considerations for Beacon Antenna 


The engineer should be aware of the type of antenna 
scan and polarization utilized by the tracking radars when 
formulating design criteria for the beacon antenna. For 
example, a simultaneous lobing, amplitude comparison, 
monopulse radar is especially sensitive to interferometer- 
type nulls in the radiation pattern when two or more bea- 
con antennas are driven simultaneously. As a typical case, 
four such antennas spaced 90° around a specific missile 
produced several 30-db nulls approximately 4.5° apart. 
As the radar looked at this portion of the pattern, discon- 
tinuities in the angle tracking data occurred each time 
missile roll corrections caused the radar to scan through a 
null. 

Another important consideration involves a practice at 
the test ranges of employing separate “skin” and “beacon” 
receivers to provide simultaneous radar display of skin and 
beacon signals. This is a range safety measure which per- 
mits transfer of the tracking function to the skin return 
in case of beacon failure at or soon after missile launch. 
This procedure, with current radars, sacrifices 3 db in the 
beacon reply due to signal “splitting” in the radar mixer. It 
also poses a question to test range personnel concerning the 
optimum radar antenna polarization for simultaneous skin 
and beacon tracking. Present opinion at AF MTC favors 
linear polarization for the radar and circular polarization 
for the beacon. After further investigation it may be deter- 
mined that this idea should be reversed, 12.e., circular 
polarization may be preferred on the radar and linear 
polarization on the beacon. Either arrangement, however, 
will introduce another 3-db loss (each way) in the radar- 
beacon transmission loop. 


Beacon Magnetron Frequency Stability Problem 


One magnetron design parameter which appears to have 
received only passive interest involves frequency change 
vs pulse repetition rate. When a second and perhaps third 
radar in a chain network phases in with beacon interroga- 
tions, minimum shift in beacon magnetron frequency is 
desired. Otherwise, the receiver AFC may fail to maintain 
proper tuning of the radar with the beacon transmitter. 
Since tracking failures cannot easily be post-related to this 
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deficiency, it deserves careful consideration in missile- 
beacon applications. 


The AGC Problems 


Simultaneous interrogations by more than one radar in- 
fers that the beacon must alternately respond to strong 
and weak signals. At times, the level of one signal may be 
near zero dbm and the other(s) anywhere from —30 dbm 
to —65 dbm. This requires special attention in receiver 
design, but anything less than instantaneous AGC in the 
beacon is undesirable for chain-radar operations. 


Considerations Relating to the Spacing of Interrogation 
Code Pulses 


For multiple pulse interrogation codes, close spacing 
and tight tolerance on the position of pulses impose a prob- 
lem in calibrating the radar pulse encoder at distant range 
stations. Multiple pulse codes are used principally to reduce 
interference from undesired radars. Close pulse spacing 
with narrow position deviation limits provides additional 
safety against such interference. However, for a radar to 
transmit reliably and consistently such multiple pulse codes, 
rigorous and precise calibration of the radar encoder is 
necessary. Previous chain-radar experience leads to a 
question in this case concerning the validity of the earlier 
philosophy of “close spaces—tight tolerances.’”’ Factors 
bearing on this problem which encourage some relaxation 
in requirements are: 
1) The reliability of precision calibration equipment, 
after exposure to a humid climate and shipment to 
a remote range site, is actually indeterminate. 

2) A reliable estimate on the long-term uniformity of 
calibration equipment at each of several widely sepa- 
rated sites cannot be established. 


Recovery Time Considerations 


Recovery time relates to the ability of a beacon to reply to 
each interrogation received. Relatively high interrogation 
rates, such as encountered with chain-radar operations, re- 
quire a reasonably short recovery time. Currently, with 
three “long-range” radars in the chain, a maximum recovery 
time of 300 usec is permissible. Future radar improvements 
to provide automatic sequencing, or interlacing of interro- 
gations from several radars, will reduce maximum permis- 
sible beacon recovery time to approximately 50 usec. 

Minimum recovery time is a problem to the extent that 
ring-around, or self-triggering from the beacon’s own 
transmissions should not occur. If recovery time is less 
than approximately 12 wsec, however, the beacon reply 
signals reflected from nearby structures (in the launch 
area) are likely to cause self-triggering. 


Beacon Transmitter Power and Receiver Sensitivity 
Considerations 


The rated power output of the beacon, of course, has a 
bearing on the predicted beacon signal level that will be 
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available at the radar receiver. To have any real signifi- 
cance, however, the beacon output power should be appro- 
priately related to the deepest null through which any par- 
ticular test range radar is expected to track. With this 
qualification, the beacon transmitter power is of most in- 
terest as one of the factors in the communications loop 
between the radar and the beacon. Three other items of 
interest in that loop are beacon receiver sensitivity, radar 
transmitter power, and radar receiver sensitivity. Systems- 
wise, the objective always is to have the quantitative values 
for those items so related that the ability to communicate 
from the radar to the beacon will be identical to that which 
will be provided between the beacon and the radar. When 
that condition is achieved, the two-way path between the 
radar and beacon is referred to as a balanced-loop. Due to 
various technological considerations, responsibility for 
meeting the balanced-loop criteria is normally placed on the 
beacon designer. He must therefore have up-to-date infor- 
mation on the characteristics of test range radars before he 
can establish optimum design limits for the beacon receiver 
or transmitter. Furthermore, he needs such information in 
order to make intelligent compromises in the beacon design 
which, for many special purposes, may be necessary. 


Receiver Selectivity Considerations 


Beacons with crystal-video-type receivers inherently 
have insufficient sensitivity and poor frequency selectivity 
characteristics. Because of their low sensitivity, they are 
generally considered unsatisfactory for chain-radar opera- 
tions. They are also considered unsatisfactory in this ap- 
plication because of their susceptibility to interference 
from radars operating in each of the primary radar fre- 
quency bands. Experience reveals that such a nonselective 
beacon will frequently be triggered by an adjacent channel 
CW signal, or by a radar operating at a frequency sepa- 
rated from the nominal frequency of the beacon receiver 
by more than 2500 me. 

Even with super-het receivers, the designer should keep 
in mind that poor selectivity will create a troublesome prob- 
lem. The high concentration of electronic radiating devices 
at a test range, operating at many different frequencies 
with all types of modulation and power levels from frac- 
tions of watts up to several megawatts, presents a more 
severe electronic environment than is normally experienced 
in any other beacon application. 


OTHER GENERAL COMPATIBILITY CONSIDERATIONS 
Data Accuracy vs Radar S/N Ratio 


In most missile radar beacon applications, the ability to 
estimate quickly the quality of the radar data is desirable. 
The simplified equations presented below provide a con- 
venient method for doing this, assuming that the beacon 
transmitted power and radiation pattern are known. These 
equations have been verified under controlled tests using a 
specific monopulse radar. However, similar equations may 
be derived for other precision tracking radars: 
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RMS angle noise error (mils): 
6 


O79 = 


(1) 


V2 fr syn 
B, 


where 


@ is antenna beamwidth (mils) 
S/N is IF signal-to-noise power ratio in the sum chan- 
nel (for S/N>6 db) 
f, is radar repetition rate, and 
B,, is equivalent noise bandwidth of angle servos. 


RMS range error (yards): 


oR = ———— (2) 
[fe S 
By IN 
where 


7 is pulse width (in yards) 
B,, is bandwidth of radar range servo, and S/N and 
f, are the same as in (a) above. 


These equations are of interest at AFMTC since they pro- 
vide convenient, slightly conservative estimates of radar 
performance capabilities over any specific portion of a mis- 
sile’s trajectory. The signal-to-noise ratio is the most 
significant factor, since it is the only variable which is not 
under direct control of the radar operator. Note that for 
tracking high-velocity close-in targets, the wider servo 
bandwidths that must be used require higher signal-to-noise 
ratios—nominally, 19-23 db—to maintain a specified level 
of system accuracy. 

Since the use of these equations presupposes a one-to- 
one beacon reply for each radar interrogation, a word of 
caution regarding them is added. If the beacon receiver 
sensitivity or beacon antenna pattern characteristics are 
deficient, the beacon will “countdown” (not reply to 
every radar interrogation). Beacon countdown directly 
affects the quality of radar data, but as a rule of thumb, 
countdown of up to ten per cent may be permissible when 
tracking missiles out beyond 150 miles. In order for the 
radar to maintain automatic track, the maximum allowable 
beacon countdown is established by the commutation 
(sampling) rate of the error signal from the radar track- 
ing servos. This commutation rate in a monopulse radar is 
usually 14 to % of the radar PREF. To maintain automatic 
track (even though the data will be very poor) the radar 
should receive one beacon reply for each half-cycle of the 
error commutation signal. 


Beacon Reply-Delay Considerations 


A choice of beacon reply-delay settings is desirable when 
performing prelaunch pad checks with the radar. This fea- 
ture will permit the separation of beacon replies from the 
radar signals which are reflected from missile test stands 
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or other structures in the vicinity of the launch pad. The 
provision of at least two fixed-delay settings is desirable. 
Furthermore, simultaneous tracking of more than one bea- 
con in the same missile will be possible if a different delay 
setting for each beacon is available. 

The range data accuracy from a precision radar is di- 
rectly related to variations in reply-delay time of the bea- 
con. The capability of current radars suggests that delay 
variations should not exceed 0.006 psec. This is a severe 
requirement when it is related to the necessity for chain- 
radar tracking and the wide differences in radar signal 
levels encountered by the beacon. In this circumstance, ref- 
erence is again made to the fact that the usual type receiver 
AGC is undesirable. 


Target Acquisition Techniques 


Two techniques used at downrange stations to acquire 
the target at maximum range involve switching to narrow- 
band radar receiver operation and the use of wider trans- 
mitter pulses. The narrow bandwidth provides an improve- 
ment in radar receiver signal-to-noise ratio of approxi- 
mately 4 db, and the wider pulses increase the range at 
which the beacon may be triggered, i.e., with same radar 
PRF. 


CONCLUSIONS 


This discussion has presented a brief review of several 
important and somewhat arbitrarily chosen problem areas 
concerning the design and utilization of radar beacons for 
IRBM and ICBM tests. It has been possible to use quan- 
titative data in only a few instances; first, because missile 
test program objectives differ; second, because of the un- 
certainties relating to the final outcome of current radar 
improvement programs. 

Portions of this discussion, however, are presented as a 
critique of current methods and techniques used in the 
formulation of beacon design criteria. This approach seems 
justified because of a continuation in basic design deficien- 
cies for a large number of the more recent beacon develop- 
mental projects. The primary difficulty seems to be based 
on failure of beacon designers to properly observe the 
“systems concept” in establishing beacon design require- 
ments. The designer is thus reminded that, first, and 
through whatever channel may be appropriate, he should 
expend all effort necessary to learn the characteristics of 
test range radars and the operational techniques utilized. 
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Design Studies for a Rubidium Gas Cell Frequency Standard* 
J. M. ANDRES}, D. J. FARMER}, anv G. T. INOUYE} 


Summary—A description is given of several studies under- . 


taken in the design of a small, gas cell stabilized atomic fre- 
quency standard making use of the field-independent hyperfine 
resonance of rubidium 87 at 6834 mc. In the standard considered, 
light from a rubidium vapor lamp is used to enhance the popu- 
lation differences between the atomic levels involved in the 
resonance and also to provide a means of detecting the micro- 
wave resonance. One of the studies described is concerned with 
the generation of requisite microwave energy at the resonant 
frequency, another with the choice of an optimum length for the 
gas cell, and the last with the optimization of parameters for 
the modulation process used. 


I. INTRODUCTION 


E present here the results of several studies under- 
\ \ taken in the design of an atomic frequency stand- 
ard. The standard is of the gas cell absorption type 
and is intended for possible use in missiles or space ve- 
hicles—applications which subject the device to stringent 
limitations on weight and power consumption, and which 
demand high performance for operation under vibration 
environments varying from moderate to severe. Of the 
many problems these and other requirements present to the 
designer, we shall consider here only three—one deals with 
generation of the requisite microwave energy, another 
leads to an optimum choice for the length of the absorp- 
tion cell, and the last is concerned with optimization of 
parameters for the modulation process employed. All three 
of these problems can be discussed within the framework 
of the functional diagram for the standard shown in Fig. 1 
and the brief description of its operation given below.* 
Optical resonance radiation from a rubidium lamp is 
transmitted through a glass cell containing a vapor of the 
isotope rubidium 87 and a small amount of argon or other 
buffer gas, the purpose of the latter being to reduce 
Doppler broadening in the microwave resonance of inter- 
est [6]. The action of the optical resonance radiation pro- 
duces a difference in population between the hyperfine 
levels of the ground state of Rb*’ far in excess of that 
which obtains under thermal equilibrium. This is the proc- 
ess of optical pumping [7], [8]. The energy separation of 
the hyperfine levels is approximately 6834 mc; when the 
atoms in the cell are subjected to a radiation field of that 
frequency, transitions between the hyperfine levels are in- 
duced which tend to reduce the population difference built 
up by optical pumping and, in so doing, to change the 
amount of optical resonance radiation transmitted through 
the gas cell to the optical detector. 


* Manuscript received by the PGMIL, July 21, 1959. 

+ Space Technology Labs., Inc., Los Angeles 45, Calif. 

More detailed information concerning the principles governing 
the operation of devices of this type is found in [1]-[5]. 
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Fig. 1—Functional diagram of frequency standard making use of 
the optical detection of the field-independent hyperfine resonance 
of rubidium 87 occurring at 6834 mc. Arrangements for pro- 
ducing the static magnetic field required and other control equip- 
ment are not shown. 


Modulating the phase or frequency of the microwave 
field a small amount at some low rate modulates the inten- 
sity of the light reaching the detector at the same rate. 
The amplitude of the modulated light intensity is in some 
sense proportional to the slope of the microwave resonance 
absorption line. By synchronously demodulating the am- 
plified output of the photosensor a de error voltage is ob- 
tained which is used to control the frequency of a piezo- 
electric quartz crystal oscillator. The oscillator operates at 
a subharmonic of the hyperfine resonance frequency and, 
through the multiplier chain, provides the microwave en- 
ergy used in irradiating the gas cell. If the control loop so 
formed is not closed and the microwave frequency is 
varied through resonance the error voltage developed at 
the output of the phase detector is as shown in Fig, 2. 


II. Microwave GENERATION 


Microwave power absorbed by the rubidium atoms in 
the gas cell is negligibly small; however, to induce transi- 
tions an RF magnetic field having an amplitude of about 
3 X 10° gauss must be maintained. A field of that magni- 
tude may be obtained with a power dissipation of 10-8 
watts in a cavity having a loaded Q of about 2000—the 
latter being a value realizable for a cavity enclosing the gas 
cell. This rather low power requirement for the 6834-mc 
signal suggests that it be derived as a harmonic from a 
crystal diode. Because of weight and power considerations 
it is desirable that the subharmonic signal driving the diode 
be generated in a transistor stage; thus the question of in- 
terest is the order of multiplication that should be required 
of the diode. The answer hinges on the state of the art of 
the components. 

One of the most recent advances has been the avail- 
ability on the commercial market of the diffused-base silicon 
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Fig. 2—The error voltage at the output of the phase detector as the 
microwave frequency is slowly passed through the atomic reso- 
nance. This curve was obtained in the laboratory with apparatus 
similar to that shown in Fig. 1. An indication of the noise pres- 
ent is found by increasing the narrow-band amplifier gain when 
the microwave frequency is off resonance. 


p-n junction diodes. The properties of one such device 
used as a frequency multiplier are shown in Fig. 3. We 
note that 40 milliwatts input at the 114th subharmonic, 
about 60 me, produces the required power. By comparison, 
if a diode of the usual type, such as a 1N23B, were to be 
used, measurements show that 500 milliwatts input would 
be required at the 50th subharmonic, or 40 milliwatts at 
the 17th subharmonic. 

Transistors presently available could be used in con- 
junction with the usual crystal diode ; however, use of the 
p-m junction diodes requires fewer transistor multiplier 
stages, and the lower frequencies involved result in a 
greater electrical efficiency. We can expect that improve- 
ments in both transistors and diodes may soon admit even 
further reduction in the number of components required 
in the microwave generator portion of the standard.’ 


III. Optimum CeLL LENGTH 


One example of system optimization concerns the choice 
of an optimum gas cell length. The optimum length is de- 
fined to be the length giving the maximum change in trans- 
mitted light intensity as the microwave frequency passes 
through the atomic hyperfine resonance frequency. The 
light absorbing character of the cell depends upon several 
factors: the length of the cell, the number of absorbing 
atoms present in the cell in a unit length, the distribution 
of these atoms among the several absorbing substates, the 
spectral character and polarization of the light incident 
upon the cell, and the atomic thermal relaxation effects in 
the cell. The distribution of atoms in the substates and, 
hence, the light absorbed by the cell is changed as the micro- 
wave frequency passes through resonance. The magnitude 
of change depends upon the cell length and the other factors. 


2Tt may be argued that further simplification of the design of 
Fig. 1 can be achieved even now by using a crystal oscillator of 
higher frequency. The choice of 5 mc is based on the apparent su- 


periority of available crystal units in that frequency neighborhood. 
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Big. J-—The frequency multiplication characteristics of a commer- 
cially available silicon diode (Microwave Associates 460E). 


Thus, in order to determine an optimum cell length, it is 
necessary to determine the dependence of the light absorp- 
tion by the cell upon these various factors. Let us first make 
some simplifying assumptions. 

We shall assume that the incident optical radiation is 
unpolarized and is of two frequency components which 
have the intensity g1 and g2 and which are capable of being 
absorbed by the F = 1 and F = 2 sublevels of the 5S1/2 
ground state of Rb*, respectively. When light is absorbed 
and the atoms are subsequently excited to the 5P3/2 or 
5P1/2 states, we assume that the states are completely mixed 
by collisions with the buffer gas atoms before reradiation 
so that the probability of atoms returning to any one of 
the eight sublevels is equal to one eighth of the total origi- 
nal absorption probability. In addition, we assume that the 
eight sublevels are mixed by means of thermal relaxation 
processes characterized by a relaxation time T, that re- 
radiated light cannot be further absorbed, and that the 
atoms at a given position in the cell remain relatively fixed 
in position during the absorption and reradiation process. 
With these assumptions it is possible to compute the light 
absorbed as a function of gas cell length and other factors 
and to determine the cell length giving the maximum 
change in transmitted light when the microwave frequency 
passes through the atomic resonance. 

For light of frequency v incident upon an atom, the 
probability of absorption per unit time per atom is 


P = mrocd(v)f (1) 


where rp = 2.818 X 10-2* cm, the electron radius, c = 3 X 
10*° cm/sec, the velocity of light, g (y), is the spectral in- 
tensity, the number of photons/cm? sec in a l-cps band, 
and f is a dimensionless number called the oscillator strength 
of the transition. 

The number of atoms per cm® in the three substates of 
F=1 capable of absorbing light is N(a: + a2 + as), 
where N is the total number of rubidium atoms per cm*® 
and ai, a2, and as represent the fraction of the total number 
in each of these substates. The atoms are not “tuned” to 


® For information about spectroscopic notation and atomic energy 
levels see [9]. [5] contains a figure showing the energy levels in 
sodium, which is similar to rubidium. 
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absorb a single frequency, but because of Doppler effect, 
are distributed to absorb over a band of frequencies of 
width App. If they are assumed to be uniformly distributed 
over this band, the number capable of absorbing light in a 
1-cps band at a frequency » within the band is 


N(a + (0D) aa a3) 
Avp 


If the intensity gi illuminating the F = 1 substates is 
assumed uniform over the band Avp, we have, for the 
change in gi over this band per length dy, 

d§,  — mrocfSi:N(ai + a + as) 
= 6 (2) 
dx Avp 


Similarly, assuming the intensity g. illuminating the five 
F = 2 substates is uniform over a band Ayp, we have, for 
the change g2 over the band App, 
dS — TrocfIoN (a4 + ds + ae + A7 + ds) 3) 
dx Avp 
The a; are functions of x through gi, g2 and also of the 
relaxation time 7. Let us consider what a;’s are for two 
cases of interest. 


Case 1—Microwave Frequency Off Resonance 


Under these circumstances a, = a2 = az and a = a; = 
dg = a; = as. By balancing the relative number of atoms 
leaving a state with the number re-entering the state, a1 
and a, can be determined. Thus, the processes of absorp- 
tion, reradiation, and relaxation must yield no net flow of 
atoms from a given state under equilibrium conditions. A 
determination of a, and a, gives: 


Ig ta 
gq = (4) 
5], + 31s a 8a 
Iit+e 
a (5) 
57, + 3le + 8a 
where 
1 g g 
Fe ae ea Gan fare 
trofcT Io Io Io 


(a, [1, I, are dimensionless.) 
Putting (4) and (5) into (2) and (3) we obtain 
dl, 312 + a) 


Ut Sip ee (9) 
dl, nee 5Io(11 + @) (7) 
dy 5h + 3le + 8a 
where 
mrocN fx 
ia Avp 


is a dimensionless quantity. 


IRE TRANSACTIONS ON MILITARY ELECTRONICS 


October 


For infinite relaxation time a0 and (6) and (7) re- 
duce to equations of the form obtained in [5]. 


Case 2—Microwave Frequency On Resonance 

We assume the microwave power saturates the transi- 
tion between substates 2 and 6. Thus a2 = ag. In addition 
@, = a3 and a4 = ds = a; = ag. Under equilibrium condi- 
tions we obtain: 


A(T, a a)(y a a) 


= S 8 
a, = a s (8) 
I, + Ip + 2a) (I. 
Figageemer ects 2 + 2a)(I2 + a) (9) 
A 
2[L2 + hls + alle + 3h) + 2a?] 


(10) 


Qy = 4 = 43 = a3 = 
A 


where 
A = 4[27,? + 5k, + Ip? + a(7Ie + 9l) + 802]. 
Putting these quantities into (2) and (3) we obtain 
dy — 1, [2hIg + 1? + (41, + 2h)a + 3a?] 


ai eee =(11 
dy [272+ 5hhlo+ 12 +a(7I.+ 91) + 8a?] ab) 
dl 


— 1,[20,? + 3hile + 3h + Th)a + 5a?] 
dy = [27,2 + 5h + I? + a7, + 9) + 8a?] 


(12) 


where y, a, /,, and J, are the dimensionless quantities pre- 
viously defined. 

The magnitude of the signal is defined by the change 
in total light intensity J, + J, as the microwave frequency 
is changed from off resonance to on resonance, i.e., the 
difference in the two cases just treated. The quantity S= 
(1, + Iz) ott res — (11 + Iz) res is plotted in Figs. 4-6 for 
three values of [1p and J20, the magnitude of J, and J, at 
y = 0, and for a number of values of a. 

The case of J1o = 1.0 and Jz) = 0.25 corresponds to il- 
luminating the F = 1 levels more intensely than the F = 2 
levels. This is achievable, in practice, by filtering the light 
from a Rb* lamp by passing it through a cell containing 
Rb*® [10]. The case I1o = J2o = 1.00 corresponds roughly 
to illuminating with a lamp containing naturally occurring 
rubidium. The case Jio = 0.25, I20 =1.0, corresponds to 
illuminating with light from a Rb** lamp. Some reasonable 
values for the quantities of interest are: 


$9 = 10° photons/cm? sec sec}, 
T = 10- sec, : 
fmt, 
Avp = 5 X 108 cps, 
N = 10" atoms/cm; (temperature~40°C). 
These values give a = 0.04 and 
a(cm) 


= = 4 for « = 3 inches) 
y Te (y or x = 3 inches). 


Fig. 4—The signal quantity, S, is a function of the dimensionless 
length y, for several values of the parameter a. Dimensionless 
intensities at y= 0 are J = 1.0, J = 1.0. 


Fig. 5—The signal quantity, S, as a function of the dimensionless 
length y, for several values of the parameter «. Dimensionless 
intensities at y= 0 are Ix. = 1.00, Io = 0.25. 
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Fig. 6—The signal quantity, S, as a function of the dimensionless 
length y, for several values of the partmeter a. Dimensionless 
intensities at y = 0 are 1 = 0.25, Iz = 1.00. 
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Looking at the curves fora = 0.05, it is seen that a 
maximum occurs somewhere between y = 6 and y = 20, 
depending on the nature of the illumination. We have 
chosen to make the cell length approximately six inches, 
using a cavity mode TEo14. As the temperature is increased, 
thus increasing N and y, the signal obtained in the lab- 
oratory exhibits a maximum at a temperature of about 
37°C. Thus, the fact that such a maximum exists is readily 
demonstrable. Since the intensities g1 and gz emitted by 
the lamp are not known with any precision, it is difficult 
to determine a definite identification between these theoreti- 
cal curves and the laboratory results. The curves indicate 
that operation with a Rb® lamp should give a larger signal 
for a cell which is shorter than required with the other 
lamps—a fact in agreement with some tentative laboratory 
results. 


IV. OPTIMIZATION oF MoDULATION PARAMETERS 


It is convenient to consider the role played by the atomic 
system as that of a frequency discriminator of very high 
stability. Departures of the controlled crystal oscillator 
from the center frequency of this discriminator are sensed, 
and the error voltage is used to correct the frequency of 
the oscillator. The process is thus analogous to the opera- 
tion of usual automatic frequency control systems, and a 
similar analysis of the control loop may be applied. An 
important parameter of such a loop is, of course, the slope 
of the discriminator characteristic at center frequency. 
Other things being the same, it is desirable to obtain as 
great a value for the slope as circumstances will allow. eee 
previous section dealt with one means of increasing the 
resonance signal and, hence, the discriminator slope. The 
present discussion considers the effect of the modulation 
parameters on the discriminator slope and the method by 
which a careful choice of these parameters can increase 
the slope. 

As a starting point we use the result of Karplus [11] for 
the signal produced when an atomic resonance is observed 
with a frequency-modulated microwave excitation. The 
general result relates the modulation rate Q, the frequency 
excursion amplitude of the modulated wave &, the func- 
tional form of the frequency modulation h(Qt), and the 
deviation of the unmodulated microwave frequency from 
the line center &, toa function F containing the time depend- 
ence of the microwave absorption coefficient. 


F =Im {i i dr exp E a dt'w'h(Qt’) 
| — (t@ + ovr = 1 (13) 


The quantity D is the reciprocal of the relaxation time T; 
it will depend on processes occurring within the gas cell 
but not on the modulating parameters and is assumed as 
given. 
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Provided the probability per second for optical photon 
absorption in the pumping process is greater than Q, i.e., 
1—1<,,,, intensity modulation of the optical reso- 
nance radiation transmitted through the gas cell to the de- 
tector will be proportional to F. For these results to be 
applicable, the microwave field must be sufficiently weak 
so that no “saturation” effects occur. 

If the microwave radiation is phase modulated at a rate 
Q witha deviation of amplitude ¢ such that its instantaneous 
phase is wt -+ ¢ sin Qt, then we may write ® = $Q, and for 
the frequency modulation function, h(Qt) = cos gt. With 
this form for h, (13) may be integrated to yield for the 
component of F oscillating at frequency Q, 


1 
Fo = es [A cos Qt — B sin al (14) 
with 
— Ie($) [Jou($) + Joa($)] 
iy a (i +O? 41 ) 
= J.($) [Jsr1(¢) — Js(¢)] : 
B= = cen ian -(sl +c), (15b) 


where J;(g) is the Bessel function of the first kind of or- 
ders, 


and 


While components of F at frequenciesnQ (” = 2,3,...) 
are also present, only the fundamental is usually of interest 
since in practical systems the signal from the photosensor is 
processed in a narrow-band amplifier tuned to that fre- 
quency. 

Note that F° has two components, one of amplitude 4 
oscillating in phase with the frequency modulation func- 
tion h (but out of phase with the phase modulation), and 
a second of amplitude B, out of phase with h. It may be 
shown that B > 0 as either b >0 or Q > 0 whereas A is 
nonzero for b = 0 or Q —0. Thus, the presence of the 
out-of-phase component arises from relaxation processes, 
as we might expect. Both A and B vanish as # — 0 or as 
@ > 0. 

We may view the operation of the phase detector as pro- 
ducing the time average of the product of its inputs. One 
input is F g and the other is the phase modulation function ; 
to account for phase delays in the narrow-band amplifier 
or for delays that may be intentionally introduced the ref- 
erence input is represented as proportional to sin( Qt + 
). The phase detector output is thus proportional to 


A sin 6 — B cos ® (16) 


which describes the general shape of the discriminator 
characteristic for the system. The curve shown in Fig. 2 is 
an example of the function described by (16). 
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In seeking the maximum slope at the center frequency 
of the characteristic, i.e., at c or @ = 0, let us ignore for 
the moment the latitude allowed by choice of © and exam- 


ine 
( ~) ( ~) 
— and es : 
OC c=0 0c c=0 


Those two quantities have been calculated for various values 
of J and g and are displayed in Figs. 7 and 8. The values of 
the modulation parameters yielding maximum slope in 
either case may be calculated ; we find 


Cor 
max 0c c=0 . 


OB 
(=) = 0.750, at = 1.732, g = 1.836. 
max 0c c=0 


— 0.770, at? = 0, g = 0.707, 


and 


Examining Figs. 7 and 8, it is noted that for values of / 
and g which result in large slope for one component, the 
slope of the other components is rather small; we conclude 
that ® ~ 7/2 or © 0 (the former giving A as the dis- 
criminator characteristic and the latter, B) will yield con- 
ditions of maximum slope provided the appropriate values - 
of J and g are chosen. On the basis of optimum slope 
alone, use of the A characteristic is clearly indicated ; how- 
ever, other factors may have to be considered. For several 
reasons, chiefly to reduce 1/f type noise generated in the 
optical pumping lamp or the photosensor and subsequent 
amplifier, it is desirable to use as large a modulation rate 
as practicable. Attaining its maximum slope at a modula- 
tion rate different from zero, the B component accommo- 
dates this condition very conveniently. 

Curves similar to those shown in Figs. 7 and 8 have also 
been calculated for frequency modulation whose form is a 
square wave. The results are similar, except that the maxi- 
mum slope of the 4 characteristic is somewhat greater at 
zero modulation rate, and the slope for a given value of g, 
in general, falls off more rapidly with increasing modula- 
tion frequency. For the B component, the maximum slope 
is about 0.5. We have not undertaken to solve for the op- 
timum form of the modulation function, nor do we expect 
that the effort would be very rewarding since there is not 
a very great difference in results between the two rather 
extreme cases already discussed ; furthermore, the compli- 
cation in producing a modulating function different from 
sinusoidal would probably weigh against its use in any 
practical design. 

The remarks which have been made regarding the de- 
sirability of using the B characteristic are based on the 
original assumption that other things are the same. The 
other things, of course, are the noise contributions of vari- 
ous portions of the system. Noise which may be viewed as 
amplitude modulation, such as produced by fluctuations in 
the lamp intensity or the usual noise from the detector and 
following circuits, will effect the output of the phase de- 
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Fig. 7—The slope at resonance of the component of the optical in- 
tensity modulation out of phase with the frequency modulation 
function. 


Fig. 8—The slope at resonance of the component of the optical in- 
tensity modulation in phase with the frequency modulation 
function. 
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tector similarly whether the A or B component is used. 
Another type of noise, spurious phase modulation compo- 
nents on the microwave signal, such as may be produced 
either in the crystal oscillator itself or in the subsequent 
circuits, may affect the separate components differently. 
When the total phase deviation at the microwave frequency 
is less than a radian, it proves possible to obtain an analytic 
representation for the noise spectrum at the output of the 
phase detector; however, the results in the form that we 
presently have them do not admit of ready comparison for 
the two separate cases. Within the present state of the art, 
noise produced in the optical pumping lamp and the photo- 
sensor seem to be the dominant sources of frequency in- 
stability. Presumably these contributions can be decreased 
as technology progresses, and the phase type noise may 
then become more important. 


V. ACKNOWLEDGMENT 


It is a pleasure to acknowledge several illuminating dis- 
cussions with Dr. R. C. Ward regarding the modulation 
parameter optimization, and the assistance of H. Tuens 
and Miss R. Gittelman of our Computation Center who 
programmed the equations of Sections III and IV for the 
IBM 704 computer. 


VI. REFERENCES 


[1] J. P. Wittke, “Molecular amplification and generation of 
microwaves,” Proc. IRE, vol. 45, pp. 291-316; March, 1957. 

[2] R. H. Dicke, T. R. Carver, and C. O. Alley, “Investigations 
of the Effect of Gas Collisions and Optical Pumping on the 
Breadth of Spectral Lines,’ first through twelfth quarterly 
reports by Palmer Physical Lab., Princeton University, 
Princeton, N.J., to U.S. Army Signal Eng. Lab., Fort Mon- 
mouth, N.J., on Contract DA-36-039 SC-70147 ; 1955-1958. 

[3] W. E. Bell and A. Bloom, “Optical detection of magnetic res- 
onance in alkali metal vapor,” Phys. Rev., vol. 107, pp. 1559- 
1565; September 15, 1957. 

[4] W. E. Bell and A. Bloom, “Optically detected field-independ- 
ent transition in sodium vapor,” Phys. Rev., vol. 109, pp. 219- 
220; January 1, 1958. 

[5] W. E. Bell, A. Bloom, and R. Williams, ““A microwave fre- 
quency standard employing optically pumped sodium vapor,” 
IRE Trans. on Microwave THEORY AND TECHNIQUES, vol. 
MTT-, pp. 95-98 ; January, 1959. 

[6] R. H. Dicke, “The effect of collisions upon the Doppler width 
hese lines,’ Phys. Rev., vol. 89, pp. 472-473; January 15, 
1953. 

[7] A. Kastler, “Optical methods of atomic orientation and of 
magnetic resonance,” J. Opt. Soc. Amer., vol. 47, pp. 460-465 ; 
June, 1957. 

[8] H. G. Dehmelt, “Slow spin relaxation of optically polarized 
ee atoms,” Phys. Rev., vol. 105, pp. 1487-1549; March 1, 
1957. 

[9] H. E. White, “Introduction to Atomic Spectra,” McGraw- 
Hill Book Co., Inc., New York, N.Y.; 1934. 

[10] P. L. Bender, E. C. Beaty, and A. R. Chi, “Optical detection 
of narrow Rb™ hyperfine absorption lines,” Phys. Rev. Letters, 
vol. 1, pp. 311-313 ; November 1, 1958. 

[11] R. Karplus, “Frequency modulation in microwave spectros- 
copy,’ Phys. Rev., vol. 73, pp. 1027-1034; May 1, 1948. 


184 


IRE TRANSACTIONS ON MILITARY ELECTRONICS 


October 


Some Considerations in the Design of the Guidance 
and Control System for Discoverer* 


JOHN J. SCHMITT, Jr.+ 


Summary—The limited number of flight tests to date has in- 
dicated that the Discoverer satellite will be a highly reliable and 
effective system from which the ARPA and the Air Force expect 
to gain valuable additional knowledge of outer space. The 
achievement of this high degree of reliability and effectiveness 
of the satellite means that each of its subsystems must, in fact, 
be even more reliable and effective. The Discoverer guidance 
and control system is meeting their specifications because of the 
simplicity of the over-all design—simplicity for reliability with 
sufficient sophistication to accomplish its mission effectively. 

The mission of the guidance system is that of providing ade- 
quate attitude, velocity, and time references to guide the satellite 
into a desired orbit, to control its attitude on orbit, and to pro- 
vide switching in the proper sequence for certain principal satel- 
lite functions. The mission of the control system is simply that 
of positioning the vehicle in the proper attitude commanded by 
guidance. The investigations that were made and those which 
are being continued in the determination of the optimum con- 
figuration for the desired performance are described in some 
detail. The present system, is discussed, which comprises, for 
guidance, an inertial reference package, a horizon scanner, and 
a computer for grammar, and for control, both a hydraulic and 
pneumatic system. 

In summary, a typical flight of the Discoverer, from separation 
from the Thor booster until stabilization on orbit, is outlined and 
the functions of the guidance and control systems are described. 


HE limited number of flight tests to date has indi- 

cated that the Discoverer satellite will be a highly 
reliable and effective system from which the ARPA 

and the Air Force expect to gain valuable additional in- 
formation and knowledge of outer space. The achievement 
of this high degree of reliability and effectiveness of the 
satellite means that each of its subsystems must, in fact, 
be even more reliable and effective. The Discoverer 
guidance and control system is meeting its specifications 
because of the simplicity of the over-all design—simplicity 
for reliability with sufficient sophistication to accomplish 
its mission effectively. It is the purpose of this article to 
describe the consideration given to this philosophy, sim- 
plicity with sophistication, in the design of the system. 
First, the mission of the guidance and control system 
must be defined. The mission of guidance is that of pro- 
viding adequate attitude, velocity, and time references to 
guide the satellite into a desired orbit, to control its attitude 
on orbit, and to provide switching in the proper sequence 
for certain principal satellite functions. The mission of 
control is simply that of positioning the vehicle in the 
proper attitude commanded by guidance. Although many 
involved and complicated techniques can be conceived to 
perform the missions defined, the approach in this case was 


* Manuscript received by the PGMIL, July 21, 1959. 
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to design a system which utilized simple, well-known, 
proven techniques as often as possible. 

As previously defined, the function of guidance is to 
provide attitude, velocity, and time references. Consider 
these requirements separately. Many designs, including 
various inertial platforms, were considered for attitude 
determination. With some of the platform considered, it 
would have been necessary to use an autopilot which in- 
cluded a body-mounted inertial reference system. The logi- 
cal next step in this investigation was to determine if a 
body-mounted gyro-reference system alone would be ade- 
quate. This would negate the high cost and added com- 
plexity of an inertial platform. In place of a detailed dis- 
cussion of satellite trajectories, it will be sufficient to state 
that the vehicle attitude at injection into orbit is extremely 
critical in achieving a desired orbit. An analysis indicated 
that the typical autopilot reference gyros alone would not, 
under the expected environment, provide an adequate ref- 
erence. Consideration was given to the possible use of a 
horizon-scanning device which would give an adequate 
reference and which could accordingly correct the autopilot 
gyros. A number of studies had been made which in- 
dicated that such a device was not only feasible, but that 
its development could be accomplished quickly and at nomi- 
nal cost. The combination body-mounted, inertial-refer- 
ence, package-horizon scanner appeared to be very promis- 
ing, especially since the scanner had real potential for use 
on orbit as an attitude-indication and control-system refer- 
ence device. 

The velocity indication requirement was met simply with 
a longitudinal accelerometer whose output drives a motor- 
driven integrator. 

Although more sophisticated and elaborate timing de- 
vices were investigated, an accurate electromechanical 
timer was chosen as a time reference. This timer was de- 
signed to include a number of cam-operated switches for 
the proper programming of certain satellite functions. If 
necessary, the cams are easily reset to accommodate various 
program changes from flight to flight. ~ 

The selection of a control system capable of accepting 
the guidance system commands and of accordingly posi- 
tioning the satellite remained. The necessary servo loops 
which included the common use of rate gyros for stability 
purposes were designed. Because the attitude of the satel- 
lite must be controlled not only during the time when the 
orbital injection rocket engine is operating, but also before 
and after that period, two sets of “muscles” were neces- 
sary. Hydraulic actuators were chosen for the engine gim- 
baling in pitch and yaw. This hydraulic control system uses 
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an electrically driven pump to supply the necessary pres- 
sure. To control the satellite in pitch and yaw at times when 
the rocket engine is not operating, and in roll at all times, 
an array of gas reaction jets was chosen. Consideration 
was given to simple on-off-type valves which supply gas to 
the appropriate jets but, because of the need for more 
accurate control at the minimum expenditure of stored 
gas, the development of proportional valves which supply 
gas to the appropriate jets in direct relation to the size of 
the required correction was undertaken. The choice of this 
pneumatic control system completed the determination of 
the components for the guidance and control system. 
Throughout the entire determination, the proper balance 
between simplicity and sophistication was uppermost in the 
minds of the designers. 

To complete the description, it must be shown how these 
components are integrated into a system to perform the 
job of guiding the Discoverer into a desired orbit and 
of controlling its attitude on orbit. Upon separation from 
the Thor booster, the satellite coasts to a point just prior 
to the apogee of its ascent ellipse. The orbital injection 
rocket engine is started at this point and adds sufficient 
horizontal velocity to the vehicle to cause it to orbit. When 
the proper orbital velocity has been attained, the engine is 
shut down. The satellite, now on orbit, maintains a local 
horizontal attitude. Just prior to separation from the 
Thor, the inertial-reference package-position gyros are 
uncaged in order that they may use as a reference the 
attitude of the Thor/Discoverer at Thor burnout. Upon 
separation, the timer initiates a torquing program to the 
pitch gyro which causes the satellite, by operation of the gas 
jets, to pitch to the local horizontal. The timer then actuates 
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the horizon scanner which maintains control of the vehicle 
in the local horizontal position through the inertial refer- 
ence package and pneumatic control loop. At the proper 
time, the timer initiates the rocket-engine start sequence 
and switches control of the vehicle in pitch and yaw from 
the pneumatic system to the hydraulic system. The horizon 
scanner and intertial reference package continue to main- 
tain control of the vehicle in the local horizontal position. 
During the orbital injection the longitudinal accelerometer 
and integrator operate to measure the velocity gained dur- 
ing this period. When the proper velocity has been gained, 
the integrator causes shutdown of the satellite engine. 
Once on orbit, the timer switches control of the vehicle in 
pitch and yaw back to the pneumatic control system. Hori- 
zontal attitude control on orbit, an important first in the 
increased performance of control systems, is thus main- 
tained in the same manner as it was during the earlier coast 
phase of the vehicle’s flight. 

The ARPA, Air Force, and Lockheed (the prime con- 
tractor for the Discoverer), as well as those Lockheed sub- 
contractors involved, are proud of our accomplishments. 
We are not, however, complacent. Components can be im- 
proved and simplified, power requirements can be reduced, 
as can be the number of switching functions. The latter 
are some areas in which we expect to gain more reliability 
and effectiveness. 

Today more stringent and demanding guidance and con- 
trol requirements challenge the Air Force and all associ- 
ated industry. These requirements can and will be met by 
utilizing the knowledge gained from the Discoverer and 
other early satellite developments, and by applying the basic 
design philosophy—simplicity with sophistication. 
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laboratory of the Air Force missile develop- 
ment center, Holloman AFB, N.M., where 
as a task scientist in the operations section, 
he is currently doing radiobiological work 
in the bioastronautics branch which is con- 
nected with biomedical programs designed 
for research into human space flight. 

Dr. Hoffman is a member of Phi Zeta, 
Phi Kappa Phi, the American Veterinary 
Medical Association, and the American Col- 
lege of Veterinary Toxicologists. 
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George T. Inouye (S’46-A’53-M’58) was 
born on November 27, 1922, in Sacramento, 
Calif. He received the B.S. degree in electri- 
cal engineering from 
Swarthmore College, 
Swarthmore, Pa., in 
1944 and the Ph.D. 


degree in applied 
physics from Har- 
vard University, 


Cambridge, Mass., in 
1953. 

He was employed 
by Technicolor Cor- 
poration, Hollywood, 
Calif., from 1953 to 
1957. Then, in 1957, 
he joined the Ramo-Wooldridge Corpora- 
tion guided missile research division (later 
the Space Technology Laboratories) where 
he is presently employed. 

Dr. Inouye is a member of Sigma Xi. 
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Aubrey B. Mickelwait was born in Hono- 
lulu, Hawaii, on May 2, 1925. He received 
the B.S. M.S., and Ph.D. degrees in 
physics, from Car- 
negie Institute of 
Technology, _ Pitts- 
burgh, Pa. 

He has had six 
years’ experience in 
the quantum field 
theory, theoretical 
physics, meson the- 
ory and is the author 
of several technical 


papers. He joined 
A. B. MickELwair Space Technology 
Laboratories, Los 


Angeles, Calif., in 1955 as a member of the 
technical staff and is presently head of the 
Navigation Section, Guidance and Nav- 
igation Department of the Guidance Re- 
search Laboratory. There he has been en- 
gaged in research and systems analysis 
work on the Atlas, Titan, Thor, and 
Minuteman guidance systems. During the 
last several years, he has done theoretical 
and design work on trajectories, payloads, 
and guidance systems for space vehicles 
such as the Pioneer moon flights. 

Dr. Mickelwait is a member of the 
American Physical Society, American Geo- 
physical Union, Sigma Xi, Institute of 
Aeronautical Sciences, and Tau Beta Pi. 
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Robert A. Park was born in Stockholm, 
Sweden on September 4, 1923. He received 
the M.A. degree in English in 1949 and the 
Doctor of Laws de- 
gree in 1955, both 
from the University 
of Chicago, Chicago, 
Til. 

From 1949 to 1950, 
he taught creative 
writing at the Uni- 
versity of Chicago; 
during 1949, he 
taught English at the 
University of  Illi- 
nois, Urbana; from 
1949 to 1955, he 
taught English, humanities, and reading at 
Wilson Junior College, Chicago, Ill. From 
1955 to 1957, he was an Assistant Attorney 
General for the State of Oregon. He joined 
Space Technology Laboratories, Los An- 
geles, Calif., in July, 1957, and is presently 
connected with a special project on the Ex- 
ecutive Vice-President’s staff. 

Mr. Park is a member of the Oregon 
Bar Association and the California Bar 
Association. 
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Major John J. Schmitt, Jr., graduated 
from West Point Military Academy in 1946 
and in 1951 attended the USAF Institute 
of Technology at 
Wright - Patterson 
Air Force Base, Day- 
ton, Ohio, where he 


attended a _ fifteen- 
month course in 
graduate armament 


engineering. He con- 
tinued his education 
at the Massachusetts 


Institute of Tech- 
nology, Cambridge, 
J.J. ScumittT — where he _ studied 


weapon systems en- 
gineering and received the M.S. degree in 
1954. 

In 1946, he was assigned to the Strategic 
Air Command where his primary duty was 
that of pilot and aircraft commander on 
B-29’s and B-50’s. In 1954, he was arming 
and fuzing project officer at the Air Force 
Special Weapons Center, Kirtland AF 
Base, N.M. He was assigned to the Air 
Force Ballistic Missile Division, Los An- 
geles, Calif., early in 1958, to his present 
duty as guidance and control project officer 
for Discoverer. 
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John E. Taber was born on May 23, 1925, 
in Burbank, Calif. He received the B.S. 
degree in electrical engineering in 1946 
from the California Institute of Technol- 
ogy, Pasadena, and the M.S. degree in 
electrical engineering from the University 


Contributors 


of California, Berkeley, in 1951. 

In 1946, Mr. Taber joined the Hughes 
Aircraft Company, Los Angeles, Calif., 
where he worked on 
instrumentation and 
telemetry equipment. 
In 1948 he became 
employed by the 
Glenn L. Martin 
Company, Baltimore, 
Md., in the telem- 
etry section. Upon 
returning to Hughes 
Aircraft Company in 
1950, he did develop- 
ment work on high- 
speed digital com- 
munication systems. He continued this type 
of work at the communication division of 
the Ramo-Wooldridge Corporation, Los 
Angeles, from 1954 to 1958. He joined 
Space Technology Laboratories, Inc., Los 
Angeles, in 1958 and is presently associate 
manager of the instrumentation depart- 
ment there. Mr. Taber holds three patents 
concerning digital communication tech- 
niques. 


J. E. Taser 


Edwin H. Tompkins, Jr.. (S’48-A’50- 
M’56) was born in Charlottesville, Va. on 
October 9, 1928. He received the B.S. de- 
gree in electrical en- 
gineering in 1948 
from Virginia Poly- 
technic Institute, 
Blacksburg, Va. and 
the M.S. and Ph.D. 
degrees in 1952 and 
1955, respectively 
from North Caro- 
lina State College, 
Raleigh. 

From 1951 to 1957, 
E. H. ToMPKINS fe was associated 

with the Rocket Re- 
search Group at North Carolina State 
College, working in exterior ballistics of 
rocket flight. During 1957 and 1958, he 
was assistant professor of mathematics 
and electrical engineering at the college 
and has had experience in electrical cir- 
cuit theory and nonlinear mechanics. He 
joined Space Technology Laboratories, Los 
Angeles, Calif., in 1958 and is presently a 
member of the technical staff, Guidance 
and Navigation Department of the Guid- 
ance Research Laboratory, Research and 
Development Division. 

Dr. Tompkins is a member of the 
American Mathematical Society and Sigma 
Xi. 


Major Henry Triwush was born in Chi- 
cago, Ill., in 1919. He received the B.S. 
degree in chemistry from the University of 
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Chicago in 1940 and did further graduate 
work there. In 1946-1947, part of which 
time he had a Lederle fellowship, he at- 
tended Northwestern 
University in Evans- 
tons LS sander e= 
ceived the M.S. de- 
gree in physical 
chemistry in 1949. 
In 1942, he joined 
the Central Labora- 
tory of Division 10 


of the NDRC, lo- 
cated at Northwest- 
erm University, 


H. TrrwusH 


where he assisted in 
research on gas mask 
absorbents. He enlisted in the Army Air 
Corps in 1942 and received his commission 
as a navigator in 1944. After a combat tour 
with the Eighth Air Force in England, he 
returned to this country for duty as a radar- 
bombardier instructor until he reverted to 
inactive status in 1945. He was recalled to 
active duty in 1948 and has remained in the 
Air Force since that time. Major Triwush’s 
principal service assignments since recall 
have been with the Strategic Air Com- 
mand, with the Electronics Components 
Laboratory, Wright Air Development 
Center, Dayton, Ohio, and currently, after 
having attended the Air Force Command 
and Staff School, with the Air Force Bal- 
listic Missile Division, where his duties in- 
volve dealing with space system reliability 
problems. 

Major Triwush is an associate member of 
Sigma Xi and a member of Phi Lambda 
Upsilon. Presently, he holds the aeronautical 
rating of master navigator. 
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Paul F. von Handel was born in Bisenz, 
Austria, on January 30, 1901. He received 
the M.S. degree from the Technical Uni- 
versity, Munich, Ger- 
many, and later, the 
Ph.D. degree in elec- 
tronics from the 
Technical Univer- 
sity, Berlin, Ger- 
many, in 1931. 

He was formerly 
director of the Insti- 
tute for Electronics 
4 Research of DVL 
in Berlin, a member 
of the German Acad- 
emy of Sciences, 
president of the Board of Electronic Navi- 
gation, and professor at the Technical Uni- 
versity, Berlin. His special fields of activity 
were studies on wave propagation and elec- 
tronic navigation, for which he received 
the Lilienthal Prize for Aeronautical Re- 
search in 1944, 

He arrived in the U.S.A. in 1953, and has 
since been working on problems related with 
kinematics, guidance and control of guided 
missiles, wave propagation, and electronic 
warfare. 

Dr. von Handel is presently chief of the 
Orbital Mechanics Division, Directorate of 
Advanced Technology, Air Force Missile 
Development Center, Holloman, AFB, N.M. 
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